


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1971 


Dynamic plastic behavior of spherical shell. 


Nguyen-Tien-Ich 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/15673 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
_ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist iL Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


" \ KNOX appointed — and published -- scholarly author. 
| inp 
, LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 








DYNAMIC PLASTIC BEHAVIOR OF SPHERICAL SHELL 
by 


Nguyen-Tien-Ich 


B.A., Saigon University (Viet-Nam) 
(1962) 


S.M. and Naval Eng., Massachusetts Institute of Technology 
(1969) 


SUED SNS PAR ETAL FULFILLMENT 
Ot tibes hOUrRE MENTS SOR THE DEGREE soe 
DOC CR Ore EH LOSOray 
AT THE 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


MAR eek 





LIBRARY | ner 
NAVAL POSTGRADUAT!: ec HOOT 


MONTEREY , Ce 93940 


Dy Avie ee ASC BEAAVIOR OF SPHERICAL SHELL 


by 
Nguyé6én-Tién-Ich 


Submitted to the Department of Naval Architecture and Marine 
Emgineering on May 7, 1971, in partial fulfillment of the 
requirement for the degree of Doctor of Philosophy. 


ABSTRACT 


A dynamic pressure decaying exponentially with time, 
whose analytic representation in terms of non-dimensional 
@Uuantities is: . 
' ps poe : (t > 0) , 
and whose peak value py, is larger than the static collapse 
Pressure p., 1S applied uniformly throughout a spherical 
shell made of rigid-perfectly plastic material. The problem 
1s to determine the final deflection of the shell for a 
given yield surface, the shell has no holes and 1s simply 
Supported at the outer edge. 


Hie cm nox ave syleld Surtaces have been used. etme 
two-moment limited interaction (HS), the uncoupled square 
(SS) and the uncoupled diamond (DS) yield surfaces. 


For each yield surface selected, the static collapse 
pressure is first determined both for simply Supported and 
clamped edges. In the dynamic problem, the shell is simply 
supported. 


The first part deals with shallow shells. With the 
(HS) and (SS) yield surfaces, solutions for three ranges of 
meessure have been obtained. With the (DS) yield surface, 
solution for one range of pressure has been obtained. 


The second part deals with deeper spherical caps. For 
all three yield surfaces, only solutions for one range of 
pressure have been obtained. 


For shallow shells, graphs giving maximum central 
deflection and energy absorbed as functions of the pressure 
Pie ccomene- pp did Maximum central deflection as functions 
of energy absorbed are also presented. 
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Global and Local Coordinates of the Reference Surface 


Pegi ab 


@z> AxXis of nee 

M: Generic point on the reference surface, which is the 
mid-surface for a homogeneous shell of constant 
thickness. 


r: Distance from M to the axis of symmetry 0z. 


z: Distance from M to the tangent plane at the apex QO. 
m: Unit normal to the reference surface, positive inward. 
ey Unit tangent to the meridian at M, pointing away from 
the apex 0. 


x Unit tangent to the parallel at M. 


owt 

: naa ee ae 
lee Displacement component in the 7 dileeer LOM 
U_: Displacement component in the n direction. 





eal 


A: A reference length. 


To: RaeeBererencen cime. 

"6 

zx Non-dimensional meridional component of the 
ate elacement. 

On 

= Non-dimensional normal component of the displace- 
ment. 

t: Time variable. 

= Non-dimensional time variable. 

A . 

=? Non-dimensional space variable. 

~ Non-dimensional space variable. 

ae Non-dimensional meridional component of the 
VielOGr tye 

= , Non-dimensional normal component of the velocity. 

2g: Non-dimensional strain rate in the - direction. 

hf « r s e ng e * 

Cg: Non-dimensional strain rate in the ty Cie ec Cirem: 

Non-dimensional curvature rate in the Gibinie Gael Ome 

: +. 

Ko: Non-dimensional curvature rate in the C4 direction. 
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Local Stresses 


Oy: Membrane stress in the e direction. 
Vat Membrane stress in the > direction 


o ,: Transverse shear stress. 
2H: Shell thickness. 
E: Distance from a point of the shell to the reference 
surface, positively inward. 
For thin shell, for which — << 1, where R_ is the 
smaller principal vagina oF curvature of the reference 


Surface, the stress and moment resultants per unit 


JNO) co, ie © 
ariel ae lil ast 
se=fouogae » x= fo oyae . = [ey gae 
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gq,6dqdq , M | O,cde 
ae ea: : ao 
o.: Yield stress. 
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Non-dimensional 


—> 
t, direction. 


Non-dimensional 


a e e 
eee OL rec Gaon. 


Y 


Non-dimensional 
resultant. 
Non-dimensional 


o,-stresses. 


0 


Non-dimensional 


Oo,-stresses. 


Y 


Rate of internal energy 


Rate of external energy 


Non-dimensional rate of 


dissipation. 


Non-dimensional rate of 


sug yon ne e 


stress resultant in the 


stress resultant in the 


transverse shear stress 


moment resultant of the 


moment resultant of the 


dissipation. 


LTO 


internal energy 


external energy 


Surface density of the shell material. 


ING IDI! Soiree jg esses: j7 stig telae) etqvoneiieLll 


direction. 


Non-dimensional pressure applied exter- 


nally in the normal direction. 


Non-dimensional applied traction at the 


e a e ® 
edge, in ene n=-d2 rection. 


se 





» | 


ae 


Non-dimensional applied traction at the edge, 
in the Hg CEOS tu 


Non-dimensional applied moment at the edge. 


Non-dimensional surface density of the shell 


material. 
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wee er ODUC TION 


The first paper on the dynamic plastic behavior of a 
beam was published in the early 50's by E. H. Lee and 
Peo. Symonds {i]. In this problem and in those by sub- 
sequent authors [2, etc.], the material 1s assumed to be 
rigid perfectly plastic and the vibrations of the struc- 
tures are neglected so that the problem of unloading has 
not to be considered. These assumptions, which simplify 
problems considerably, do not introduce important errors 
in many cases, Since on the one hand dynamic energies are 
sometimes much larger than the elastic energy, and on the 
other hand the active duration of a pressure pulse is often 
men Smaller than the corresponding natural time period of 


mime Structure. 


Theoretical investigations oe the dynamic behavior 
of axisymmetric plates and shells were initiated with a 
paper published on circular plates by H. G. Hopkins and W. 
Prager [3]. Various boundary, geometrical and loading con- 
ditions have been considered: H. G. Hopkins and W. Prager 
solved the problem of simply ssoeetad circular plates 
igen inttoOrm dynamic load OF rectangular profile [3j; A: 
wu. Wang gave a solution of simply supported circular plates 
under impulsive loading [4]. Recently, M. F. Conroy pub- 


lished a paper on simply supported circular plates due to 





IVS 


dynamic circular loading [16]. The problem of clamped 
circular plates was treated by A. J. Wang saa H. G. Hopkins 
Pominioulsive load [5] and by A. L. Florence for a uni- 
form rectangular pulse load type [6]. A. L. Florence also 
solved the problem of clamped plates subjected to a 
rectangular pulse load type uniformly distributed over a 
central circular area [7]. The problem of an annular plate 
clamped along its inner edge and free along its outer one 
and subjected to impulsive load was examined by G. S. 


Smapiro [8]. 


In all these problems, the maximum displacement was 
Supposed to be infinitesimal so that the strain maee 
velocity relations ~— linear and only the bending moments 
are considered to be important. N. Jones examined the 
problem of finite deflection of simply supported circular 
plates under impulsive load and developed a theoretical 
procedure which retains both membrane forces and bending 
moments [9]. The solutions compare favorably with some 


experimental results on simply Supported circular plates. 


Although the behavior of circular plates has been con- 
Sidered in some detail, there are few solutions available 
for the axisymmetric behavior of shells. P. G. Hodge 
meucdted the wrigid plastic behavior of an infinitely long 
cylindrical shell reinforced with equally spaced rigid rings 


uniformly loaded with a rectangular pressure profile [10]. 





Ie 


Another paper by P. G. Hodge with a numerical solution, 
examined the influence of blast characteristics on the 


final deformation of circular evi indricalesiel ts. pli]. 


egos eee ne em CWONCeNOr cd gyre ldesurraee 15 Oc 
crucial importance. In most of the problems dealing with 
plates, the Tresca yield surface, peeaner Of LesmeeotIechey, 
Mess been adopted. However, for shells the Trescasyield 
surface becomes: quite complicated. Although Tresca yield 
surface for axisymmetric shells has been developed by E. T. 
Onat and W. Prager [12], 1t has been used rather as a 
reference surface. To make shell problems tractable anal- 
ytically, some simvlified yield surface must be devised. 
In static plasticity, the general theorems of limit analy- 
Sis allow one to bracket the exact value of the collapse 
pressure and to evaluate the degree of accuracy of an 
approximate yield surface. P. G. Hodge and B. Paul have 
considered the influence of approximate yield conditions in 
the dynamic behavior of a simply supported cylindrical shell 
subjected to uniform pressure pulses having various pro- 
files [13]. Although they made some useful suggestions, no 
general theorems corresponding to the limit theorems have 
been established with respect to the predictions of approx- 


imate yield surfaces. 


The spherical shell which has two equal curvatures is 


fiche soteliegeonetry that was Considered. “R. San- 





oe 


karanarayanan obtained the solution to the problem of 
Spherical cap subjected to uniform dynamic load which decays 
exponentially [14] or has a rectangular profile [15]. MThe 
results obtained are involved and a numerical approach had 


to be used to examine the question of admissibility. 


It has been observed in Eharticity that problems are 
usually simpler when analyzed as shallow shells and that the 
Poles OOeatioad May Give Some andications Of Ehe cCoOrrespon— 
ding behavior of deep shells. This study starts with 
dynamic plasticity of shallow spherical shells. Analytical 
procedures have been used as far as possible and three 
kinds of yield surfaces have been used, of which two only 
give distinct solutions. Then general shells have also 
been examined. The solutions obtained are complicated, for 
in the cases of plates, except for [9], and of cylindrical 
shells, there are two general stresses and one velocity com- 
ponent, for which there is one equilibrium equation; the 
other two equations can be found from the yield surface and 
the flow rule. In the case of spherical shells, there are 
Five general stresses and two velocity components, for which 
there are three equilibrium equations and the other four 
equations can be found from the yield surface and the flow 
rules. It is obvious that difficulties increase with the 


number of unknowns and the order of the equations. 
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PeeeeGenetal Relations 
2.1.1 Strain rate-velocity relations 
For rotationally symmetrical shallow shells [17] 


loaded symmetrically, the relations between the generalized 
| 2 








strain rates and velocities wnen it is assumed that y' <<l, 
; max 
y' being the slope of the meridian of the shell, are: 
a oe 
Cit a eae CZ cll) 
re 
Kg = -h = (2 see, 
gy =y' = y"w (2.1.1.3) 
Ky = -h(w' + y"v)! OD Tcl 7h) 
; we. Cn) 
Wa etl oe a. = oe 
6 
and h = Aes (Zl ee) 
O 
S 
Np Gan Gees ce 
n 
0 
s, 
v 
Fig. 3a me 


Positive Sense of Moment : EoGqee oD 





20% 
In this case, the reference length A is the base radius 
R of the shell: 
A = R 


mivem, = eoOnme2 . ai. 5) page 719: 


15 


#H (2a ee) 
2R 


Zee we Guta baum, Cua Ons 
Using the results from [17] for shallow 
symmetrical shells loaded symmetrically, and with some mod- 


ifications to include the inertia terms, the rates of 


internal energy dissipation and external energy input are: 


l 

dst = | (mee + ae + MK 4 + “PS CL 
1 a ae Le 

eer = | [x (p_ - yw)w + x(D 4 - v)v] dx + 


[x(Sw + Ayv) - hmy(w' + y"v)1 Ly 


p. and Ds are respectively the normal and meridional 


n 


components of the applied load. 


e e 9 


With Cys Cyr Kye Kg tepaO MME ee cali so ld) Celgene so ace. 19 


aia applying Ehe principle of virtual velocity, we obtain 


the following equations of equilibrium: 





xS = PRU ar ss Mo] one + . yn’ yx 


(xn 4)" - Ng - y"(xs) + P 4x XxYV 


y"(xn4) + y'n, + (xs)' + pix = xXYW 

where 2 =w' + y"v is the rate of slope change or angular 
velocity. 2 Raswarysense Of Gotation Contram, Co the sense 
Of positive il eniGmenie | cera ‘ yn?Q, = ~i(-24) where i is the 
non-dimensional moment of inertia of unit area of the shell 
with respect to the reference surface, represents the rotary 


amertia. 


Tf the load is applied only in the normal direction as 


a pressure load, and which is our case, we have: Py = Q, 
P= Dp, and if the rotary inertia We neglected, then: 

xs = eb Ceti - Me] C2 oi eoea) 

(xn 4) | ol = y" (xs) = xyYV (2 lee) 

y" (xn 4) + y'n, + (xs)' + px = xyw (Za) 


For a shell without a hole, simply supported at the 


edge x = l, and made of isotropic material, we have the 
following boundary conditions: 
At x = 0: 


From isotropy consideration: 


s = Q 


Ge 
I 
© 
= 
lI 
© 


or there will be 
a hinge 





Wee 


ier ote ee 


For simply supported edge, not free to move: 


m, = Q 
o) 
w= 0 ’ y=0 or there will be 
a hinge 
The initial conditions are: 
t=O , v=ev=0 ,we=we=0o0 


Poe ViIeversuriaces 


E. T. Onat and W. Prager have developed a yield 
Surtace for a rigid perfectly plastic material that obeys 
Tresca's yield condition [12]. The complexity of the equa- 
tion defining the yield condition is such that considerable 
difficulty may be expected in the solution of all but the 
most trivial of problems. Therefore simpler approximate 
yield surfaces have to be devised. Three approximate yield 
surfaces will be used in this study. These are: the 
two-moment limited interaction proposed by Hodge [17], the 
uncoupled square yield surface which is a simplified version 
of the former and the uncoupled diamond yield surface 


eeoposed by N. Jones. 


= a 
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ghee ete) | Fig. 4b 
Two-Moment Limited-Interaction Yield Surface 
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Fig. 6a Fig. 6b 


Uncoupled Diamond Yield Surface 
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The Hodge and Jones yield surfaces allow interaction 
between force and force or moment and moment but consider 
interaction between force and moment to be of limited 
importance. These yield surfaces and the planes that define 
them and the normality requirements on the associated gen- 


eralized strain rate vectors are presented in the following 


tables. 
T AOC) 2 4 Jlexoek 
Two-Moment Limited-Interaction Surface [17] 
Face Equation Strain Rate Vector 


(a9, Cb, keane Kg) 


li Shi ae Wel Oe, Or Oy) 
2 De = |] eCOy oe 0} 
3 “Nin, eal as ele O10) 
4 “Ny = I u(-1,0,0,0) 
5 “Ay =) u(0,-1,0,0 
6 Meats = 1 u(1,-1,0,0) 
7 aby = ] LCOr Oe so) 
8 my = I WOr Oe 7) 
9 TMG tM, ap u(0,0,-1,1) 
0 —Mg = J Weer =) 
ae ig = Tt COMOROS 


p 
B mM “Mo = | ee, Ole) 





Face 


Face 


Baller al. 3. 2 


DReoupled ogqvare Yield Surface 


Strain Rate Vector 


Equation (0, by, ke, K 4) 
np = 1 u(1,0,0,0) 
al eeO es yOe.0 } 

=<. = all u(-1,0,0,0) 

“ny = u(0,-1,0,0) 
My = 1 e700) 
my = ih Omron a) 

“lM, = i u(0,0,-1,0) 

aS =) 07 O. On 1) 


AWatoyl fey 4 2) S846 


Uncoupled Diamond Yield Surface 


Strain Rate Vector 


Equation (E,, by, i K,) 
any = ] jes Gales SOC) 
Bern aa (a 0.0) 
-n “Ay =] u(-1,-1,0,0) 
Meo = 1 u(1,-1,0,0) 
mg +m, = ik UChr 0 elk) 
ey ugar = l u(0,0,-1,1) 
“M g-M 4 ae u(0,0,-1,-1) 


W070 p= 2) 
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See © 

Peg static collapse pressure of a shallow spherical 
shell structure according to the Tresca yield 
Sie pace. 

Poy = static collapse pressure of a shallow spherical 
shell structure according to the two-moment 
limited-interaction yield surface. 

Pog = static collapse pressure according to the uncoupled 
square yield surface. 

eae = Sigeiee CGOlllevecie cesses sree ibng (Wo “Uli ubelereybloulere! 


diamond yield surface. 


Then according to the theorems of limit analysis 


(e.g., P. G. Hodge [18]), we can construct the following 


bounds: 
ON ss} Poy < Por < Poy (2. 1233) 
0.309 Pos ss Pom = Pos (2.1.3.2) 
O68 Pap < Pop £2Rp (2.1. 325) 


These bounds correspond to complete inscribing or cir- 
cCumscribing of the yield surfaces. The interval between 
Ene lower and upper bounds can be reduced if only part of 
the approximate yield surface has to inscribe or circum- 
scribe the exact Tresca yield surface. It will be seen, for 
instance, that for the uncoupled diamond yield surface, the 
direct stress resultants of a shallow sphere under the 
medtle cOllaose load have a value of -1/2 approximately. 


Then, with n, = -l1/2, n, = -1/2, the uncoupled diamond 


0 p 





PAE fa 


ieolG@eeurtace can belmade to circumscribe the correspond 
mag part of the Tresca yield surface by a multiplying factor 


Sees instead of 2 as stated in (2.1.3.3) page 26. 


2) deep EScontrnulrLes 
From dynamic consideration, it can be shown that 


S, n m, must be continuous whereas n m, May be dis- 


o” 9 | gc 
continuous [17]. 

For velocities, we observe that the basic 
assumption that straight normals to the reference: surface of 
a shell remain straight and eras the reference surface 
during deformation implies that shear strains are negli- 
gible, and therefore w cannot be discontinuous. However 
w' and v may be discontinuous {17]. When such discontin- 
ulties occur, we can see, from the expressions (2.1.1.1)- 
(2.1.1.4) page 19, that the strain rate vector has the fol- 


lowing direction: 


eg 0 
‘ Ne Ax 
S V 
Mie 1D tho il) 
k = 0 
: 
k -~h (Aw '4y" Av) line of discontinuity 


Ax 1s the increase in xX to get from one side of the dis- 
Semcinuity liane to the other and AF is the jump in F cor- 


mesHponding to Ax. 





mace yianto Plasele Response of Simply Supported Shallow 
Spherical Shells made of Rigid Plastic Material that 


Obevoweiewiwe- jomene Limiteda-Interaction Yield Surface 
The applied dynamic pressure is assumed to have the form 
of an exponentially decaying pulse: 


P = P e t/To 
O 


Or, with non-dimensional quantities: 
p=pe. 

We will solve the system of eet 2a alee ae de 
(2.1.2.3) page 21 using successively the three yield sur- 
faces defined in (2.1.3), page 22 In each case, we begin 
by determining the static collapse pressure, both for simply 
Supported and clamped Supports. A solution is then devel- 
oped for the dynamic response of a shallow shell with simple 


Supports when the pressure 1S applied dynamically with a 


peak value larger than the static collapse pressure. 


mee SGhEaAELC Colne seme > Uaes 


When the pressure 1S increased slowly from zero 
so that the inertia forces are negligible, the equations of 


Peeeitorium (2.1.2.2), (2.1.2.3) page 21 become: 


ee h Cm) ' amp utay DR Sileal) 
(xn) Zaliaes y"xs — 0 C7 ae) 
Wg (xn) + y'n, + (xs)' + px = 0 (2°20 le 3) 


The boundary conditions for a simply supported 


edge are: 


r=. tee Si = «0 (2e2.1 74) 





Zon 


Because of the boundary condition at x=0, a Suitable 
yleld regime will be made of one of the ae je ee are oS Sie shgte 
one of the faces 7, 8, 0, a. Because of the boundary con- 
dition at x=l, only faces 7 or 0 can be retained. Moreover, 
Since the pressure iS applied externally, we expect n, to 


v 


be negative and m, to be positive, according to the sign 


® 
convention in Fig. 3b page 19 Ese rie winimiioc Ot “caGices 
reduces to two: regime 4-7 or regime 5-7. It turns out 


that 5-7 1s the correct yield regime. (Fig. 4a,b) 


With regime 5-7, we have: 


ec Se Sa 0, (282-155 Sa < ¢ 
m= 1, Ky = 0, (2.2.1.6) Kk, > 0 
With ny eae (222 lee) reich tt 27 eles) ocicion moe conten 
el Va ne Waids! =O) “(2.2.1 2) 
SI ae OE a (>. <2) he Be a o> oS OE 2) ats) 


0 


Solving this system and using the shallow shell approx- 
mictton: y'* << 1, we get: 


xS = - = px? see ne eG 2 Ihe JUG) 


Ny = - 1 + y" (5 px? - xy") (2c 2a eee) 


By integrating (2.1.2.1) page 2lbetween 0 and x, we 


obtain: 


x 
SEL = ho = px? + | Ey Cicemne 2s eee) 
O 






30. 
Peon ly ouoDOtLeca Snel is 


At x=l1, 7 = 0 for simply supported edge. Thus: 
1 
oa Sy as Gila 6 | SAG SOIR 2627s eres) 
O 
BOura erat molcdee, 7.8 — 07) and wey tind again the col— 


lapse pressure of circular flat plate. 


With this expression of p, we have thus: 





16 
x 
R 
e 
R 
inspy i 
1 
mee= -—- (3H + 3 | Ey 'dé) x? + oven OSG HED) TU a a 128) 
O 
i 
no, = -1 + y"[(3h + 3 | Ey'd&)x* - xy'] OSI A 2 sy) 


e) 


1 x 
Pe ee) ede) x? & =. | cy"ae (O<xel) (2.2.1.16) 
QD h hx ‘ a ee 





Sule 


The flow rules are given by (2.2.1.5) and (2.2.1.6) 
page 29. With Ey and Ky CTO) ih el ea pT) Sys Sool a0 2) a eet!) 


mage 19, we have: 


v= y'w (ool) 


w' + y"'v = A = constant C232 eis) 


Solving, using the shallow shell approximation and with 


the boundary condition at x = 1, we have: 


wos wo a) DD 1 DEO. 


and V 


tl 


y'wo (1 ~- x) (Oso raiee sb) 


From these velocity expressions, we have: 


e = -y'w_< 0 since y' and w_ are positive 
Oo - O 


x [a 


[-w sa y'y"wo (1 See) 


For a shallow shell, y' is small and y" which is approx- 
imately the curvature of the shell is also small, thus 
y'y"(l1 - x) can be considered as small compared to 1 and 


therefore, we have approximately: 


ered 


w > 0 Since w > Q 
o- O 


k = 
The solution is consequently kinematxally admissible. 
We assume the reference surface is smooth so that y' exists 


everywhere, particularly at the apex x = 0. Because of 


Symmetry, we have then: 





324 


For most practical geometries of cap, the equation of 


me middle surface can be written in the form [17]: 





_ on se 72 xo 
Y n 
For y' to exist at x = 0, we assume n > l. 
To determine c, we use the fact that at x = 1, we must 


have y = Z2/R, where Z is the total depth of the shell. 


Then: 





(op ae 


IE 
2 


and the equation of the shallow shell becomes: 





ye = x n>l (ee 22) 
7,2 
with Ga.) = 19,7 Re SSL 

max 

From (2.2.1.22), we have: 
y' = n= nt 
x x 

_ az n a n (ES Neal 

and [ey d& =n R 3 Ge a pee 


Wiel EnNecemexpressions, (2.2.1.14) 7° (2.2.1.15) and 


(2.2.1.16) become: 





Eyae 


_— Borie 2 Ne oe 
ES = -(3h + Sarna =) x +n gx CO” SI ae La a!) 
a. 2 = ey Ge Sein Ee a) aa) 
9 R nti R R oBo de 
(O < x < 1) 
Ree te 2 | x?) (one = < 1) leoEy 
6 Ree ae 3 ame 


From these results, we have: 








lon. fe n-l 
aa se as om an — 
m a= 2x + oo rape (nx 2x) 
me. 0). = Oo and 
o) 
ee ee eee ee ee ee 
ra - h n+ 2 h ntl pts 
m" (0) 2 fe BOwe iia 2 
Thus there iS a maximum at x=0: 
MghY) 5 a oad 
iY = Weitere Se ANG 2foOm: 
aoe x 2 n Z 
m Sahin Oo 
h ntl R 
ieee eee atic ae 
ne i Dw 
h ntl R Se Be h ntl R nae 
1 n(n-2) 2 
as a ie 
(2a eee on) 


‘> 
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then x ae aia we Pomimega= lve for 0 -< x Gail, mM 


decreases monotonically from 1 to 0 as x varies from 0 to l 


and is therefore admissible. 


ook =< nm < 2 £aen mg 60) goes to positive 


infinity, which is not allowed since me must be maximum at 
x = 0.% 


MiSs Mib ie MSGMINthoOmulcuM@t acdmesSsaipile for nm < Z. 


For Ng, “ie Ievekvieus escent (2 ela hye 


' = Z - a pil 3n Z 7 2 Z n-2 
gy = BO Mt [(3h + -27 ®) 2(n-1) 5 x ] 
Homer, We have: 
IL 
3h + 207 2 fae? 
aus © O. were O < see< ee 
a 7 7 D ia aky 2 
R 
*For mg = 1 to be maximum at x = 0, it is necessary thet 


ay be negative in the neighborhood of x = 0. 


ee 
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aie 
Sig | ve 
3h + neo @ 
Z < ol 
2(n-1) R 
there will be a maximum at 
= , Ib 
Sol 
2 3h +r ntl R n-2 
an 7, 
2(n-1) R 


To be admissible, this maximum must be negative or from 


Wee. !.24) page 33: 


Ge aan Shak a Zt ine 
n(n-1) R on [3h SE =f oy n R x, ] <7 lt 
oust. 
(eee Sree n 
2 mp ish + car Re! *, <1 a) 
Since x < l1 and 3h + 3n_ B < 3h + 32 we see that 
n— n+l R R ° 


we see that the inequality above will be satisfied if 


Mee) BE 
aa yaa [3h + =. < 1 


Meecoraging EO Kraus 7125, wage 25),as a rule Of thumb, 
a shell is assumed thin when its thickness is everywhere 
less than one-tenth of the radius of curvature of the refer- 
ence surface. According to an analysis by Reissner quoted 
by Kraus [25, page 229], a shell is considered as shallow 
when its maximum height Z@ is less than one-eighth of its 


base diameter 2R or: 
Z/R < 1/4 





Ore 


From a practical viewpoint, we would expect a shell, 
even shallow, to have a maximum height Z at least equal to 


its thickness 2H. Thus, we may assume that 


and: 
Zi; 2 ee Ze 
O[Z(3h + zz! = 0 i ) 


ama the Condition 


n(n-2) 


Z 32 


on page 35 is equivalent to 


Z 


n(n=2) x 0(=5) <1 


2 


For practical value of n, we can expect that this 


inequality is satisfied. 


meex = 1, we must Rave also: 


a Nas Nop Sa0l 4 Cle 
0 < n(n-1) = [3h - n tne 2) =] coe 
The condition: 
Co Mina, = [3h - nin 2) =] 
1s equivalent to: 
1 iil ols) 5 3 





Sie 


ana it 1S satisfied 1£ (2.2.1.26) page 33 is. Theremaining 


Send Gaon, Which is 


ek n(n-2) Z 
n(n-1) R peia! ares R! alk 
ramieleije SatiStieca if we have: 
Sadao) Sie 
R aeees 


2: 


Seee = h is of the order of ar which 1s considered 


negligible compared to 1, the condition above is expected 


to be satisfied for practical valuesof n. 


Thus the solution is kinematically and statically 


mammssibole for shallow Shells of the form: 


y=ox r 0 Medea” 
iemecenartion (2.2.2:1.26) page 33 which is: 


ig (eZ) 
ata R eck (232.2126) 


few satistied, 


The static collapse pressure is given by (2.2.1.13) 


page 30 , which can be written as: 


_ _ eat A, 
Dit eee 6A + aan CSE NAR 5 LAT) 
For a shallow surface of second degree, n = 2, which 


may be a shallow spherical, parabolical, ellipsoidal or 
myoerbolormdal Cap, condition (2.2.1.26) is always satisfied 


and the static collapse pressure in this case is: 

















38. 


p = 6h + 4 = (252 wiles 
Be Clamped Shells R 


ieee Ssiel] 2s clamped, the boundary condition at 


x = 1 becomes: 
w = 0 
v= 0 , wio= 0 or there will be a 
hinge 
m 
w) 
8 
9 fj 
mc —_ ne 
3 8 





We Would expece w! +0 at x = 1, and there would be a 
hinge. Fromthe considerations on discontinuities in 2.1.4, 


page 27 , we see that this is possible if [| =f 
We therefore assume: 
On <5 Xe eee Commer S17 


Ue Xo ig eco eso, 











oo elo ors0 = Xs us regime S27 


Soe. 


For a spherical shallow shell, the value of n in 


W272.1.22) page 32 is 2, and we have: 


yae x? 
fei tlle y' = 2% X 

n 28 
and | y R 


(25 Zale» 


(ARPES! C,) 


(22 oe) 


With this, we have from (2.2.1.10) page 29 


emieme2.2.1.12) page 29: 





eis Blo 
xS = 5 (p R) x 
| — a 6 ON 2 
“ihe = 1 + R (p R)x 
_ edge _ 42, - 2 
the = ] ran (p R) x 


oe 


2s 


eZee 


For the velocities, the solution of 


(2.2.1.18) page 29 , without using the boundary 


at XxX = l, yields: 


ana: 


He Zales) 1) 


2 sles) 
Papelera) 


2234) 


(2.2.0 Deel ee 


condition 


2 a S)) 


oe) 





40. 
Peco licen tom u < x < 1 : regime 5-5 


With ng = -l, the equations of equilibrium (2.2.1.2) 


Meelis) page 289 still yield: 


eel _ 42, 2 
xs. = 5 (p mn) x (2.2) 
(Ome = 1) 
ngs Z _ 42, 2 
Way = 1 + R (P Rx (2277 he 33) 
(O;-seexn sas) 


From the relation corresponding to face 8, which is 


mm, - m4 27 ‘NMfablef2= is.) 


and from the equation of equilibrium (2.2.1.1) page 28 


we have: 
m' ~~. a ( BENE) 2c 
o~ x 2h ‘P R 
Integrating between u and x, with m (a) = 0, we have: 
= x _ 1 ie ee lag tee a 
m, = log 2 - a (P aa) Ox u2) =(2.2.1.35a) 


= x _ 1 ae 2 ee ee 
m, = 1 + log aa (p a) oo) CL 5 SS EY) 


For the velocities, we have the flow rules: 
ey = 07 Om ew thnmme2 lel) “Dace 26 


= y'w (egies? 


<e 





4l. 


and K + Kg =—N—wote rw oecomeSs, With (2.1.1.2) and 


@ol.1.4) page 19 


(w at ve aoe 1 
w' + eG XxX 
and which gives: 
up oe yy ON (252) ieels) 


solving Cras 7) caiG) waza eco) bORaX. w and v and 


using shallow shell approximation, we have: 


W 


B log x 5261, 2) 


v= 2x B log x . BOAO) 


Bemeoundary Matching 


We have four unknowns: A, B, p, u. To determine them, 


we have four conditions: 


(u) = 0 


Cena 


ie) en 


p 
? 


(Gai) eon elmudacy, Oo wat x= u 


Cra) 


(iiii) Continuity of w' at x = u 
Y 


From my(a ) = 0, and (2.2.1.34) page 39 , we have: 


Oo’ 


35 2s ne (a eal) 


G 











42. 


Peon nei) ———teecaiica (2.2.1,.354) page 40, we have: 
4h 4Z 
p = oe (1 - log u) Ey (22.20 2) 


From AW Oma ncemticnmncm.2 021.35) page of and (292 2en oO) 


page 41, we have: 


e- 


je) ere) tea Wo ( 22 ees) 


From Aw' = 0, and Glee cS heise Viocpe 2, 2. ills aie) Seis ste 


pic 2.2.2.1.39) page 41, we have: 
B/u =A OD ere 


Solving for A and B from (2.2.1.43) and (2.2.1.44), 


we have: 


O 
at Sedoa ay (22a dies 
we 


and B= - T - logu (22 a Lee) 


With these expressions of A and B, the velocities 


become: 


Paemet 22.135) and (25224 .36) page 39 


. x 


i wo tl = Hl Sloe wre (0 —— Sah, (2 ee 
we. 2a [pe os nk 
eo ze Wx Glee oa aE CCS 2 oy) (25254 8) 





page 


have: 


2) 2 . 


ae 


He@iignc. 2.59) and (2-2.1.40) page 41 


ae log x 
ao Wh i lee Ue Se 2 ee 
Serre Cae ge SLOG) ox 
V= r Te ieee u Ses 12 eee 
Eom 2.2.1.49), we have: 
W 
Ree) oe 
lth = iexe; 1b) ‘ 
W W 
and at x =1,w'=- ees eee 0 Gs = 
a eo cat : jhe) 
iiWietews 2 pange Cine lerae sy = Pe wa en kK <7 0 vais 
required 


The value of u can be determined from (2.2.1.41) 


41 and (2.2.1.42) page 42 by eliminating p. Thus we 


25 = 5 =- log bie (e2en 7 selena 


From this equation, we find: 


u? 


R 


eS 31 and 


R 


u Oe oa C2 ee oe) 


Then P. 1S given by either (2.2.1.41) page 41 or 


1.42) page 42 


_. _ 4%, 6h 
pe Ps R u2 





44. 


4, Admissibility 
1. Kinematic Admissibility 


a. 0 < x < u, we must have e, < 0 and k, > 0. 


ob 0 


Wen Mawr rom. Zee) ane (207 kaa ee 


ey aie 22 /RT eo 0, Sees Wo > 0 
and: 
kg * : rf @ NESE Wes ai 
yo WS ae SEI KS ie Isic: @ pon r= 0 0X! k > 0. 


p 0 


Hissiienew | 22k oo) eainGen(2 7g oo ee pagemsa 6 


we have: 
7 ah - on <u 
and 
kg gO tym) = tee > 0 


When k. is admissible, then k, is also admissible. 


6 p 


Thus the solution is kinematically admissible. 
ii. Static Admissibility 


a. Ny, we must have -l1 <n, < 0. 


8 
BPMOn (202.1. 3 5) pegeuwo7, and Using > from 


(2.2.1.41) page 41, we have: 


2 
ng = -1 + 6h = x (DEA Sey 


on 


a 


_ 
7 
> 
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From this, we see that Ng increases monotonically 


and at x = l, we nave: 


pomene tea? aR 


which is still negative since h= 1S very small compared to l. 


eye OU 
p 

POOR se <oeue S iene race Olas ov Sale 
going( 22 ees 4 pace oo lence lcasnca omen com 


(2.2.1.41) page 41, we have: 


(Olea) ) (2.2) 


From this, we can see that m, decreases monotonically 


® 


mae to 0 as xX increases from 0 to wu. Thus m, 1S admis-= 


Y 


mele in this interval. 
Hoists < x)= 1, we muse tavye =i) < ail ee 


Huon (2 .2elnooa) page 40 andiucingp frome les) 


page 41, we have: 


x) 3° s=u- os 
MH = log ei a annua Gs Oe = 0 
2 24,2 
mt=ss- chee as <0 See ee 
d Xx u u-x 
Thus a decreases monotonically from 0 to -l as x 
increases from u to 1; m, is admissible in this interval. 


® 


The admissibility of Me follows, as a consequence cf: 
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fm Conclusion 
The solution is an exact one for the yield surface 
eweeted, and the static collapse pressure of a shallow 


Pemerical 1s given by (2.2.1.41), page 41 


_. _ 421, 6h 
eo Gee 


menu: Getermined from (2.2.1.51) page 43, which is: 


ae = LOG u? 


G 


and which gives 


2? = 0 S88 approximately. 





47. 


22 ee amt ee Kesponse Of Shallow Spherical Shells 


eee ew merece aC Raniges 0s DOs Pels= sine 


If the peak value of the dynamic pressure is 
sufficiently low, the whole shell will collapse in one 


regime: regime 5-7. (Fig. 4a,b) 


lee sxe bisakena: 
VELOCITIES - The flow rules and the boundary 
conditions are the same as those in the corresponding static 


peoplem. From (2.2.1.20) and (2.2.1.21) page 31, we have: 


w= wo(l - x) (Do Pil, 20) 
ae oe 7 
VES ie wox(l x) Ole 
Wo is the velocity at the apex x = 0 and is now a function 


Of tT, and y' has the expression (2.2.1.30). 


SDRESSES = FrOMsene “equal ioram equations 


freee 2.2), (2.1.2.3) page 21 and with n, = -l, we have: 


p 
-l - y"xs - Ng = XYV (7nd ela) 
-~xy" + (xs)' + y'n, Ta OPS Any C2 25 7e be?) 
Moreover, from the flow rule e, = 0, we have: 
v= yw or 
ee A Sy (2272 ls) 


Sonnac eee lal} 2 W252 Ae), Wee ee oe Soe eclinG, 


Peeehe (2.221220) page 31, we obtain: 


48. 


ese — i — 4Z ie “5 xe — a 
xs = 5 (p x + yw =r (2 t2i eee) 
and 
a PIA Si: a ee ou Sn eco 
Diy 2 al ay eee eae) cP Cae, aaa mae 
(2,.2% 206 ceo 


PeOniecouatlonror Equallibimium (2, a2) page Zi. wate 


My = Watcitteahating Ereom O71 FO x, we heave: 
= = I _ 424 2 ws oe 
Ms = : En \P R)* + eh (Se 5 ) (22 2. oleae 
From the boundary condition: x = 1, cn = 0, we can 


determine yw: 
en 8) Ded Oe ae ee mee 


mreeqrating (2.2.2<1./) with respect €o6o tT, where 
p= Poe and using the ‘pede COndlELon: <1 = (05 Sw 2 = 70; 
we have: 
") 


yw = 2{p, (1 -e 


: = (Qe25 lee) 


FINAL RESULTS 


Wiieh 2. 2.251 7) and Uzeacdscel, 6) vabove..7e Mave 


meen (lsc. co eanad. (252.052 apace. - 


yw = 2[p (1 - eT} = Sgree i) (Qe aon 
2 2205 (ree) = oe) eles). (@.2.251,10) 
oes al Uo P. poe tele 

a = ~l, Me = 1] 
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momo e2.2.1.4)page 22, with- 


we have: 


ON 


ej = (p 


2 


MOM 2 eves se 5} 


Gilee 


=, eee 
n, = -l + 2 | 


6 


MEO. 2 oc eo) 
6h - 





We 
TH Se one 
ace es 210 =epe) 
ee a (252°. 2 Sica) 
page 48 
5 je = Ser), 
> 2 ) = 3 
: x* + 3 (P Pp.) x ] 
b2. 2 eee) 
page 48, 
(omnes) oars 
== x? 7 —— roar 22) 


2. Admissibility 


1. Kinematic Lomas Solely, 


We must nave e 


i S710) seta! Kp > 0. It has 


been seen previously (page 31) that these conditions are 


Seeastied for Wo 


p.(1 - e 'f) - pit, = 0 


eS = 
At T = Ter where T 
W = 0, hence Ww = Vv = 


fs 


PamOh, Cig econ (45 2-2u le 3) y oeciees & 


2[p (1 - a oe p.t] > 0 


is determined by: 


(26, 2592) 
S 


GO and the Mervlon ceocoes. 





a0: 





Fig. 8 


Graphic Determination of T 


I Ora eVe Ramis sina ba ty 


a) m,,condition for admissibility: 0O<m 


) 


Bion (2.2 .2e).13) page 4, werderiver 


= 


m', = - gf 26h - (p - p.)] + 3(p - pg)x } x 
(22 eo leo) 
and 
m", = - g { 2[6h - (p - po) ] TOD ee eco) 
(2 62a LS) 
At x = Q, 





Dele 


To be admissible, it is necessary that there be a 


Maximum at x = 0 or 
(po; Pp.) - 6h 
Brae = << 0 Gus 
p x=0 3h 
Lea De < 6h 
Since p = Poe decreases monotonically, the above 


Mmeejgality will be satisfied for any 1 > 0 af it is at 


fe—- 0 OL: 
ae 6n (2p 2 = la) 
If Po 7 Ps . Chenmewe.) 1 = 20>) Km ) x=07 Oo; Siow 
oo) i, = = = < Q and there is still a maximum at 
od x=0 h 
x = 0 
iaeKi(2.2.2.1.15) page of at can be seen. tmat mg = 0 
for another value x Of: ae: 
2 6h - (p - p.) 2 6h + (p. - p) 
xXx = - 5COOeOoooOOOOOO CeO OOS cS EO 
m 3 p- Pp 3 Pas ae 


(Zee2 «2 yleaieey) 


which is the position of the minimum of m 


o: 
Assuming (2.2.2.1.17) above satisfied, then: 


IE) Ot pa? 0, X. 1S negative and M ¢ decreases mono- 


mM 


fomacally in the interval (0,1) from 1 to 0 and is therefore 


admissible. 


Pee eo eye 1S DOSLEIVe. FOr Museo be. admissible, 
S 


Q 


Xm must be outside the interval (0,1), otherwise mg will 


mi 





DIpae 


decrease in the interval (0,x_) from 1 to a minimum (m yee, 


mie ancrease from this minimum (Me) in C6. Zero etic 


mneervyal (x ,1). Thus this minimum (m,).. , if it exists, 
m | d° min 

will be negative, and make Mt not admissible. Thus, we 
must have: 

Sac 

. mM _ 
eee trom (2.2.2.1.18) page 
12 hae en a= S30 
a ee 


qa 


Since Poe decreases monotonically with T, this 


inequality will be satisfied for any T in the interval 


[O,t-], ede Se al Tp. Moh ee} © 


12h +p.e f= p. > 0 (2.2.2.1.19) 


S 


where t. has been determined in (2.2.2.1.14) page 49 


b) n,, condition for admissibility: —l<n,<0 


Q/ 


fom (2.2.2.,1..12) page 42, we have: 


ee 4 _ a = 

Ty Se ( [Siw Bo = je) & BUS <1.) } 

(2a 2.2 ey 
and (n' 4) = 0 
x=0 
Since (Ny) 9 = -l1, there must be a minimum at x = 0 or 
(n 9) x=0 2 10) Wola hoe se Sfeqbel sae 
=i 
Poe Pe. Ss Lee 





5) 
Since Poe decreases monotonically with tT, this 
inequality will be satisfied for any T > 0 if it is at 


fe— O, Or: 


Be Se, S teas @.2oPoN 20) 


From (2.2.2.1.20) page 52, we can see that n', = <0) 


for another value ae OF 32: 


3) = 3G) Sak) 6h + S(p. - p) 
“n 8 (p = Pp.) — 8(p. = p) 


(Cae Pras ray Pe) 


Mach 1S the position of the maximum of ng. Assuming 


We2.2.1.21) page 53 satisfied, then: 


heap — Pe De AE) 


‘ is negative and n, increases mono- 


6 
tonically in the interval (0,1). 


ie) — P. < oF x SRDOS tive. se x < Or : 
6h - 3(p. - p) <0 OA ENS BE!) 


ren Ng Maso eles mun ects x HONORS ede acon (22 2 ale iL Z)) 


page 49, is: 


2.1244 2 
(n,) = -1 + 5 pléh + 5(p, - p)) x2 


n 
0° max 


To be admissible, this maximum must be negative. 


Sieiee ee it I oe will be negative if 
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roy) 
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te 
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—T 
Since p = Poe Ceehedses MONOuem cclely an it liaam: 
Pers inequality will be satisfied for 0 < Tt < Te, iar 
at T = Tee 
z eg R 
6h + 3 (p. Poe ae 7 (25, 2 eee) 


maemi2.2.2.4.19) page 52 we have: 


p.- eS leer 


S 
Therefore: 


6h + 5(p. - p e ‘f) < 6h + 5 x 12h = 66h 


O 
prac Onacition (2.2.2.1.24) above will be satisfied if we have: 


66h < 35 (2.2.2.1.25) 


Since h is of the order of a few per cent, let us say 


10%, and - is at most, about 1/3 for shallow shell, we find 


that, for these extreme values, we have: 


As 2 Bue etl 
ae 
35 = 9 
and 
Bain < 35 


Thus, we can safely assume that condition (2.2.2.1.24) 
page 54 is always satisfied for thin, shallow shell. At 


x = l, we must have 


-l < GaN ee) 0 (oka 





Doe 


Peemmc2.2.2.1.12) page 49, this condition is 


18h + p -— Po 


For the first condition, we must have: 


18h + p - DE > 0 


is eeceonGLtlen willl pe Sacieried for 0 a ee tT LE Di oe fas 


at T = Te Or: 


18h + pe 'f - p. > 0 (2 ney 


S 


Me emCOMpaAre this Condition With Condltron mt2.2q2. lao 
page 52 , we see that condition (2.2.2.1.26) will be satisfied 


iommes. 2.2.1.19) is. 
For the second condition, we must have: 


18h +p -p. < 35 


Mars COnNGItion will be ‘saktistied for 6 Sol as Ts li eS 


ieee 0 Or; 


R 
ih +p - p. < 35 Cae ieee yy 


From (2.2.2.1.21) page 53, which is: Darra Ae TeaZins 


we have: 
Seer Pe ares aS) be hy ais 
and (2.2.2.1.27) will be satisfied if we have: 


19.2h < 35 GAD 





oO: 
Comparing this condition with (2.2.2.1.25) page 54 
Memeo esee that condition (2.2.2.1.27) is satisfied since 
Mee. 2.1.25) is. 
Bee eee On Cl Wis anenie 
From the previous study, we have seen that: 


wes acmLlSSibDlLe Lor: 


> 
De Be 6h (222. 2 ee ee acien om 
12h + pe f - p. > 0 (2.2.2.1.19) page 52 
“Np 1s admissible for: 
ewe Oe 16 222 2 Dace 
18h + pe * =p. 2 0 (2.2.2.1.26) page 55 


We conclude therefore that the solution is staticaliv 


admissible for: 
Po ~ Po < 1-2h (2.2.2.1.21) 


and 12h + pe “f - p, > 0 QePa A, kei) 


3. Maximum Central Displacement 


The displacement iS maximum at x = 0. fIts 
maximum is reached when { = oe Pncegratinouls . cee eee 
page 48, from 0 to.t, with yw (t=0) = 0 and using relation 


We2.2.1,14) page 42, we obtain: 


— a = 2 
ene = 21 (p, P./T, jOney i ] (2 eles Oe) 


Sf 
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Generalized Stress and Velocity Distributions 


4. Energy Absorbed 


The energy absorbed Babs by the deformation 


Sumene shell is: 





t Jee 
E — | | PV GGG at 
abs n 
Oo 0 
. oo Ww 
vs fomeniowiermale Veo lLocitey: on zee ss Te 


ado is the surface element: dg 2Qnrar = 27R*xax 


2 


N 
P is the applied pressure: p = = Pp 


t is the dimensional time: ft = a and dt = T.dt 


Die 





Overy 
Then: 


; Oe ly 
ae = 277R Ne [ pat | sexo 





with w from (AEG EVAN ERE) | Sjexelonay Hsrem cshe\e! yore One we 
have: 
4uR?N T¢ IL 
E = p Gomme Glee ail) e tle ‘dt (1 - x)xdx 
abs Y e) - O Pe z 
2m7R*N, Po” T T 
= ea ou, ad Gai ee ne = ee 
Habs 3Y 2 o = PoPs i Ces sea (E 
(2 2 ee Ne? BD 
4. Final Displacement Components 
Witin( 202220 la? dapage So mend ome eierOn 
Pyome2.2.1.21) page 31, we have the following dis~ 


meacement distributions: 
a) Normal Component: 


Le ZNWS) Da) Wg D tae one il =] 
OA 


b) Tangential Component: 


_ 4Z ae 2 7 
Ye = 4—[(p, - ea eee @. Saks pee (ll x) 
(22 eel 0) 
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Final Displacement Distributions 


5. The Case of Rectangular Pressure Pulse 


PES 8 ics evovo ise lie sh jorechers tis OU Meier Wena as) 16 = 


Gime Variation: 


Oe ti < Le p= Po 


From (2.2.2.1.7) page 48 we have: 


ee 2 (Oe) (0 is) 


[A 
=a 
[A 


A 
loa 
eas 


ee 2(po - p)t (Os 1) 





S10 
and: 
Wo eRe Cry 


Ves OR aor ty (ea) 


Then we have: 


yw = 2(p. - Py)t (1 - x) CORES Se oes) 
YV a ie = Dn) Xl =o) ‘CO aumeree <2 3) 
R O S ote pe 
and: 
yw = 2p. = PT?) (1 - x) (Cr emir) 


_ 42 _ 
YV = “p(P, Pot) x (lee) (i ee) 


The solution is kinematically admissible for 


Po 
Ge < too where T.=-— 
— — f£ f£ p 
S 
aa Tes the motion ceases. 
The stresses are: 
pimeeerrom (2.,2.2.1.12) and (25252... 13} 3 pageud9: 


i 


o7 < = 2 (Pp, = P.) 
-L+ = | ——_,-—_— 


y) oo = 3 ‘ 


a ae (Pg - Pg) x2 = oes x 3 (osc 
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ies Static admissibility conditions are: 
OS Sy ant oes WA Aey 00 se 1) 


pi@en corresponds to (2.2.2.1.21) page 53 and 


Po Ss de2 a 


meeereecorresponds to (2.2.2.1:19) page 52; with (2.2.1 .722)— 


z a7, 
ee 


mers tatter condition is equivalent to: 


HS 6h 
Spewitn h = H/2R = 7 a escemml 231.0 1.516) 3 
3 


Thus, for the solution to be admissible in the case of a 
rectangular pressure pulse, the height of the shell must be 
smaller than three-eighths of its thickness. This shell 


could therefore be considered as an imperfect plate. 
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Becac.. Otedium P fe IRE : as 
0 Ger eligem Iaice Jka Zin cs Po ees 6h 
A. First Stage of Motion: 0 ig a | 
Ge ee ee ee 


When De a Daa ik Boy, Sele eseumciareyce 


fWee2.2.1.21) page 53, we have Ing) 9 < 0. Therefore n, 


mould decrease from -l, which is not possible. 





We first assume the following yield regimes: 


Oe regime E - 7 
i le ese ae ae i 
The solution was found to be statically admissible, u 1s 


determined by: 








ara 


Meee decreases monotonically from its initial value Uo 


eae 1 = O: 








xe) Omen) —  t 


SG 
i 


1 ® 
However, this solution is not kinematically admissible: 

a = a . bd 
while e, remains negative for 0 < x < u, e 


0 > 


Meee ~ X < 0.8u and becomes positive for 0.8u < x < u, 


is negative only 


which is not admissible for regime E-7. 


Tt turns out that the correct yield regimes for this 


stage of motion are: 


Oye 2 Uy: regime E-7 
< < 

Uy - x _ Us regime 4-7 

un < xX < al - regimes sa, 


Mv SOlLUtIOnNS are then: 


eos XS uy: regime E-7 
ny = Ng = 1, €y < 0, &y < 0 
m, = 1, k = -h(w'ty"v)' = 0 (2.2.2.2.1) 
With ay = Ng = -l, the equilibrium equations 


me. 2.2), (2.1.2.3) pagel? become: 


-y" xS = xXyYV 


(cs) = KV Sy! soe 





64, 


For a spherical shallow shell, we have from 


fee. 29)43(2.2.1.31) page 39 


Z Zo 2a 
Vis ome va _— Ros ae = 
Then 
7, es 
~ z eS) = XYV 2 1c ee ee 
2 a7 ee 
goo aie (= = 1) ) ae ay ree eee 
Let xs = Y, and solving the above equations for V and 
w, we have: 
a PMY Y 
UE Te So (Zn ee 
ae o. 42 5 
ON Stee ae sa (24.222 2) 
From the flow rule kK =O dm (25202, 2.40) eacgemas 
we have: 
: : : DeLee ae ae 
wht y"y = So : | 
x 


where C. 1s independent of x, but may be a function of T. 
Dirirerentiating (2.2.2.2.6) page 64 with respect to T, we 


obtain: 


es 27, nee 2° 
! Sauer ao 
ywl + ae S. C22 SP aay 


WineeErome2.2.2.2.4) pagenc= anew yw) seoyedtn tenons 
mereima (2.2,.2.2.5) page 64 with respect to x, equation 


(e2ac.2./) becomes: 





2 4 
27 Pinas, (cas ah ees ae 2 =  (C 3 

> aa es a 42 ee C22 ee) 

22 

WEth € = R*, we obtain the following equation: 
Z @ 

eo me Stee ya oe (22252 =2me) 

AE dé pone 

lee 


With the following conditions at x = 0: 
aie SS Se) corel elec ctiecmt (eee 28) 


page 64 and since p and yw Slig(S wee Inolabieres 
at x = € = 0 yy 2350 


PHemMOMocenNeOUS ecuartonwon (2a202. 200). 





has a general solution Y,, expressible in terms of the modi- 


H 


immed Bessel functions of the first and second kind of 


order 1: Ih, and Ky [19] 
Yo = E[c,1, (&) + C4K, (é) ] (2aee 2. ee) 


ime Weonskian of Ely. and EKy ee Cetanead oct 


aa 5K) 


_— 


W(EI,,&K,) = 


qd 
(gI,) Be eo 


00. 





a8 
Using the Abel's formula and using the approximate 


expressions of To? Tye Kor Ky given in [19], where Iie Ky 
Mere w modifi Gd Bessel £Eunctions of the first and second 


kind respectively, of order n, we have: 


E(1,K, + IR) = Gil DO Wh 
ana W(éT,,6K)) Ss 3 Cea 2 2. 2 eee) 


helene Enis value of W(S1,,5K,), and using the method of 


Wemiacton Of parameters to solve equation (2.2.2.2.9) page 


65, we obtain: 





é S 
—_ O =o 
= SZ27R? a, | nk, (n)dn 6K, [ ni, (n) on | 


cf c, (éI,) x a C5 (EK) 


iMitzegiating by pares: 


E 2 | ane 

| mean = - SK. + [x (man, (cimee at = -K,) 
é é | oem 

| ni, dn = a _ | Toman, (since ae 7 I,) 


Preenioting that E (I,K, = TK) Se ee eae. 


C 


é é 7 2 
O 
Y= a799K -£ + EI, | Kodn + &Ky | ran] PU otek 





Ce 


To determine Gy and Coy we have the initial conditions; 


since lim [EK, (8) ] = anal) 
E>0 


E g 
and em | K dn seme lestann| | —1Oo sien 
GOO Sa). 0) 


= lim [-n log yn + ihe = 0 
E+0 O 


we have: 


C E 
ee ee eat 
—— | 1 + E (I,K, + K,1) + eI. | K dn 


OQ, 
He 


d— ~ 823/7R3 


E s! 
= BK. | Iodn | = fait 


O —_ 


Oey with (2.2.2.2,11) page 66 


Qu 
ae 


6 E E =| aes 
O 

ah Sa | K_dn - &K | I an | ee 

dé 8t3/R3 L Oo} 0 Oo} 2 eC: 

( Dea G7 .. 2. 2) 


iene ei) 0 for any eli nnte svaluemer ort 


therefore: 


Sse 
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E E S = 
Q . 2 « 
| ne jaan | SC Shi Meee ek | is en tae ee ui 
(E = ag (Qa.2 50 eo a) 
meeme(2,2.2.2.4), page 64 we have: 
C g g ; 
ee O _ _ a7, ee 
a ram | “d+ ty [Xan + [roan | az C1, 


(22a aa ee) 


Mreom02.2.2.2.5), page 64 we have: 
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ia >  “R * 27/7R | x, | Koen FG | to2n | sae 
Me R en SL: Ss) 


Maeegrating (2.2.2.2.14) and (2.2. 29255) eage oo =e ren 








ieeepect to t with, at tT = 0, v=ew = 0, we have: 
: Cy Ay aes 
yv = - D77R =a + I, [ x, (ndar + Ky [toiman | - B Cit, 
Mame =), (0 < x < u.), § = 2 x DAP De) 
ee. me 1 tf on, ps It I 4 R e 6 e 


+ or CI ODD Aa) 





Ooe 
iWiewmmicial conditions: 
VG SO) = wilt =O) = 6 
will be satisfied if we take; 


a) = 0 Oy Poe) 


J 
co) 


eT) (2 aa 


From the equilibrium equation (2.1.2.1) page 21, with 


My = 1, we have 


eV :¢ 100 ge ame 1G ae 


d 


Paerineecma elie trom  OPtOuc ,) Wi cin Wer ON 27a eles) 
page 49, we have: 
: 


Mee ne + oo, 1 SD 8 + 
ee «6827 Re 2 


E n 
ni, in) | | xotmeny | dn 


E n : a 
+ | nk, (n) | | To(ny)any | dn + Cy [ary (man, 
O 


0 


Or, after some integrations by parts: 


os g 2. é 4 
=>) oe | 2, Mie | Ke Ue [ Tan 
: (n,)d d (Di D7 a a 0) 
6) 5° Lge F 
° - _ 22 
CE oo x) 
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- | Toi) | | Kin any | dn | +- Ci nijdn 
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Siallow Shell Approximation 


Z 
“With y' = (<2 Se) ee wes have strc wiselo teaa, 


limited series expansions: 


2 


e- 
T(§) ile ae 


S 
2 


Ie 


. Z Fe 
K.(E) (1 +45 +S 


a (Y 5 + log 4 


Where Ye 1s the Euler's constant 


Es : 
Va= dim | ) 7 =2 Og n | = 0S 2 


n> p=l 
ee & 
mee) - 35 * i¢ 

_p ee 7 ee Sn ds _ g 
- (n)dn = E(1 - od ne gi (4 - = jee 2 
5° 2) Sli eee YR Bees Wes otis SING, 
- Gallia = i ig: 
‘ge n)dn = & T < 


trem these results, we obtain: 


1 oe 
Y = xs [ ne ©, + 3 COX } (22 ee sc) 
: 22 LL 1: 
NV Se ae MON Se 5 Cae (252i), 





Jue 


* — 2 e 
Yw = pia Ag ac = ede a7 Cy + Cox 
CL AS) 
C ev has 3 il \ 
SER SES Zz pe © % ey ss (22232 ncr 2 4) 
Zz) US XS UD? LEEKS VININGS = 7 
For tnis regime, we have Np = = lag, My = and. thie, st lew 
rules are: 
ec, = 0 +v'i- y"w = 0 (22 28 ey 
K = 0 >w' + y"v = A/y (22 2 
where A is Piae Pendent Of K, and May be artunctloneo ee 
Eliminating W, We Obtain, With. 24,7 2 cel 
wee. i. 3l): 
(yv)" + (A8)2% (yo) = A (2.2.2.2.27) 


mnemoeneral solution Of )(222.2-2.27) os. 


- oA DY . 
Sere B cos = x + C sin = x (2.2.2.2.28) 


Mie, £roMm (2.2.2.2;2o) page.) ©) 
Mie — ee OS eat x re B sin <2 x (2a Bee 


B, C are constants of integration which may be 


iairctilOons of T. 


Hacwmequattons of Motion (2 allng.2) 4) aoe 3) 


become, with a -l and y" = 22/R: 





oe 


(xn 4) - £2 (xs) oy (Zee Oe 


22 : t = 24 — yyw 
Re txN 4) aS) eee =p a= XYW (22 a aD ae) 


Pen ieaclOn Of XS.from <{2.2.2.2.30) and (2.2.,24223) 


yaeids: 


DSS 
R? (xn 


27, 27, ‘ ai * 
) = - —(p - Rx ie ee a2 i] 


nie | - 


(xn 


Y Y 


Mewang yvand yw by differentiating (2.2-2.2.28) and 


Wee 2.2.29) with respect to t, we have: 


A 





r ZL ZZ 27, ss 22 
(xng)" + (Qe) " Gm) = - lp - Rx + g77R + B cos Bx 
ee Za le ee: 22, a 20 
+ C sin 2X + Rx lc cos xX ley Sale) Rx) 


OO Pe SI.) 


mine, General solution of (2.2.2.2.28) 1s: 


255 22,4 BR gin 22 A _ p- 2YR 
a =. COs R* + FE sin R* + 82 37R3 227R ~ 
mee Viel kien e 2 22 


Deanc & are constants of integration which may be a 


MNCL OneOL TT. 


Paomaccuatton (2.2°2.2,30) we Nave nen: 


> Ventre ae ee, ee MBS AS tele ys 
= qZ2/R2 Ag? /Re* hel: 3Ces R* Dys2n R* 
+ xx COS xX sx" Sin =x C2 Zed ee 
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x 4Z 
where Q= p - pr OF with P. Er OMe 2 ee el ee See 


Ee 2 Secs 1a ah OIE IT OX, SS) ) 


PEO ene wequidibrium equation (22152 5a) 





with Mp = 1, we have: 
A (XM gy) = ho > xs and 
a8 = Q - A 2 EB ah 
an) = 1 s> ox 422 /RE* 822 /Re* a DZ7R 82D RX 
D 27 C 2x 27 
* omy RO tees = oy sna Se Sa 
Y4z2/R%2- 2 . 22 
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_B ___2x on 2hy — 42?/R'x?~ 2 O45 22, 
Dee 7 Re 87 3/7R: =e 


(292 52) 


epollow shell approximation 





Using the shallow shell approximation y'* = 


we have the following series expansions: 
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Dim ie ME bE 
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Sin —x = is 
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foascipetmese results, we obtain: 
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yw=C- af Bx (DeDLDED LED 
ei ee AZ 27 2 eee 
ee yy RO RP oe Lk ee 
AD 4 Se =x dO c= 
SZ Re DAR 2 Oar 
a7 2 oe Ve 3 lL 7, 7, 3 Ve os 
Oi: A ENS 
2 ae 7 3 
eb Ate 
Bammer era a (Cae RZ jh ae oe 
Me Re 47,7 /R? 
Bee 2h Gey ee (2.226252 2400 
6 Rr? R i Linear ae 
aS core ee eee me oe rr 
64,8 + 82eb - 228 =e-E - @ A+ S4-D 
hxm HE ot — ] x - x’ 
167 7R’ Lay Re Bo Fe 
sia - 42° 8 x84 Dg BG - 6 2B) x! 
6 R2 24 R R 
C2522 G2 
3) ae CS I icterep inva. ai /. 
For this regime, we have: 
Op Ss t Mp = 1, 
exc tne flow rules are: 
Cn = 0 and = = Q 


With the solutions from the low pressure range, and 


Sbserving that in this case, we cannot use Ehe boundary 





eae. 


Pemart1ons at x = 0, we obtain the following results: 
e 7, e 
Y= — «CKX (x - 1) (2350 eee ee 
yw = K (x - 1) (220.0 oy 
eee RO. Se 
Np 1 R [G h 5 + We | 
(2a ae oe) 
xs = G@-h- S52 x? 4 3 x3 (22m 45) 
_ de iL. ae _ oe we K + O 3 ws aye 
Rea = 75 K + 7 QO G + Gx iG ae + 1D Kx 


C2 p22 Op 


G, K and its derivative K are constants of integration 
Wiel may be functions of T, O is known and has been defined 


mmee sf. 2.2.35). 


4) Summary 


Met us integrate (2.2.2.2.35) with O(0) = 0: 


7 4Z 


= pa ws = — 
ie Gras = p, (1 e -) (Dp. 6h) t 


Q=pjil-e 
C2e222 4a) 
With this notation we have the following results: 


mee =< xX < Uu regime f= / 


TPs 


MaMEOuVeCc 2.202 521)—-(2.2.2..2.24) : 
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ake alee 

St Se LR Ci +. 3 Cox) (2 ee) 

. 2B aS 3 i = ete 
AXM 4 = hx + 3 Re C\x + 19 Cox (2a 

w =O + eae Cx (Q222 4) 

Y\ a 1 2 ORE oe 

1 ee I Se Wg: 

Yv = > * (52 C, +3 Cox) (23a 
b) U, SX < Uy: regime 4-7 

Ng = =] ; Mg = 1 


mmem@iceom (2.2.2.2.3/)-(2.2.2.2.41) 








se c aeies 2 3 
1 a pe a a D ca peep ibe 
n R R R 1 R R A 
ee XS sy x 
p BF /R 22/R 22/R 
lleeZ whee ag 1 Te Oc Z 4 
oS (3¢ aa Bi): sae 17 PR a D 2 = B) x 
C2222 25.2 330) 
47,2 ret 973 ie 
ora 6 A+ D 2 
2 3 
oo R - R xX + > (Cc - E) Nee 
Ae Ro A <7 : 
4) 97 3 a - Zz es ; 
ae be! a ee (Qe 4) 
b gee : ae z ¥ ae 
6 sai + 855) - £533 a5 ii - Q A + = D 
hxm, = ——————__-— + [ h + = J] x - ————_ x? 
‘ eZ AR Wig Le Whe Ae 
ei G- 422 8 xt e+ Lg 2G - 6 2B) a 
6 R2 24 R R 


(22 e222 





(wae 


yw = C - <2 Bx (epee cs) 
P48 22, - og 
vv = 377 + = (CoS = Bx) (202 ee 
e) Uy SX S PaGeCime so / 
a = =| “Sy My = aE 


paemepenrame(2.2.2.2.42)—-(2.2.2.2.46), we have: 


ny = 71 = iG ~ he 5 x? 3 = Kx") (2.299.208 
er ey ee AZo Oty = oe (2.2.2.2.45) 

Se See ae ee Ge nee ee: 

Sah = iC K + é QO G + Gx aR ae x” + is Kx 
(D.D 29 AG) 

° DUP 

YV = “Kx Cae (Zee ee) 
yw = K (x - 1) DDD DAs) 


epeeBoundary Matching 
Th 12 unk etme ee: 
ere are unknowns: Cy, Cyr Ar “RB 


47,2 Se 167 ° a0 m ae 
RZ IF S35 le EK AG Uy, Uys 


To determine these unknowns, we have the following 


eemarctions of continuity: 


Ti, ) Om “Se mall Sa eies ens fle tear 


bo As ea 


Bye COntInuity of n 


mons [{l/}. 


Pmcome mul cy Of yw because Of Ghee soul elon ee tees ners 





oe 


Strains are negligible [17], and also from the fact that 
yw must be continuous in Tt, otherwise the acceleration 


would be infinite, which is not possible because the load 


applied is finite. Since yw ewes em eyeienioial (ene Uy) and Uss and 
Uy and Us aie EUNCEVONS sO T,ete 1S Mecessary thee vw be 
continuous across xX = Uy and x = Us- 

em ecOncinulty of vv: Here vv must be continuous with 
respect to tT and thus must be continuous across x = Uy and 
xX = Us. 

weconeinuity of yw, yw! may be discontinuous. From 


fe yield surface and the direction of the strain rate 


weetor, this discontinuity can occur if |m = 1, which is 


4 | 
Mmeremchne case here since at xX = vy and. a= Us 1 4 | 


Therefore yw' must be continuous across x = Uy and x = Us. 


lutea << Fie 


Thus we have six conditions of continuity at each 
boundary resulting in 12 equations to determine the 12 un- 


Pmowns listed above. 


Solving and using the shallow shell approximations 


whenever possible, we have (Appendix Al): 


C = : De ORR BE 
SS = ee) ( ) 
= -e) - SARA Agate 
C= 2[p. (1 e -) pt] nee ) 
a = T ; 
k= -2[p(l =e) 7 pgt] De 50) 
24 ‘s hal as 
R B= 2ip,(1 -e ) = pat] QP 228) 


- a =e ee DD Pen ae 
Cy 2[p, (1 eo =) pt] ( ) 





does 

















47,2 a —, 
az ©, = Po(l - ee) - pot = 6ht CLO SS) 
oy COs 27,2 ye =—T Dee 2 
IN C(l + yun?) = -2Ipy(2 - ee") - pot) () +32 us*) 
(252 2 2) 
cr 8 Sarr. ee 2 we 
D) = =o Dee = = (Cu. c+ ICs U >) + 0(2> : 
(25-2522 woe) 
ea 472 ee ee 2 7,4 ee ee 
Ey = 2? Oe = nyl6(Cu5+ 2Cuju,)u; ee ever Q)u; + 5Cut}u, 
(2 22 ee 5G) 
bE; ee 3 
z Sg Sa Ola 
iL Re R3 ae aa, * , 
= mem Se 6(6(CUL + 2CULU,)U 
~ 4(3C - Q)u2 + 5Cu*}u? O22.2, 37) 
ae Bec 
ee 72 
G= h + 0 (S5) (22 ioe) 


Uy and Us awe fo "se sawn iIn Nee sEiee ce 


epee oo 


S 2 
5 c uy zu (222. 2e5o) 


It 


, : ae ; ; 
6 (Cu, + 2Cu ou >) (ay + Us) a (ate Q) (a, + usu, + Us ) 


A 2 2 = ( Ty DEC 
+ 5C (uy) + Uy ) (ay +. Us) 0 DE ) 





80. 


Thais system is equivalent to a first order ordinary 


Seererential equation, 


Bee scolved numerically. 


mer inal Results 


fatten che unknowns thus determi 


U5, we have the following results: 


ao 5 x < u,: regime Peaee (Gismilce isi 
ie ny = 
My = ] 
6h - (p - Pp.) 
Sh SO a), 
= ]- ans 
a én 


yw = 2[p (2 - ee") - 


Being hi giv non lne at aie me 


(Sample output in Appendix A3) 


and 


ned, - i seers yon Uy 


4a,b) 
~} 
el =p (2a ce 
= \p ~/% oZe4.4. 
pep 
ee oe: 
an oo (2 een) 


pat! (2, ox): (2522 2S) 


<7 — 
- 22 6éht-[p, (1l-e eee ; 2[p, 4 a) pt] ; 
i Rk 2 = 3 
CU PO PASTE) 
b) u, < x < Uy: regime 4-7 (dame nse ile JS) 
Ng = Sly Me = 1 
Z 2 Zz Z 
are > (p-p.) ~6h 3 ujtusu,tus uy, 
p Ieee oe 3 u,tu, u tu, 
=) Ge) aoe 
_ —p [x* (u,° = U5°)x* + u, U5] OOD OS Sep 

















of ae 


eee sOme Manipilations, mM, Cam be wreittem inet 


d 
zOrm ° 
(X-029) (ae x) 
eZ ll 2 > 
SE ee a rece cai ieee ena = —_ = s 
ae i ER : (p Pp.) > edt aS: P. EP Alo) 


- (p-p.) (u,tu,)] X= [2(p-p. = ede |) 


mee P—p.) (usrus)] a1 2 OA SAMA 65) 0)), 
PvP s Ii 2 u 3 er Se Sa 


yru 
S = same as (2.2.2.2.61) 


mM, = same as (2.2.2.2.62) 


fomemcame as (2.2,2,2.63) 


° Bags 


_ £4 ee oe eee oe 2 
io - 5 Ip, e ) pot] (> aa aD Us ) (225 252m) 
c) UU, <x < l: regime 5-7 
ag = -l] , My = 1 
= ITE Coa ce AA a 2: 


Se= Same as (2.2.2.2.61) 
Smee sais i 2.2, 2.2, 602) 
Yw = same as (2.2.2.2.63) 


- 42 oe oe Z 
YV = = i p,@ - e ) pot] > aed ll - (227 ae oo. 


With these results, it has been verified that all the 


Semcimelons Of COntinuity are satisfied. 





Seas 
teomrave also Verriled that the equations of equrl-— 


[Meron are all satisfied within the limits of the shallow 


Sa@elil approximation. (equations 2.1.2..-2.1.2.3.) 
meome (2.2.2.2, 22) which is 


oe 3 a 


” C 1 


and with C ana 0 from (22. 2.22425 lane F022 ee ace 


—T 

; Sip. (1 -e )- pot - ile 2 ae : 
25 ane era ee ney ee ae S u 
0 ee oe 


(2 5-2 eon 2 re 


Let: 


—T 
p,@ fey) Spas 2st 


f(t) = 2 en eee 


—T 
po (1 ee p.t 


Iemma Wrote -(2..2.2.2.09)0 im Ene comm: 





Do iy 2 5 Ee a > 
ee goal QE 3 
1 
inet: 
1 
¢= {2 ee ce (2.2.2.2.71) 
1 
Waene we Gderive: 
LS Cr 
Uy = Threat (Qo, DD5ID, 
ie LS! Goats ew aoa o> 





S22 
mmaeehe solution of the system (2.2.2.2.59) and 


Mn2.2.2.60) 1S equivalent to the determination of the 


mmc tion g. 


From (2.2.2.2.73) we can have u. in terms of f, f, ei g. 


2 


p@estltcuting into (2.2.2.2.60) the expressions of Uy, Us in 


obtained by differentia- 


u 


me? .2.72), (2.2.2.2.73) and a. 


mon of u.,f(t) and therefore f(t) being known, we arrive to 


a. 


the following equation to determine g: 


[p, (1 ~e@ ‘) - Dee 1.2ht] [p, (1 Se by gel (Se - 5)g = 


oo + 5) x { (<5) (poe ” = Pp.) [p, (1 -e') - Pot ~ 
~ 1.2h1] (llg* + 5)/g - (<5) (poe ite Pea a 2 [p, a - 
-—) - p.tl(g? + g + 1) (6g? + 10)/{g(gtl)]- 1.2hp [1 - 


~ (l+t)e ‘jg (2 ep ona 


To determine the initial value Sy Of the auxllany vation 
1able g, we assume that the applied pressure p is main- 


tained to the constant value Py ELOM Ome Toe and is allowed 


meomaecrease exponentially for tT > Wee As Teltuiten sds @ote 0 Eee eee To! 
g will have a constant value Jo and therefore g = 0 and 
= 0 for 0 Ss. Us BRCOM 2 62 o02 ee ee) oe nee p=p.e | and 


p. (1-e") replaced by Py and pot respectively, g= 0, and 
= 
Heeang that the last term -1.2hp, [1 - (1 + t)e (‘J comes from 
f which is zero in this case, we have the following equations 


to determine oe 





84. 


Milgy 5 


Jk 
Ta 'Py P,) (Pp 7p,~1. 2h) 3 
G’t95tl) (692 +10) 
Go fgot1) 


bo} 
© 


which can be reduced to: 


2g say osiN)) N(fote ys eal) ~ 3(g +1) (l1lg *+5) aa) 


ieee coCarranging: 


_ eo ae n = 
(g 1) (12g, Og. 10g. - 5) = 0 


O 
(272 ee re) 
In the limit, when Te > QO, Chis Vargument “Tse s er 
fed, and the initial value of g is to be determined from 
Wee 2.2./4a). The solution oe or 1 corresponds to WU, = Us 
mmemivas been found to be unfit. Therefore the initial value 


ee mmene solution of: 


1 Ja. - 9g * - 10g. - 5 = 0 (22 es 
Telit Nas one real root: 


g = 1.494 (22 ain) 


With the initial value Jo of g determined by equation 
: ' e e + 
Wee 2.2.75), the initial value To of g at tT = 0 from 
W2e2.2.2.74) will have the form - Using See ah rg meets 


mee to find the true limit of Jor we eb ranma. 


Ona. 


Vere 2 ‘% Be — 2 q 

3 ‘ Po ; (6g5+10) (3g 6g. got5g5+5) 

O 7 _ - Sie Gu. ce Vso : 
(D. p.) (p, PD. We Al, 432g- 414g) 4672g- 900g°-320g +50 





er a ee ee 


(25 Zea ele 


With the initial values of g and g thus determined, we 
@meenen integrate equation (2.2.2.2.74) numerically. With 
@ Known for each value of tT, we can determine the corre- 
Bpmeenaing values of u, and u. by using equations (2.2.2.2.72) 


1 2 
Pena 2.2.2.2.73). (Sample output in Appendix A3) 


Hemomocerve from the expression of £(1) am (22272-7770) 


felt: 
Fig Po(ire pet = Lighw 2 pp eer 
O S O S 
from which we have: 
De Oe tae Ze 
BOs ee re A C2e 22a eee 
awe i 0) =- T° 
De. Ok, eee 2 Ways 
Voy = LOL ee a (222 ee ee 
= 2 
Py Pe 10 + 69 5 


From these expressions of Yj and Ung, we can see Eons 


et * ae nd 
Yo and Ung are Desi bivenon ly er Daas IS Agvercberel cloves © iecLid 


again, by a different way, that for epee Sete 1.2h, the whole 


Shell collapses in one regime. 


We can also see that the value Ty OnE Niger Cse a hlgakellay- 
. 0 





cH opr 


meaetermined by: 


SEL be ee 
Pott e ) pt 


os 1.2ht, = 0 C22 2c ee 


ce tr = Ths Phe col Laps ee meciiie 5+) Or Rene hn eeiacm: 
5 io PPopweads to the whollershel 1 andythe ~f1 vs eect ade 


eeemOotlon ends. 


Pmmmoecond Stage of Motion, T. < t < Tt 


It ie 


Om t > 7 the whole shell collapse under regime 5-7 


ee 
as in the low pressure case. The expressions for the 
stresses: Dg S, 1s and for the acceleration Of the cen= 


tral point vie remain the same, because we have the same 
Meeieterons OL equilibrium (2.1.2.1)-(2.1.2.3) with the Same 
Known stresses: n, = -l, m, = 1, and the same boundary con- 


v 


qitions. 


Mmmsme nave. £rOm (2.24 2.15511) = (22a lls 
Ss Ex x?’ (2.2.2.1.11) 


as 24 (8 (LL ees ene @ 
1 + [3 (p p.)x 5 (p p.-1.2h)x*] 


R 
OPA eA) 


= 
ran 
| 


6h - (p-p.) ; P-P. 
i a a mt GI 


3 


(22 2 = ie ey 
pigeee&rom (2.2.2.1.7): 


YW = 2(p-p,) = 2(p,e 9 - PB.) 





Trem t low cules remain 


Se 


the same, wand from (2. 2. 25eo) 


eee. 2.1.10), we have: 
Yw = yw (1-x) 
gor aoe ae 
YV = RW Ox C1 ~) 
To determine YWo: we integrate YW from tT, to Tt: 
; : 1 T ae 
ist) = co Une | YW At = 2| (poe - p.) dt 
iL Bal 


yw (t) - ywo(t,) = 


Momeni 2.2.2.2. 


ywo(t,) = yw ot, ) == 


Therefore: 


e = ial 
yw ot) =: 2{p fe 


W 
YW 


With this, we 


yw = 


AZ aa! 
= “a Ip, f1-e a pv] pelle 25e) 


2Ip (el -~e') - pee 


63), we have: 


2Ip (1-e °}) De 


qa 


6 ae pets ts 7 + 2Ip.(1-e ' = oh 


(1) = 2[po(i-e Pop] 


have: 


2{p (1-e *) = "pth tee (DAD Ae ane) 


(a2 elon 


Sieh 


These expressions are exactly those obtained in the 


mit wor Low pressure (2.2.2.1.9) and (2.2.2.1.10). 


Heeecan COneluae that for 1 7 T the expressions 


alle 
Rumerie Stresses and velocities are the same as in the low 


pressure case. 


By the way, we note that tnese expressions are also 


maiese Corresponding to the interval u, < x < 1 an the first 


Ve 
Seage Of motion. Thus after spreading to the whole shell at 
i= tT the solution for this regime remains valid after- 


led 


wards, which is what we expect, and which ensures the contin- 


feiy Of the solution with respect to T. 


At T = Te Qete rina ned bi 


re - pot, = 0 (220.0) eon 


wie 1S the same as (2.2.2.1.14) we nave yw = vv = 0 and 


meme motion ceases. 


We can visualize the relative portions of Ty deter- 


mmea by (2.2.2.2.79) and tT, determined by (2.2.2.2.80) 


£ 
above in the following graph: 





ene) 
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Profile 
(a) | 
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Lh ee ey S| NE Se ee a a . 
el TE 
—T 
Raia 2 (alo) 
oar 
S 
Po 
| 
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of tT) & Te, | 
f c | 
(b) << | | 
| Z/aeN I i Pe te, Te ne, Nein 5% 
TY] We 
al tam ie ) 
tan (po) S 


tan ~ (p_+1.2h) 


eens Arter Bal 
Sercsamissibility of the Solution 





Ba ciemene 





m Kinematic Admissibility 0k. 


ayo = tT < T) 
Oe) XS Uy: regime E-7 


CeonavetOns £Om Kinematic admissibiii ty: 


©, < 0 (De 2ece ac 
C4 < 0 (AZ 2.82) 
ao (2.2,2.2.83) 


feeomer(2.1¢1. 1.) and with w and v from C2 2 2 ae 3) cana 


Wee. 2. 64) respectively, we have: wa 
Sls (Ses )eigee = 1a AG | 


wee 24 «4/8 = 
Yep = —— Creme! See ke p.t]x 5 


R 
(2 2 eee) 
2) 


pemee for tT < Tye paelem eer iyees i ee) en) 


(ise ") se oa POs 


ep will be negative if 


(2.2525 2.55) 


ae 2 lols), we have, simnlarly: 


=-[p, (l-e*) - p.t]x - 2[p. (1-e"*) sag 1.2nt) | 


OPA, Peete, 


e 


p 


will be negative if 


= 
Dp (l=e j= plotline 
x < oe Oe Oe oa in 2 csOre 


0 = 
Do ise) ewer 





oe 


Parca ll 2) with w and v from CP ere ON) ars hig t | 
(2.2.2.2.64) respectively, and using the shallow shell 


approximation, we have: 


eg a 
YK, = Ip, (1 e ) pt] (232525 2) Oe 
From this expression, we can see thatyK, 1S positive 
momeO < tT < T¢ amMGmcnuscee 2 e242 6c. oS) es Sat LsSt lcd awe Compa. 


memes. 2.2.2.85) and (2.2.2.2.87), and noticing that 
we xX << U,, we see Ghee 90 cl, CO maa: MOimomn 212 eee 1) er ctintal 


(2.2,.2.2.82) are simultaneously satisfied if we have: 


-T 
p_(l-e §) - pt - 1.2ht 

E ad (2.2.2.2. 89) 

po ti-e ) - pt 


Pa 

A numerical approach has been used for P,P. having 
Eeeevyalues 2.0, 4.0, 6.0 and for the ratio of (Z/R)/h having 
Mgemvalues 4.0, 8.0, 12.0, 16.0, 20.0 (which is equivalent 
to the ratio of depth/thickness having the values 1.0, 2.0, 
mermerea. 0, 5.0). It has been found that condition (2.2.2.2.= 
S94) 2S always satisfied; thus, for this range of Tt and x, 
the solution is kinematically admissible. (Appendix A3) 

a) 3 


< x <u regime 4-7 


1 Ds 


Goines ealenos! Meola S er ieleisieser eye iin Cletadlent | suey. 


S00 eee ee) 
Ko Ze ee Se) 
With w and v from (232260 Oo cate seo 2 OG) 


respectively, we have: 


ae 





ae 


0 Ea S x 


we have e. < 0. 


Since x S Uss 6 


With w and v from (2.2.2.2.63) and (2.2.2.2.66) 


respectively, and using the shallow shell approximation, we 


mmc Again: 


ee We 
Wp ace UO Vc 


Q ) - Pou! 


which is positive. 
MitG eli ehe imeerval (U,,U5), the solution is admis- 


Sible provided Uy and u. S20 S) se walieisl Svea elgvelie 3) << Uy, 


A numerical approach has been used to solve for Uy, and Us 


and it has been found that they do exist and are in that 


< ite 


weeder., (sample output in Appendix A 3) 


a2} Bigg Se S 1: regime 5-7 


The solution for this interval is the same as in 
the low pressure case. It has been proved that the solution 


is kinematically admissible for tT < Tes 


Therefore the solution is kinematically admissible for 


et < 1 


ey oe regime 5-7 


1 ig 
Lie wsOlUETON fOr eh le Stages On Ilo Mon iso Ne eSalle sdlc 
mice tor che low pressure case, TIhas solution nas been 


Suovedeto fc kinematically admissiole; LEhneretore tie 





Jo 


peiidecion 1S kinematically admissible for 0 < t < Te. 


ye Static NCU OL SS 1b Oa IL aL 97 


a) on 


4: 
Sone ee Vor vada Sater iiey se ng 0 
Cj tele AS Uy and Uy S XS ci a = ~] 
lo (Ge Uy SOS Uss we have 
Meee 4 2.2.2.2.65b) : 
Ce Gee Ca) 
ng = -1+ 32 = 4 px) D222, 52) 


where: 
ee (p-p_)x- — eee ea) Sie) eee 


il ao 
u tu 


_ [2(p-p. = 2h (p-p.) (uj, tu,) ] ; 


meea.'S a polynomial of 2nd degree OE eS ES ei 


epex) = Ax> = Cx - D 
Peon this expression Of a we see that a necessary 
Semaicion for ay to be admissible is that P(x) be positive. 


Paeco ice tne Gum Of the roots of 


ll 


Pix) O , we have 





aS, 


The numerical approach which computes u, and u 


I 2 
fieso evaluates P(u,), Plus), Beye! 
ma S 
(DD1) = ee a 3) 
- S 
‘CORI |) = os) = i 


It has been found that for the range of pressure: 
ee < eee. 6bh, and for the height of shell ranging from 
Meco 5 times its thickness, (DD1) and (DD2) have same signs 


and P(u,) > P (uy) = (0h. 


Since (DD1) and (DD2) are both positive or both nega- 


ave , Uy and u.~ are on the same side of S/2. We have thus 


2 


the following disposition, according to the sign of A= p-p.: 


ene) 12 eke) 





S) nae 
a a 
pase 2y ise |Site (Semele le: 
A = p-p. > O and P(x) = 0 A = Din Deans Grand P(x} = 20 


has imaginary roots has real roots 





mime in any case, P(x) > 0 for u Ss a 


A P(x) 





1 oe 
The maximum value P(u.) Of PG) iS settee wo coe aaa 


From the expression of ng i (2 02 eee 2) emer yier 


(x-u,) (u5-x) 


Ey 
Sa ae Sey ad 6 R os Vee eC) | ae 
Mie can be proved that; 
(x-u,) (u5-x) . 
x Dee WM = OD 


Since 0 < uy <r) 


and ee) de == JP (0) 











96. 
we have: 


7 2 
Gow fa: 3 h 


See = 


Z| 


Srunce Z/R and h are smaller than 1, we see that 


in.) im lee thus fn, 1S admissible for reqime 4-7. 
max a) 


b) a 


Semaitiqn fOr admissibility: 


OBI s << ek 


p 


m, has the same expression as in the case of low 


p 


pressure range, and it has been proved that M4 1s admis- 


Sepa le for: 


p -p_ < 6h (CSomaa fi Gite 2 ee 


and 12h +pe f£ - p. > 0 (condition (2.2.2.1.19)) 


S 


Sencvtion (2.2.22.),19)) has been founda mune uecIa serene 
satisfied for the range of pressure: oes oe 6h and for 
fmeeneight of the shell varying from 1 to 5 times its thick- 


Mess. (sample output in Appendix A3) 


Therefore m, 1S admissible for Seg oa Gia. 


g 


Ge) Ng? 


Soncwerone for admi SSioaliey. 


= 1h 


Kits 





ne a 





OF 


From the expression of Ng I 02.2377. 0 pea scan 


write Np mi senNemLorn. 


Ny = -] + BESS x S(x) (Dire 2 
R 
where 
S(x) = S(p-p_)x - 2(p-p. - 1.2h) — (2.2.2.2.94) 
3 Sc 2? S e e @ @ @ 
Hoe, 2 ~l, we must have 


0 


Ss) 


Wit e on = OF S(Xii is monotonically anereas ing 


memetron Of xX, and S(x) 1S positive for u, < x < 1 if 





2 
S(u,) Lee 
2h?) (2 See Ieee Sel wees. 2 1. oaths 6 
= 2 Z'Po Se IE kd 
(DA 2) 
The maximum of ae OCCULrS 4b 2 - =" 1) Vand esea ewe. 
athe es: 5 Oe Ouro 
n _- = ee =) = =) = = ———s 
a) aa x Lar eacnes P.) 5 (DP. ee | es Rn! 2 + 3h] 


Since Dag 6h, we have 


g a a 
(201 5) eae +2—(ht+3h) = -l] + 8he 


Since h and “ are small compared to 1, we see that, 


Wwe nave: 





Oi 


ae tt sh oe O, S(X) is monotonically decreasing 


[iimM@enlOn Of X and S{xX) 1s positive for u, < x < l if 





Z. 
peo) 1S: 
SO = Sie Ee = +t Diy) a0 
Be Pd 2 ‘PvP. ; 

or P-P. 

47 le, 2) 

6 

Ox P~P. tse, ee 0 


Siem ode cieases) With we Ehis second) Lions satis. Ted 


Maramnme< 1.3if it is at T = T 


i. L 


pe 1- py + 18h > 0 


fmmweoconmpare Enis CONndLE1OnN WiEn  Cconaieaon. (2.2225) le 


aniG MOnieece Ema T ¢ > Tt) ane re ties 


i - > poe f, we can see that it is always satisfied if 


fe 2.1.19) is. 


From the expression of n, (2.2.2.2.67), we have: 


6 
Beers se — Bless. < Ll. 2a) es 
8 R S S i 


‘ —— 
and n 9 O 


Slee. a omen) 
S 
1@) ao es 


n 8(p_-p) 


and Ng has an extremum whose value is: 
n,) = -l + ay aon wananiatialat: Was 
8° ext etal 6 


( 





pee 


HHhis extremum is larger than -1. For this extremum 


to be admissible, we must have also: 


(no) < 0 or 


ESE 


5) (ile Ziodae Saye) 
AS 9 S 
Rx x nna. pe “le ec, 


< Teler= Ria < i. “Enis Condieien Will ber scatvst ecw: wo ene ce 


5 (alee nee ero) 
oe x ge Zeal (2.20 oe) 


PRemmcOnlaGrttton (222.2... I9te we deduce: 


2 = as 
cage cans Parp< <a 


With this result, we have: 


57 ; 5(1.2h+p.-p) ’ eon 
R 6 -- R 





iMecetrore condition. (2.2.2.2 .96). will be sation ec mee 


22h x 


| 


jams "CONGItION-aSsS precisely <conGdwelom W242. can co coin 


Wemxenow that it is satisfied for a thin, shallow shell. 


The numerical method which gives uy and Us, has also 
mena that condition (AD2) defined in (2.2.2.2.95) is 


peti sfied. (Sample output in Appendix A3) 


Poets = fey 5 has the same expression, but now 


meno extremum at x = O. To be admissible this 








TOG 


extremum must be minimum or: 


(any) ar ZrO P-P. See An 
Since p decreases monotonically with T, this condition 


mee be satisfied for T > Ty teak cee Seechices e = Ty)? 


P-P. < PLe - p, < 1-2h (252.2520) 


Five e yo 


Relative Positions of cm and Ty 


The function po (1-e *) has a derivative egual to ae. 
From the law of the mean in mathematical analysis, 


there is a value T, Ores con Ua Se bel ehiever 


S| 
py (1 e ) 


il 





JC ee 


PGOlmn 2 sno 19)» we. have: 


Se ee 
Pisce =) (p, + 1.2h) 7, 


Therefore: 


n= 
Poe p. + IRE Wan, | 0) c Ty) 


: : prado [7 : 
Since TS Ty and since Poe decreases monotonically, 


we have: 


p e' < poe 3 nae sa OA Ztoh Oia; 


=e) _ 
Poe Des eo a 


and condition (2.2.2.2.97) is always satisfied. 


Thus, we have found that Np is admissible for 


) SS Tp. 


Seeonc lusion 


The solution is found to be kinematically and 
Statically admissible and is therefore an exact one accord- 


ing to the yield surface which has been selected. 


Peeeviaximum Central Deflection and Energy Absorbed. 
Since yw has the same expression as in the case of low 
pressure in both stages of motion (0 < T < T, and 


i S tS T ¢) the central deflection and the energy absorbed 





2 es 


will have also the same expressions: (the difference is 
that now PoP, Can go up to 6h) 


MaOMmm Ze 2226.62) a ) : 


- _i 2 
ims 2[(p, Po), 5P.T ¢ ] C22 cee ew 


Pee@imenez. 2 2 1. 26): 


27R°N Die 


2 _O | Se ee? _ =u 
mabs| ay qe SB Baht = Usrwse 98! | 


feeeiial Displacement Distribution. 
1. Normal Displacement Component. 


Beticeotrating (202.2252.63) wrth respect “co 1. eua 


observing that at tT = 0, yw = 0, the final normal displace- 
ment component yw, can be written as: 
= _ bie 21 ¢7- 
we 20 (pb, Po) Te 7 Pots ] (1l-x) (Jee ee) 


2. Tangential Displacement. 


The expression for the tangential displacement YV¢ 


depends on the position of x. 


meu and u are the initial values of Uy and uU., 


10 20 
respectively at T = 0, we have: 
Mer eo tes integrating (2.2.2.2:504) am respect to 1; 


20 


we have: 





IOeie 


aS. c el eae 
ea RIP. P.) Ty 5PoT¢ pre ata x) (2526252 29.9.8) 


b) YI as Ys? USING Bene velocity fronme( 2.2.2.7 .6>) mean 
(Zee OO ee re nave: 
TE 
ae x2 a! ; 2 2 
Ye = > I Ip, (1 e -) pt] (2x x U4 nea 
az (‘ft -T 
+>. Ip, (l-e = p,t]x(1-x) dt (en BE GS SS) 
Ux9 
where Th ES the. Value Of Geter waked uy = UE 
us (TL 4) = Ox 
TU Can be found numerically. 


Bes X < U 


Oe 


using thes velocity Erome (2.2.2 .2.049e 


Were? G6), and (2.2.2.2.68): 


we have: 
v, = 225 * 
YV ¢ R°}. 


,(6nt-[p (l-e ‘)-p.t]  2{p,(1l-e ")-p,t] 
O = O S 
{2 a ae a «far 


22 "2 = 2 2 
22 [p. (l-e )-p,t] (2x-x “Us, alae 
pa 
4z{'x1 iis 
+ aa ip, (1-e7*) -p_t] x(1-x) dt 
oz 


OP arse ey, 





104. 


Ty and Th5 are the values of t for which we have 


respectively: x = uy and x = u,: 


uy (tT 4) = x 


eet 


Z “2? 


These expressions of YV_e are rather complicated. 


Since YVe 1s not important compared to YWer approximately 


Oi 


eh 
YW ¢ R 


and Z/R is small for shallow shell, -we do not evaluate YVe 


numerically. 


ty 











LOSs 


peeec os iedlum High Pressure Range 6h. < Pe aes tee Ah 


Emme eat stage Of Motion, 0 < t < cae 


When eee ep 6h, then from the expression of me ein 
Wee. 2.2.62): 
(m AD eae = 0 
Depo. 6h 
See 3h ae 
foeeehn means that there 1S a Minimum at x = Q. Since at 
x = 0, le = 1, this also means that ne will increase from l1, 


which is not possible for rigid perfectly plastic material. 








106. 


We therefore assume these collapse regimes for the 


Shell: 


The solutions are 


eo =< Xx < U 


fgem (2.1.2.2): 


aie if LOM 


xXYW 


Il 


YW 


where 


ewe 2.3 )h 


(p ~ <2) x, 


[A 


ene ae elo. O defined 


yw 


< 
< 
< 


< 


then: 


regime E - 


(2 2a ee 2) oe 


Q 


regime 
regime 
regime 


regime 


Cee 


0 Mace oeen Cerined imeem. 


leak 
Bay 
4-7 


Dat 


ae wen. 


stiol, (02) 5 aera anes) 


4) 


2.47). 





LkOi7 


Zu SKS UU regime E-7 


O Te 

The equations of equilibrium and the flow rules are 
Biemcame as those corresponding Co the anterval 0 < x < Uys 
first stage of motion, of the medium pressure case. 


iiemesorestie resulting equatlons@rsthe s/Stem Mme berm. 


meer = XS 1S the same as in (2.2.2.2.9) which is: 


d*y dy O 
gee” & ae OY = pasa © 


But now, we have only one boundary condition: 
at x = ul, Y = Oy since Ss: =A0i-aG: x 5— Uo and s must be con- 
PeUOUS . 


The exact solution has been found (page 66) to be, with 


C mow GdliEferent from zero: 


2 
C, E E : ; 
i sre |-8 + oI, | Koa te SK, | to%| SE C,(éI,) + C.(&K,) 
_ 22 
g = zm x (2 2a) 


Expanding and using the shallow shell approximation, we 


obtain: 


Vos 





27,2 = > le: 4 
R4 O 


Using the boundary condition: 





¥(u.) Oe 
we obtain: : 
ee _ DF, oe 2 7 tL oe 3 
€5 = w Cus 3 COU Je 2s oe a) 





108. 
With this, we have then: 


a : 977% « 2 5 [we 3 : 
oo = De C, (x Un ) + z Cots Un VW HZ 2 cd) 


From equation (2.2.2.2.5) we have, with Y from 


(2.2.3.5): 
a on E E 
yy = ONG cu 422/R2 j oI. | K jdn- KG | Idan | 


22 
R 


ee 2. ee 
C, (61) = Rm © (EK) (2. 2 ee 


iu R 


Expanding and using shallow shell approximation, we 


obtain: 


AZ* » AZ ose 
Ra C,x ++ R27 C5, == lore Rx 





xyWw = Qx + Cox" + 


where Yp 1s the Euler's constant defined in page 70 


With from (2.2.2.3.6), we have: 
3 
oom Mm pe 5 Z 
Yw = 0 + Cj ll-3 Re — (Yp + FOG Rr x ]x 


Dye Z, 
[ee Re HUE + log a ] 








Since ° log = is small when = is small, and 


= 0. 5772 that lect 42? “o | milous) 
DI 
and = aad ee 


Gere 
R (y,tlog =) with respect to l. 


Tnerefore, we have: 











Oe 


= OT pz ©} + Cx We? 223.9) 


This expression can also be obtained by direct sub- 
Pmelevon Of: Y from (2.2.2.3.7) into(2.2;2-2.5), and which) 
is what we will do for yYv and m4 since they can be obtained 


mapalgebraic substitution or integration of Y. 


Byetmetedraring (22252.3.9) with resoeck "tou qeezem 


0 to tT, we have: 


w=Q+or-C, + Ca (22222103) 


iiemmnatial condition (yw) = 0 will be satisfied if we 


@ 
choose C, (0) = C.(0) : 0. 


From the equilibrium equation (2.2.2.2.2), we have: 


Diane 2) Ce eee 
Ree. Mes ee ee 


1S = 
(222 enn) 


Integrating with respect to Tt from 0 to T, noticing 


that Uo oma LUMCELON Of Tye we Ooballn. 


T Dae 
mee 2c | 22° : eee 22 Lg 2 
YV = za | 2 C,x + 3C x = | ge & + 3 C oY) Uy at | 
G22 S12) 
As before, the initial condition: +t = 0, yv = 0 will 


be satisfied if we choose S10) = C. (0) = Q. 





Paciimemosequd Libr lumeequation s( 2.0.2 a1)e wi times 


we obtain, after an integration from Uo EOux 1 le nese le 
O 


_ 2 72 ee 3 _ 2 3 
Seu = hx + 3 pz C1 3u. x + Zu, ) 
dh ae ih _ 3 & 
ae Fe Cox ale Xx a 3u, ) C25 Dees Send 3) 


Pye, SX SU 


1 eh Bute 4-7 


a) 


Here again, the equilibrium equations and the flow 
rules are exactly the same as those corresponding to the same 


interval, first stage of motion, of the medium pressure case. 


tmererore equations. (2.2.2.2.37) 7 (2.222.2.4)9) were 


Ea valida in this case and we have: 











. A = B DE Zs 
Yv = 5a7R + pr * {Cc - R Bx) Zee a) 
yw=C - 22 Bx (22.952. 38) 
Vag Ns le § Diag ty 
a Le ee © ee 
Sos 9 82°/R: 22/R 2 22/7R 
1 Wf oe 47° ee 3 ERG Q7 3 a Ze b 
Vy iy NIG ime eee eg DS be eB) 


C2 Dee 2) 





i De Re 











47,7 ae : 973 - 
Ree A +t az D Ae oe ee ae ; 
ee Ay 2 / Re Aye R2 
3 ee . 
ae (8 a D - 6 = yi See 
(2) 
Ye < oe 2 ae Sine 
Lom + 85D, - 22% a E- 0 Ao = D 
hxm = = 7 aE [ h + —-—— ] xX - -Orororo x 2 
¢ 16Z°/R* 4Z°/R* BASSE 
+i (Cc - 42% i ee ee L (8 Z° Deeee & B) x" 
6 R2 24 R32 


(292 5 22 le) 


4) a 1: regime 5-7 


Again, the equilibrium equations, the boundary con- 
ditions (at x = 1) and the flow rules are exactly the same 
as those corresponding to the same interval, first stage of 


motion, of the medium pressure case. 


Therefore equations (2.2.2.2.423-(2.2.2.2.46) are 


applicable here, | 


5) Boundary Matching 


2 7. eae e r 
A, pr 2B: oF e / 


There are 13 unknowns: le vk 2 


: ; e) ile 
a7 Sempehy - = 16% 5 . 
Meee’ R° 7! RR? *! Re Fr Gr 


pe: 





lNoeretnag that the conditions of continue sor on and s 
have been used in the determinations of the constants of 


Macegqration and the continuity of n, 1S not in question at 


p 


x = UL, we have the following conditions to determine these 


unknowns: 


GBs a 


eyeecemtinuity of N yr S, M4 ee Uy and x = u 


by conelnuity of v, W and w' at x = u. and x 


1 = Uy. 

The continuity of w' is required since at x = Uy, and 
xX = Us, jm | < 1 and since a discontinuity of w' is possible 
only when jm | = 6 

eecontinuity of vy andw atx = Uo > 

Here, at x = Use hy = 1 and a discontinuity of w' is 
allowed; thus we have 14 conditions of continuity: 6 at 
each boundary x = Uy and xX = U, and 2 at x = Uss resulting 


apparently in 14 equations for 13 unknowns, and the problem 
seems to be overspecified. Fortunately, the COnGlt ONS ere 
continuity of w and Rat x = UL can be satisfied by one 
Single equation and we have therefore 13 equations for 13 


unknowns. 


Solving and using the shallow shell approximation, we 


obtain the following results (Appendix A2) 





Cc = =e OD, SM 
<2 3 =€ CPE ORGS, 
K = -C (2 2.2). 32.06) 
C. =e (QD 2.3527) 


sz C, = C-Q (DAD as eis) 





PAIS 


a a 2 
Ngee ocuen = u,”) QD. 18) 


2(1-u,) ?(1+2u,)Q-12h 
Chet ee Oe (Qu) 2a ceo 
(l-u ye Eee ) 
O © 























Gerla : seas , 4 Sieh ee 
ees ac CaO) Uae = C=O) Une cone a 
Mee, 248" rae 3C-0 3 SG ee, cee 
pz fs (Cuzt2cuju,) uy cmcel 4+ Tat} + (C Q)unu, zulu, | 
(2 52. rae) 
ee _ Q7 3 ee _ 27,2 ee 5 °° 
DiS ay - BG = Sa (et + 2Ca a.) 
Die 7 CDOT . ee 4 4 ss 3. 8k 
Me -Taay | TC CO us + 5 (2G-O) = cae 
+) al 1 
C202 ee eo 
2 7 ore DU? PRE 2) Coton ae, 
E — 5 E-Q = A ee { > = 2( E04 
1 R (u, u,) [2 (uj, tu.) ] a Oo 
A ea ee 272 soy Qs 3 
+ 3 (2C Q) u- Cu] R2 nce Q)u- xC us] 
(222.2 6323) 
ia 162°: SZ 2 2 2 
+. — oe oe 
D 2 S a a. Sai 
AZ“/R 472/R (u. u,) [2 u, tu.) S 
- 2(C-0) u? oa 2 (0C-0) uo eeu" ] + 2 4 eb) ui-au"] 
O é O Oo R2°3 O O 


(Zee ee) 











ie 


ae _ ii er 2+ 2 1 C-26F 1h 
G=h D Be Q) 3Cu Ju. 212-(u,4u,)1 Se 
- 2(C-G)u2 + £(26-9)u? - Gu*] 
O eS) O O 
(ALE eae oF 
We observe that when Us QO, expressions 


Wee 2. 3. 22) -(2.2.2.3.25) will reduce to expressions 


Wee 235) —(2.2.2.2.58) of the medium pressure case, 


The equations for us are: 


G0 )*@eu ) = 2222" = op 
O O e 
(22 ee ous) 
ae 12h 7 
(cu ue) (1+u,) = p_-P_+6h (aE 0) 


(222 123 a ole» 


Meme@ean also solve (212.2:32326a) for Uo (Ade 
u, = 4-3 cos (TE) Don DT 


g = cos se 1 (ee 6G) On C2. 2p Oe) 


= 81h : 
» = cos? { ce (t = 0) 


(2 aoe OO 





| Hi ow 
0 and © have been defined in (2.2.2.2.35) and 


(2.2.2.2.47) respectively. 


With 
F(t) = pone eine 
Q 
12h 
F(0) = a ae 
Po p+6n 
and 
_ <n 2 
P(u.) =e u,) (1l+u.) 
equations (2.2.2.3.26a,b) become: 
Ee) = Ca) (22% 2g See 


end the variations of us can be visualized in the following 


graphs. 





Grapnic Determination of.u (tT) 


Ege, lt 





Teak eie 


(Nia) Soe Tor ue = QO, with Tt. determined by: 
E(t} el hic 
—T 
Pp, (1 e 0) - (Sa Se 6ht =) (eee ee) 


Heeeana Uu. are to be determined from: 


Ih 2 
3(2+u_)u* - 5u? - uv? 
is 24+ 2 28 Sh 
iL 
and 


ae 9 °° ne ee 9 > 
6 (Cugt2Cuy5u,) (u, +u,) 4(3Cc Q) (ay tuju,tus) 


7 Se SE eee 
+ SC (uj +u (uj +u5) + 12(C Q) ue 8Cu 0 


5) 
OS) 2) 

We can see from (2.2.2.3.26b) that if Pasa. = Gia, 
u(t=0) = 0 and thus find by another way the upper limit of 


the pressure difference P,P. for the medium pressure case. 


At t = 0, the initial value Us of Uo 1s given by 
fees. 2/) and (2.2.2.3.28b). Lt 18 worthy to note that 


hoieoa Given, Yoo depends only on the pressure difference 
os ee 

Tomumibeagrace (2.2.2.3531) and (2.2.253-32)) Gumericalks 
we note that (2.2.2.3.31) gives U4 when uy and u, are Hla aye 


Therefore, it would be convenient if we can solve 


Weems and (2.2.2.3.32) fos Che To do it, we 





ee) 








EAS ee 


Meererentiate (2.2.2.3.31) and using a Gon ele 2 elo ee 


Do 


After some algebraic manipulations, we obtain: 


ae an Ss A 
36hp, [1 Olea a ; (uy ua) Uy 
It =o) - 3 ae 3 
QO (1 u,) (1+3u,) [uct3 (2+u_) uy 10u; J 
: 1 


a 2 3 ae cP ey: 
(1 u,) (1+3u,) [ug+3 (2+u,) uy 10u,]0 


Es oe Z — 2..2 
Bae uo) (1+2u_)Q l2hjuju, 


2 
(us +uju, 


oe ilies 


, Pea +5) uy 
+ 2[(1-u,) (5+10u+3u_) Q- 36h] 


_ 2 = 2 a RF: 2 
x [2 (1 uo) (1+2u,)Q 12h} (uj tus) uy, 





G-uc 
3 a a) 2a neue = Gal 
Zhe. uo) (3+2u +u2)Q aes 2u.) | Sonn 
iL 
(232.2557 


Thus Uy) and U5 will bescomputed samen (252. 29 oor) ecu: 


meee. 3 -32a). TO determine Ehe anitiall values Yio and Yoo 
of Uy and U5, we use the same argument that has been used in 
the medium pressure case to determine the initial value of 


@f go. Then we find that u and u are to be computed 


S 10 20 


O 


from the following equations: 











Lise 


—_ ae 2 —_ 
oi) ee Uso? (1+2u. 5) 12h] uj us 50 


u 
_ 2 ye We YAO SOY DO 
+ 2[(1l Ugo) (otlout3us.) OQ. 7en I ptinn 10 


S, 


~ F512 (1-u 0) 2 (14+2u,.) 0 -12h] (u 


2 Zz 
00 Ui gt U59) U9 


Mew 
. 7 2 2 ee ¥ 90 10 _ 
Zale Yoo) (3+2u tue) Q. 12h (3 Moo) 1a) Fu, = 0 


and 


3 _..3 
3(2+u, 5) ae IUs 4 Uso Ys 


Su 


20 10 


Q is the value of 0 At TL ]—s07 cGerined. aa 


eee 2.2.35). 


With these initial values, the expression for ay 


will have the form 0/0 at T = 0. Using the L'Hospital's 


rule, to find the true limit aoe of uy at T = 0, we have 
found: 
Py Hg Yj 91449) 
Pio. 5 \iliaest eee 
Qo Or? 10” = 20 
- ee 2 : 
where with Sio = Uz gtU5 and Soo Ui 9 + Us) we have: 





a eo 


H(u (64+5u2 


= A 
(lt2u. t+3u ) 10 ) 


oo! M19 429! Be 2 La D6 


s 2 3 2 
2(S5tl3u.oti5u,. +3u,.) x (So +Uy QoQ) YYZ Q 


2 aes 2 2 2 2 
+ 2(345u5 t7uc. 3U5 6) Yi 9 + (1 Uso) yz Uso} 
3 2 
39 >» ,5 10 , ; “10 
Me 50) SF 104207 MiotT G7 (At 7U Gg t3uGg) (20g ti) 
20 20 
and 
eau u u..) = (14+3u__)%(14+3u__) f[llu,.u?. + See bie 
eo 10° 20 OO OO ee gee Ce 0) 
: 25 vs : ee 
F 1LOUT pYo9 7 ST Yo + be¥Got 4 QU 29 7 Mg! (S499 + FX 
2 
5 “10 2 3 2 2 
+ ra .. ~ 2(2+7u +3uc.) (4) pF3Uj Q¥ogt2U, QU5_, 
u 2u? +3u..u2.-5u- 
Fee ee eer eee 
Udo 3 OO ©O smo 


With these initial values, we can then go on to 
integrate numerically the system of equations made of 


Wee. >. 31) and (2.2.2.3.32a). (Appendix A4) 


6. Final Results 
With these unknowns determined, we have the following 
meSuLtsS: 


vigee a x < U regime E-H 


O° 





uz Or 


Peaonret 242.2, 544): 


fegom(2.2.2.3./)% 


a 2 ae 
: ee Uy) (1+2u+3u2)Q Aa er 
xS = —————— (x 6-0) 
2(l-u_)*(1+3u_) me 
Oo Oo 


: - 
- 2(1-u,) riiah “oes a3) 
3(1 uy) (1+3u,) 


(Zee 2 een 
REO (2.262.353. L3): 


(1-u,) ? (1+2u,+3u2)Q-12h 
“hxm, = hx + (x 2-3 x a 2us) 
6 (a) : (1+3u) 


2(l-u_)*(14+2u_)Q-12h 
= +(x" -4u3x ale 3u°) 
Was ol) 2 (leeesee 
(22s oo) 





alae 


ME@meGl. 2.25.10) : 


Q 





w= (1-x) DDD eos) 


feeemn (2 «2.2600 )2) : 











: 2x*3u_x*4+u? 
a 22 O O 
ih, ys) =) ———S (22.2.3 306ea0 
O 6X 
; : (x-u_) 7 (2x+u_) 
NS eee 0 Eee AP AGS) 
i Jil 
O Gx 
e) U, S$ X S$ Uy: regime 4-‘] 
Ny = i My = 1 


xn, : ees. 5Cu, u, 
3(u,t+u,) We 
; 2Cu-3(C-Q) . Aaa Z 2Cu-3(C-Q) ax 
3(uj+u,) O IY R 3 O 


oan 4 SF ee ee 
: 2 | QMuytusu,tus) Cu. (C-Q) 2 
R 


5C oa 2 

- ==(uftusL) + —----- u 

3(u,+u,) oe es 3 (u,+u,) 0 
12 ee 3 Se Z Lo 


(Zaz ae one a) 











ee 





cu = Aesey eee 
= er Aa _ O ; O 
a x + ix u,)u, x) 3 was 
1 ee 
_ — oi nae oes 2,6 (tu, ) (x+y) | 
u, Tu. 
CP Pris cree 15)) 
eet nally: 
2/R (x-u,) (u5-x) 


= —f{[ + — 
Ny 1 


3 x 
(1-u.) (1+3u.) . 


y) 2 s@ 
( : (l-u,) (34+2u tus) Q-12(3-2u,)h 
O 


3(u, tu.) 


Ge ae (> Ouro eS ein 
© © © 
= ee ff (ui tu.) xtu_u.) 
iL 2 LZ 
3(u,+u,) 


+ 2 [2(1-u,) ?(1+2u,) Q-12h] (xtu,) (x+u5) 


(22ers oe® 


xs, hxm, and yw remain the same as for u 


Y 


<< Se 


O IL 


Ewoment 2 ~242.7205 7) 2 





Wee 2c ee} 





d) 


and 


ic 


> 


Dies 


1236 


Uy SX 8 1: regime 5-7 


D = —], My = | 
Peeom (2.2.2.2. 44) : 
n 2 ees ee 30-09 eee oe 
9 R 3 2 2 3 O 


(222. Dries 


= 22 /R 4 = 2 aoe | S 
eile Tey (gl2 0244) (1+2u,) Q-12h] x 


eon 2 Dee 2 
ae: uo) (5+10u +3u-)0 SUSel) ox 


+} 


ov} 


tg) ARG) Oe 


(2522-2 oe 0)) 
XS, Tesi and yw remain unchanged. 
ieOnunde 12.2. 2542) 5 


* - 24 9 = 1 
ae T-0, x (1~x) Ze 2 6) ena) 





u Homgevem OY (2.252.3.27)) ane 1 Zope 2oc,)) 


uy and u, are SQOlutiens Of (222502 sree ome 


DERE 


* 





L24e 


it has been verified that the solutions obtained 


aiitsony the Conditions of contimuity at x = Uor X = Uy, 

x = Us and the boundary conditions. At x = Uas the slope 
vw 1s discontinuous and thus the curve x = u. is a hinge 
circle. It can be verified that across the hinge circle, 


the component Ky of the strain rate vector is positive, 
ema therefore satisfying the normality requirement, since 


aa — Ui, im, = l. 


O o 
It has been also verified that the solutions obtained 
feetsry the equations of equilibrium (2.1.2.1)—(2.1.2.3) 


within the limit of the shallow shell approximations. 


At T = Ge Ceceinmameds Oye 2eau. Oe eO us = O and the 


first stage ends. 


Pomme eond Stage Of Motion: 7T << t < 


O Ll 


At T = [ee a QO and the shell yields under 3 regimes 


aS in the case of medium pressure. 
When us = @> then: tf remr(2 2527 5.2) em ener 


ee os z 
Re 2 (Oe O11) ee aa aes 0(2) 


which is the same as (2.2.2.48} of the medium pressure case. 
Since in both medium and medium high pressure cases, the 
knowledge of C will determine otnee unknowns involving 
stresses and acceleration, then these unknowns will have the 


Same expressions and so are the stresses and accelerations. 











Ron 


To obtain the velocity, we have to integrate the 


meee lerations from ue 5 @ aes ae EG 00 meer 1) the accelera- 


iy 


tion yw comes from equation (2.2.2.2.5), which is: 


i> 


t t 
yw = Let p - Ba US) 4 p - ps + 6h 


Wiheieesm hiaGm2s2. 22.08). we have: 
yw = 2(p-p.) (1-x) 


Mmeecrating from i to tT, we have: 
e . e = —To_ —T - : . 
WC eayw te) 2Ip fle o) Se taiaio) TI le) 


where yw(t)) is the velocity at T = To: From the con- 


Pau cy OL yw with respect to t, we have: 
aie) = yw(t) 


and from (2.2.2.3.35),we have, with Ue cay) =e 


yw(t) = Q(t>) (1-x) 
ama, trom the definition of 0 J 1262252 oe: 


yw (TO) = Ip. (1-e" *9) - pot, + 6ht,] (1-x) 


Thus, we have: 


cE 


yw (T) = rom lacy Oo) - Date + 6ht_ J (1-x) 


~ Jip. te e287) De Ge Vp) 





a7 On 


Peaotmecguatt@mn(2.2.2.3.30), derining Tye we have: 


ena chus : 
; st == sae == : a 
yw = 2[p (l-e ) pt] (1-x) 
which is the same as in the medium pressure case. 


From the linearity of yw and trom GhevcConttmii ny yeu 


yw and yw! at x = uy and xX = Use 


pon for the whole shell as in the.medium pressure case. 


we obtain the same expres- 


For the meridional velocity, we observe from 


Wee. 2.3.14) that at Tor with u_ = 0, we have: 


and from a result just obtained above, 


re) - Q(t.) = 2tp.(1-e" °°) - Pavel 


From: 


C= 2(p-p.) fort? To: we have: 


C(t) - C(t.) = 2Ip (e o-e | 


) - oe Ga ee 


Mmempcom oe ceen that the conditionvos Comermurey Ot june 


WeMocmEmeamten respect to 1 is reduced to thac or C7) fhus 





Zee. 


C(t) = CAT 5) and we have: 


C(t) 


lI 


_ —To6 a = oval cia § - = 
2[p 1 e ) Pot, + pole e -) p(t EAM 


It 


r . -T 
eC) 2p. (l-e ) - pot] 
which is the same as (2.2.2.2.49) of the medium pressure 


SaSse. 


Since all the unknowns involving in the different ex- 
pressions of yV are all expressible in terms of C and Q, 
the former has been found to be the same and the latter 
memaros unchanged from its initial definition in (2.2.2,.2.47) 
we Can conclude that the meridional velocities for different 
intervals (0,u,), (U5,U5), (u5,1) have the same expressions 
as those corresponding to the medium pressure case. Thus 


(222.2.2.61)-(2.2.2.2.68) are applicable here. 


We observe also that all the stress expressions in 
this second stage can be obtained from the corresponding 
ones in the first stage by letting us = QO. Thus the stresses 
and accelerations are continuous with respect to T as ex- 


pected since the applied load is continuous. 


u, and U. are to be determined == rOmumlcecncecooo nau 


Wee 22 60) with the values of Uy Us at T= tT, as initial 


@onGitions . Sinema ct ae 





JE oye 


+ 
weave, trom (2.2.2.2.59), at t = he : 
yo Se 
By BU 
Woich is the same as (2.2.2.3.31) with us = Cree T= 1s 


PNILGIG) 2 En Bie Sj 23) aS aeeyeltvelee| ric) (PNA A eG i Rein Be 


make u. 0. Thus u,,u, are also continuous at tT = 7 


At tT = Ty determined by (2.2.2.2.79), we have u, = un=0, 


and the second stage ends. 


Geeebiaicad Stage of Motion. 


After T = Ty the whole shell collapses in one 
regime 5-7 as in the low pressure case. The solution is 
then the same as the second stage of motion of the medium 


pressure case. 


PNET, De = -]l, M, = ty and. (2522 Se ee cic e: 


still valid. 


Dv Admissibility of the Solution 





1. Kinematic Admissibility. 


< < 
a) 0 es STs 


Om << eUl regime E-H 


O° 


For kinematic admissibility in the interval, 





eS 
we must have: 


ke > 0, k,> 0, @ 


: d — Q 
Mabe yw from (2.2.2.3.4) and yv 
EeOnimet 2 scl. 1s J:) * YE, mes, 


wesc t < Te, we have e, < 0. 
From (2.1.1.3), similarly, we have: 
Yey = y(v'-y"w) = : £4 em © 
agit Zeal... 2) 8 we have: 


; : te A 
Tee a wy es 0 


pnb 7 trom (2.1.7.4); 


we 
tt 


-h(witty"v) ! = 0 


il) us SX < Uy: regime E-7 


For kinematic admissibility in this interval, we 


must have: 


(De 
AN 
© 
ve 
|v 
© 


aE 


Wenmemoel.1.1)-(2.1.1.4) aqgaanymeand waiehey 7 amc. 


meen (2.2.2.3.36a) and (2.2.2.3.35) respectively, we have: 


- 22 Q 1 


Yen = = reer} 
s) R 1 Uy 6X 





ee : 2 3 
[ 8x 3(2+u_) x “- Us ] 





LEE 


To have Ye, fee MUSE have: 


P(x) = 8x °-3 (2+u,) x*+ue out. (u = oa 
(SD 2 De ore 


It can be shown that to have P(x) negative in the 
interval (Uord,), it is necessary and sufficient to have 


either P(u) or P(u,) negative, depending on the case. 


From (2 2 we have: P(u) = ou cele) Wil Ga 


1s negative since 0 < us 1, and 


- a Di aes 
P(u,) = Buy 3(2+u)uj+us 


Pemngut te aerinitionof we im (2.2.2..35.51 ane ede 


2 


some transformations, we have: 


a oe? 
P(u,) — 3u, (U5 uy) 


fiers negative since 0 < U, < Us. 


Therefore ve is negative. 


0 


For Ye ys we have, from (2.2.2.3.35) and (2.2.2.3.36a): 


° Zane 


= lL 3 eu 3 
Yeg = zr I-0, rates (nese 3(2+tu,) x un! 


To have Ye, < 0, we must have: 


= 7 2_..3 
See = 10x” 3 (2+u) x us < 0 CO <P GES u,) 


(Qe 2a eo.) 





ies A 


Peomectore, 16 Can be shown that for O(x) to be 
negative in the interval (Uo,U,), it 1S necessary and 


Sufficient to have either Q(u.) Or Q(u,) negative, depending 
emmchne Case, 
FeeOmer@cnzZ. 215s 45), we have: 


— 2 — 
Q(u,) 6us (1 un) 0 
and 
= 3 tees 
Q(u,) = 10u, 3(2tu)uy Uo 
It has been found numerically that Q(u,) is negative. 


Therefore Ye, is negative. (Appendix A4) 


igienq Kgs we have, from (2.1.1.2) and with vw and vv 
feet 2.,2.2.3.35) and (2.2.2.3.36a) and using the shallow 


Shell approximation: 





Poe =Sisaec 0 > 0 and Un ek 


With rk, ae Og Yeq <0), Ye, < 0, the solution is kine- 
matically admissible in the interval (uj,u,). 


111) U, S$ X < Us: regime 4-7 


For Kinematic admissibility im Ehis anterval; we 


must have eg < O and ky 2? 0. 





Ren 


As before, with the shallow shell approximation, 


memnave, from (2.2.2.3.35) and CD Ses 








- 4 Be we i 
YK oo I-a, which is positive and 
; ey’: 
a ee See 
Yeg = R [-u, = which is negative 


ae) eS es Se regime 5-7: 
For kinematic admissibility we must have: 


Ko ? 0, e, < 0 


With the shallow shell approximation, we have: 





ky * 


which has been seen to be positive. 





For ey we haven Ero meme mew os caiici( 222 25 oemlnee 
s Wa 7e 
Yeo ~ ~ R T-uy * 


Pueeteis negative. Thus for 0 < Tt < ae except for con- 
dition Q(u,) mei we heave been tao) emuemy eri aide eer 


that the solution is kinematically admissible. 





ee 


meer. 5S T< T 


O 1h 


For this second stage of motion, it has been found 
from pages 81-84 that the solution is kinematically admis- 
feeeme if Condition (2.2.2.2.89) is satisfied. This has 
M@eem Verified numerically and found that the solution is 


kinematically admissible. (Appendix A4) 


at. < T Su 


aL ean 


For this third stage of motion, the solution has been 
mone, tO be kinematically admissible for 0 < T < Te where 
Te 1s determined by (2.2.2.1.14) and which is the instant 


where the motion ceases. 


The relative position of T can be seen in 


of ‘1’ te 
the following graph: 


(p.+6h) T (p.+1.2h) tT p_T 





Tr 


tan (p_+1.2h) 


Pelative Positions of To! Ty Te 
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feet atic Admissibility 


K< < 
a) 0 Se 


am 


p 


Iie Gerda cions fOr) admissib1 | ty. atsom Oni a ame 


® 


om 0 < xX < Use MM, = l and is admissible. 


p 


For uy < x < 1, we have from (2.2.2.3.13), 


(x-u_)? 


(2. 242050659 


fomen1s £Oorm, we find: 


mae) =sO Sand in ae ee 


oO oO 


Thus Ms has an extremum at x = US. To be admissible, it is 


necessary that this extremum be a maximum or mM, Uo) << 10), 


From the expression of m,, we have: 


p 


m'(u_) = [(1-u,)C-9] DD D2 OS 


O 


Sr T= 


From the expression of C in (2.2.2.3.20) we have: 


7 mo 
ms (Uo) 7 (TESS) (alee |) (O (QT-Q) 


With the expressionsof Q from (2.2.2.2.35) and of 0 


maou 2.252.2.4/7), we have: 





lease 


Let us consider the function in the square brackets: 


eit) = "wee Slee °)-. We have df/dt = -te ‘. Thus iy am 


1s negative for t >0 and f(t) is monotonically decreasing 


Hemet >-O. Since we have £(0) = 0, we conclude that £(71) is 
negative for tT>0Q. Therefore OT-O 1s negative for T> 0 and 
ma (u,) = 12/[(1-u,) ? (14+3u,) Q] x (OT-0) is negative as re- 
quired, and my (Ug) = 1 is a maximum. 

Or fe to be admissible, we must also have: 0 Ss <a 
BO u ae eh hOMeEne Gxpmesonen lor 7 Da (2 22 ee 
we have: 


ae x Po_| 36x44 [38a 5-4 (G-9) J (x? 4x2) } 


t 
> ee 


(2.2 225346) 
We can write M4 in - tEhe f£omm: 


a 
= = O oe 


; - 2 2 
™m C(x u,) (3x +2uxtus) 


d 12hx2 
- 7 2 ; ; 
+ 4[0-(1 us)C] (x tucxtus) 


Let us consider the expression inside the curly 
meekets., in the second term, we will £Eind (= (lene! 
Beettivye lf we Compare it with (2.2.2.3.45) and note that 


we have found ma (ug) £ 0 fand bhe Second wesc tor aS els 


1S positive. 





Leh ee 


in the first term, the factors involving x are 
mesative for x > Us and the coefficient C ae been found 
numerically to be positive, except for a few values of T 
closed to ie for which t= 0 and for a peak pressure close 


to the limit value Po. determined by 


i 
= — ie 
Por (D. in ae 


But these few negative values of C are also close to 
zero so that the expression inside the curly brackets 1s 


always positive. (Appendix A4) 


Therefore m! is negative for 0 < 7 < Ts and Ue Seas il 


p 


and a decreases monotonically from 1 to 0 as x varies from 


Uy two t!. Thus m, is admissible. 


p 


i) ng 


ime Conditions for admass ibiieee ee se No << 10% 





mor 0 < x < Uys and Us < x < l, we have Ny = ~l 
For Uy SUS ES Uy, we Have Erome(2a2.72580 7 oe 
(r=) ieee) 
Z 1 2 , 
Ng =~] + . = P (x) with: 
Py oe 2g eo 
ISS a Cx? + IT C(u,+u,) 3 tex 

= 0 ee on Per us 

uf cps oD) a5 ‘atu, 7 ale Ae ON 





Poe 
and the necessary condition is that P(x) be positive 
Ox U, S$ ¥ < Us. 
As in the case of medium pressure, we have found 


numerically that Uy and u, are both on the same side of the 


midpoint S/2 where S is the sum of the root of the equation 
Oe, 
and that 0 < P(u,) < P(u,) 


P (x) P(x) 


a= 


Fig. 13 a,b,c repeated 





Iss) e 
The maximum value of P(u.) is about 6h. An analysis 


analogous to page 95 allowsus to write 


Z Z 
ae) St ai Px) | 5 ees 


n 
d° max 


From this, we can see that n, is always negative, and 


> 


therefore Ny is admissible. 
Hee) No 
The conditions for HOV IL SHS) Lia I ahieNe reste) =| all Ny S Or 
Hoa 0 Roe 5 U5, Ng = -l and is admissible. 


ore Uo S&S 1, we have from (2.2.2.3.40a): 


22 


a = -—-] + ro R(x) , with 
BS oe HOA Liouony 2 2 ; 
Rae Cx aan + S[(Cc QO) x Cu lus 


It has been found that R(x) is always positive for 


u, Ses ca and its maximum value is about 6h so that we 
have: 
-l © Ny Ss Om ois Us <x< 1 
and therefore Ng is admissible. (Appendix A4) 
b) Eo, SSNS 
18) eae 
p 
Hiesconaltions for admisstbvlity ve. ss 7 < 0. 


The discussion is the same as in the case of medium 





WS 


pressure and the numerical results obtained are the same, 


> 


as far as the signs are concerned, that is, (DD1) 


u,-S/2 
and (DD2) = u,-S/2 have the same signs, 0 < P(u,) < Pus), 
and the maximum value of P(u.) is of the order of 10h 


so that (Appendix A4): 


(x-u,) (U5-x) 


cZ 
ae R Nears sar ag a ‘ ee 


oy} Un 


(Ny) ax = 


Pitch 1s 


NO 


5 


Hele si on 2 


AN 


Since h is of the order of a few per cent and Z/R is 


smaller than 1, we see that (n is still negative and 


nee 


the solution is admissible. 


et) m 


p 
MicecOndlte1Ons £TOr aGmnsoaion iis jc meee wa ole 
It has been found that nN is admissible (pages 51 
eameryo2) if: 
P-P. < 6h (AM1) and 
Zhe Sh) mes 0 (AM2) 
Since p = pe ‘decreases MOnOcOnlea ly conc le ron 
(AM1) is satisfied for tT > ie busca oh ames hcl ersvic Bc cin — Te 


Thus, we must have: 
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p, (l-e 





Graphic Determination of Ea 


(Decoy eye) 27: 


where a mo Gite BOOtwO @ COM ation hia 2 2 or nee 


As before, we note that the function po (l-e7") is 


; ; ; a 
continuous whose derivative 1s Poe From the theorem of 


the mean, we have a value Ta Ot Fe tiie yater aL (0,7) 


such that: 


p,(l-e ‘°) 


zo 





141, 


Therefore: 


ernie IB = 
Poe m = (p.+6h) 


a ; 
Slice T. < i and Poe decreases monotonically, we 


have: 
p eg 8O < p+ 6h Oia. 
O S : 


—-T 
ea OS a 6h 


pee (aM?) is satisfied.., 


For condition (AM2), again since p decreases. monotoni- 
cally, we see that it is satisfied: for Ga Ae eee if it is 


at Tes Or: 


(AM2) = 12h + a > 0 


5 


where Te is the instant when the motion ceases and is 


Beeesmined by (2.2.2.1.14). 


A numerical approach combining equation (2.2.2.1.14) 


with equation 


will allow us to find T ¢ and the upper limit Por of Py fox 
Piamem 4AM2) = 0. The value of Por. depends on the depth of 


mae Shell. (Graph G.9) 





ran. 
Aa.) Ny 


From the study in the medium pressure case in 


Beces 96-97 we can conclude that Ng is admissible if con- 


Seeerom (AD2) in (2.2.2.2.95) is satisfied. This has been 


@one numerically and n, has been found to be admissible. 


6 
Cc) TT) sus T 
For this: stage of motion, we have of = -l, mM, = 1 
and we have to examine m, and Ng 
From a result in page 10lwe have: 

6 We ae I oh 
O S ° 

Sememat nh. = -l has a minimum at x = Q. The other condi- 


0 


tions remain the same, and we can conclude that the solution 
is admissible in this stage of motion if condition (AM2) is 


Seeisried, 


Beeeonc lusion 


The solution has been found to be kinematically and 
statically admissible and is therefore an exact solution 


according to the yield surface which has been selected. 


BR. Final Displacement 
1. Normal Displacement 


<< 1 | ie ies eigte seve fene 
a) ee Sty where -_ e i" 





Eas 


at t = 0, we have: 


T 
tS 


eee: | ywdat + | ywdt 


0 ie 
O 


Peo e 2S. OS) ance (2.2.2.2,.63),.we have 





1 
an @) e 
= _ Q = =u es ae Ps 
YWe (1-x) | | I-a, dt + ae Pot; 2(p, P. 6h)7, 
O 
- p_ (t2-1*) ; (2225 2,3047 a 
S ie O a e a © 
meeoes X < U 
— —- 00 
Hea ws tne Value Of € fo. White teel | ee, 
ox om Ox 
then we have: 
“ox : ie ; "¢ ; 
ge. = | Ywdt + | YwaT + | YwdT 
0 oe oa 


From the expressions of Yw une (2.25.2 50e a 


ere? 3.35) and (2.2.2.2.63), we have: 


ie ic 


OX, 0 0 
Y We = | Qdt + (1l-x) | | I-u, at + 2(Do-Po) Te 


c 
S OX 


ee? 
-2(p.-p.-6h) tT, - p(t; Te) : 


Fe Aes re 4 9)), 
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2. Tangential Comoonent 


Let Yio = u, (0), Uog = u. (0) and meee ie and a 


determined by: 


Ue ‘TO? 7 
Ce ee ae 
Uy (To) ics 


If all these quantities exist for a point x, we will 


: < < 
have: i ee To. 
< & 
a) Yog SX S il 
Then, none of the quantities T T iG exists. We 
Ox ee 
have: 
ts . Te. 
Ve = | y ViGloa yYvdt 
O 45 
O 


fleoin (2. 2,2.,3.41) and (2.2.252.04)) we dave. 


T 
a Gh se) "ODay + ipe=p.) a 
We 58 I-u, See Fe 


O 
- -p_-6h - Sp_ (12-1?) 
ie Pg — 2” s TE ive 


(Qo oe a 


Igy) aa ears 


io eo} 20 


ie ti ses such Enat To, < Tor then we have: 





aS. 








a aL T 
2x , O ; i 
es yvdt + VV GT ee ares ie 
; tox ie 
PGOMONGZ 624253699), (202,283.41) ana (2.2.2 765 eve 
have: 
t i 
oa ry O ° Bs 
i = | | — eae x cla ne me so ee) Ca i 
Oo O 
O oe 
47 - _ iz wn ae 2 2 
= Rel x) L(p, al be (Py Ps 6h) 7, Woe cae 


mie Li x is 


Ws 
Eeaom (2.2 
have: 
aE 
mi 4 
ie  R 


O 


ane 
1 


Te ae aa 
O 


(22 Ceo 4 OO) 


Suemy tiat cP , we have: 


. Oe 


a a5 i 


O @ 2X @ f 6 
VMACUn Se | ty ViGl aiemct: | yvdt 


e) 


1 


c ax 


e) 


$253.39), (25.2.2.2.66 aman 2.2 eo 


Dn “ox cna 9 3 
Lay | ip, (l-e “)-pgt] (-u5t2x-x*) dt 
z 
O 
+ aes lige : eee ae _ (Vee) 
R 2 Se oe Lf Pot ox o 2x 


(1/2) p_(tE-12x) ] 


(Qe  seacc) 
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2) Conclusion 


From these examples, we can see that for other values 
@eex SUCH that x S Uy: the cases to be studied become 
numerous with the coming of the value of Ty) into the 
problem. This lengthy analysis will be tedious, without 
much interest, since ae all these cases, a numerical 
approacn has to be used to compute the actual values. 


Moreover, the tangential displacement is of the order 


0(Z/R) compared to the normal displacement. 


The most important displacement is that of the apex 
of the shell, which has been determined in the following 


Pieemand has been Carried out numerically. 


FF. Maximum Central Deflection and Total Energy Absorbed. 
1. Maximum Central Deflection 
The maximum central deflection YWo is given by: 
Te . 
Y Wee | yw dt 
O 


From (2,2.2.3.4) for 0 < tT < T, and, from the expres- 


sion of yw Ore ees T S$ Tp, we have: 
a —T = 
an = p,(l-e ) Die 6ht (ONp ca Tt) 


Micyeiuem page 126: 






Oe. 


yw 


- ecg are 
ae 2{p e -) Pt] Cc aa < Te) 


With these results, we have: 


_ _ on 2 de = 
a 2 1(p. Po)T, SPT] + (p.+6h) (l+57 ) Tt) ee 
(2 ae ees Sao) 

Maas result can also be obtained from (2.2.2.3. 46bD) 


mie xX = 0 and T = TT, 
Ox O 


2. Total Energy Absorbed 


From a result page 58, we have: 


Te, I 
be 2 , 
Be = 277R Ng | | pwxdx dt 
Oo Oo 
A J} Te 1 
- 2 
eye = 27R No | | | pwxdxaa =: | | pwxdxdtT J 
Oo Oo T. O 


0 


Or 0 5s) t < T,. we have, Prone 2 wee oe 


w= Q (O < x < uo) 


pia Eiom (2.2.2.3.35): 





- O - 
Yw = i-u, (1-x) (UO oa) 





ong ae Se Te, we Nave EEO 0 a7 a Oe 


yw = 2[pj(l-e") - pot] (1-x) (0 < x < 1) 


With these results, we have: 


a Uo i) 
abe = 27R°N. [ | pwxdx + | pwxdx ean 
O O "0 
f cl, 
4F | pwxdxdt | 
18 0 


ne 
which gives, after integration with respect to the space 


Warlable: 
2uR°N | *6 Py a = . 
2 | anes a - = _ 9 
Fabs 3Y I, 5) Ip, fe e ) (D. S)n}) ele J (1tu tus) dt 
‘f ae 5 -2T —T 
ai i, lees Ip. (e gS: )- pte A ley 
O 


Since Uo cannot be expressed explicitly in terms 
of tT, the first integral has to be evaluated numerically. 


The second integral can be Performed , and we obtain: 
i 
27R°N O 
O O = Ti 
= — = e -e 
Fabs 3Y | Z [Pa | 
O 


~2T) 


= 2 i 2 ae 
- (p.-6h) te ] (1tu tus) dt + 5 Be [({l-—e ) 


Melee oO) 2) = p_p 


aoe [(l+t Je “{l4t je °£] 


(Zier aero ) 
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earls 


Poe Dynamic Plastic Response of Simply Supported Shallow 
Spnerical Shells made of Rigid Plastic Material that 


teewe sence encoupled Square Yield Surtace 


Peo ieeocatl GC COllepso Prescures 


yt 
Lo =Le 





The yield regime for the whole shell is 4-5. 


imeoie the shell is Simply Supported, Enis ww) eesube 
in the same equations and the same boundary conditions as in 
Peceron 2.2.1 for simply supported shell. Therefore the 
Same results will apply to this case, and the collapse 
static pressure P. 1s 


Pare Om i 22 2... 2 7) for a shallow snell whose 


equation 1s: 





152); 
— Vi Fr 2 
y = R Oe , N >? 2: 


6n Z 
S n+i%l1R 





eomechat m, 1S allowed to be negative, the condition for 


? 


admissibility becomes: ( a Se Wit ele 1 veneer: 


ac 


n(n-2) 


L 
ae eo nal R < in way 
Se Beem (2.221428) page for a shallow shell of 
second degree: 
as a 
Paya 6h + 4 R 


2. if the shell is clamped, expression (2.2.1.12) 


page for a shallow shell of second degree becomes: 


a = i oe 
hxm, = hx @ (PD mn) x 


which yields: 


- yo eee 
ft 1 eh ( Rp? x 
At x= 1, Ms =-l] or: 
7 i 42 


With this expression of p, the stresses are: 


ny = —j, M6 = | 
s = ~-6Nnx 

no = -1 + 12h 2 x? 

m, = 1- 2x? 





The flow rules remain the same as in the simply 


Supported case, and we nave: 


Wo= wo ~- x) 
e 7 Dy e Z 
NE = 2 wox (1 x) 


From these results, we can verify that the solution is 
statically and kinematically admissible and is therefore an 


exact solution. 


¢ 


2.3.2 Dynamic Response 


Depending on the pressure difference Daa oe the 


yield regimes are: 


a) 4-5 for the whole shell, resulting in the same 
equations and the same results as in the low pressure case 


Seepoeeeton 2.2.2.1 and is valid for 0 < Paar = ee ine 
1! 38> 5a © 1c Uy < X < uu. and 


a5 £Or eee aia 1, resulting in the same equations and the 


Be@estesults as in section 2.2.2.2 and is valid for 1.2h < 


bye C-oF for 0 <x < u 


Po ~ Pg < 6-O0h. 


O 


c) C-E for 0 < xs uy, C-5 for u, Ss xs Uy, 


Beemeor UU. < x < u. and 4-5 for u, < x 1 resulting in 


1p = 2 2 


the same equations and the same results as in section 


O 
< 


Meee 2.3 and is valid for 6.0 h < Da, “eee Gop Wieleice tele 


to 


value of A depends on the depth of the shell and is not the 


Bemicmas tn Section 2.2.2.3; since the condition which me 
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has to satisfy is now different: in this case My 1S 
allowed to be negative. Therefore condition (2.2.2.1.19) 

which is: 


S- 
poe p, + 12h > 0 


is now not necessary and has to be substituted by the one 
which expresses that the minimum value of ie at xX = x be 
itmeger than -l. 


mrimom (2.2.2.1.18) page , which is 


_ 2 6h + (pg - Pp) 
Hes) sees 


we have: 
OF Digs 
4, ys, Mamie SES ey ae 
SU sia a 18h “m 
and we must have 
6h + pp. - 
nn a 
Leh m — 
6h + P. — © ; 
a mn = 
Since x < 1, this condition is satisfied if we have: 
6h + P, - P 
ith 


“> ae 
Poe Pp, + 30h > 0 


iiiiameondlEton Will be Ssatasi1ed forge a T > 1f it 


i_At Tt = Te. Therefore we must have: 





ie\oe 


P 


cae § 
f 
x p. + 30h > 0 Cee 


If we compare this condition with the one which 
Seeeeoses that ng be targer than -l at x = 1 in (2.2.2.1.26) 


ame Which 1S : 
pees & pn caaien Mc 
O S  — : 


femeadi See that condition (2.3.2.1) will be satisfied if 


eemertion (2.2.2.1.26) is. 


Thus condition (2.2.2.1.19) page 43 in the case of two- 
moment limited-interaction yield surface is replaced by 
Bemeition (2.2.2.1.26) for the case of uncoupled square 
yield surface. This later condition gives a limit value 
Por. OnE Py higher than that given by the former, and which 


can be found from the simultaneous system of equations: 


Gia ee Saye Pt, = 0 


Pon ie 


ee De an) abielal 


I 
© 


PoLS 
where ive and Por are the unknowns. 


This system has been solved numerically and the values 
Ep computed for shells whose depth varies from 1 to 5 
oL 


~rmes their thickness. 





ee. 


2.4 Dynamic Plastic Response of Simply Supported 


Shallow Shells made of rigid-perfectly plastic material 
Miceeooeys the uncoupled diamond yield surtace. 


wee otacic Collapse Pressures 





A. Simply Supported Edge 


The equilibrium equations are the same as (2.2.1.1)- 


Peril. 3) . 


The boundary conditions are also the same, which are, 


Hem a Simply supported edge: 


x = Ll: My = 0 - (2D dlee) 





As in section 2.2.1, the examination of the boundary 
conditions will allow us to narrow the choices on yield 
regimes and to select the correct one wnich has been found 


to be made of faces 3 and 5. 


For this yield regime, we have: 


* 





ae Sai! (2:aeelecsl) (Tab Vem 4310.3) 


Ms = | (Zeal (Table 2.1.3.3) 


J Eo a 


Pieenmec2 fli (20228. 3) ana (2.42121) (2.452 ave 


have the five equations needed to solve for the 5 unknowns 


stresses: 


S; Np, Ny Mas M4 « 


It has been found that: 


At x 


mrem this 


Teen this 


s s(p - <4) x (2. 4.1.3) 

_ - = z = (p- Ao) x? (2.4.1.4) 

: - Fi : z a2) x? (2.4.1.5) 

5 - a (p ; <2) x? | (2.4.1.6) 

bon Be fe = SEZ (Dy, Ava 

2 8h R 

= 1, we have Ms = 05 (Or 2 rOn. Be2 nore) 
O=5- oe (p - 

equation, we have the static collapse pressure 
p=p, = the (2.4.1.8) 


value of p, equations (2.4.1.3)-(2.4.1.7) become: 


Ss = -2hx 7 (2.4, 1.9) 
ng = - (1/2) + hex? Or ail) 
- _1L_ +4,? 
ny = 5 hex (2 ee ad 0) 
a 22 
a = oy (1-x~*) Oe a 
mail nee 
M 5) (l+x-) (24 1s) 





Dror. 


Poom Enese equations, and observing that h and Z/R 
feewomall and the product h(Z/R) is small compared to 1, 


memean see that the solution is statically admissible. 


fee crlow rules for regime 3-5 are: 


eg = 24 (2 274 5 eee) 
Ky = ky (2:54) ieee) 
With re Kgs e ys Ky from (2.1.1.1)-@.1.1.4), equations 


(2.4.14) and (2.41.15) yield the following ones in terms 


of velocity components: 

w= (v' - smh x (2.4.1.16 
w' +y"v = (w' + y"v)' x (254s) 

Beem ehnese Equations, and observing that for a shallow 


emer) OL second degree from’ (2.2.1.30), (2.2.1.31); 


y' = y"x = <2 x, we obtain: 
v = Ax (2) 2A aie), 
w= Bx? + w. (ae) 1. 18) 
Using the boundary condition of x = 1, we find A = 0 
and B = “Wo. JHeybks; = 
v = 0 (Do 20) 
w= wo (1-x?) (2a) 


With these velocity expressions and from (2.1.1.1)- 
Ceerierveanwicn vy wand vy" from.-(2.2.1.30) and. (2.2. lee) 


we find: 





Loe 


e = e = 7 22 ® 2 
e = aa = 5 we (1 x) 
K 9 = Sa = 2hw 


Since we Peo ieee gem a eens ly ,ewenmanic: 


one ce O; Kg = Ee Eas 


me ovseene Solurion as also kinematically admissible. 
Therefore the solution obtained is an exact one for the 


yield surface selected. 


B. Clamped Edge 


If the shell is clamped at the edge x = 1, the boundary 


conditions at x = 1 become: 
= 1, w= 0; w' = 0 or there will be a hinge 
Since a hinge is likely to develop at the clamped edge, 
we assume the following yield regimes: 
ime < X < UU: regime 3-5 


The system of equations for this regime is still 


pragma: (2.2.1.1)+@.2.1.3) and (2.4:1.1), (2/421.2). 


Therefore the results (2.4.1.3)-2.4.1.7) remain valid. 


Wee E£low rules (2.4.1.14), (2.4.1.15) still apply, and 


therefore equations (2.4.1.18), (2.4.1.19) remain valid. 
eee ex <  eegime 3=8 


For this regime, the system of equations to determine 


jnme stresses are (2.2.1.1)-(2.2.1.3), (2.4.1.1) and the 





Gor. 


equation corresponding to face 8 of the yield surface: 


My - m, = 1 (Table 218 3) (2 Ae) 


Since s, n Ny are eG De Vaeeermineauerom (2.2 aloe 


op’ 
wee. t.o), (2.4.1.1) as before, equations (2.4.1.3)- 


(2.4.1.5) still apply.To determine m, and Ma, we have 


Y 
(2.2.1.1) and (2.4.1.22). Solving this system with s from 


wee. i. 3), we obtain: 


= 2S 1 ey 2s. 2.2 
M4 = log . ah ( a) (x a) (ZG 2) 
ee Sn Ways ee Le (po = 27 
Q 7 J an ‘P ee (2A leery 
where me(u") SO eros ts =o ronmseiemmace lhe tc Ms must be con- 


mimnousS at x Wye 


The flow rule (2.4.1.14) still applies from which equa- 


imeem) (2.4.1.18) remains valid. 


Using the boundary condition at x = 1, we have: 


je 0) elgkel ee = 0, (Oe S15) Oe oo | 


The flow rule corresponding to face 8 of the yield sur- 


maeee 1s (Table 2.1.3.3): 


nS a5 Ka = 0 (24 se 26) 
With k, and ky from (2.1.1.2), (2.1.1.4), and with v= Ol, 
and the boundary condition: w= 0 at x = 1, we have: 


w= Clog x OO 57) 





tS 


3. Determination of the Unknowns 


There are 4 unknowns: B, C, u, p, for which we have 


4 equations: 


Meeconcinulty of M4 at x = u expressed by the equations: 


@apetrom (2.4.1.6) 


= 2 ae ot ee (2.4.1. 28) 


mine COnelnuity Of wat x = mG whleny Leon 2.4.4 oe 


gives: 


Bu? + We = C log u (224 te 


ims eCOntInuity of wi at x = Lyees Lace w' can be dis- 
a | = J, Wile as Pmnec. the weaceenere, 
Deeeherentiating (2.4.1.19) and (2.4.1.27) with respect to 


continuous only for |m 


x, we obtain: 


2Bu = = (2.4.1.30) 
ur 
Mme hie boundary conditvem at x9= 13) for w' to be dis- 
continuous at x = l, we must have nha = -l. Then from 
(2.4.1.23) we have: 
ie ioe eee 
Sees! ee an = eee) DALE Bi) 


Sonny ice eee) ane 2a lO) Oe es and (Cyc 


obtain: 





MGe 


W 
= -_ oO 
2w 
— = ee ee eerie @ BEE ae ae >. 


Sei wo cencany wt, Zo) and (2.4d,1e3)) for u andi, we 


Socain ; 
4h 22, 
NS Do Se eae (2.4.1.34) 
and 
i 
"oO Bom Weve; tl (2 deeb sey 
which yields: 
u = 0.64 approximately (2.4.1.436) 


With these values of B, C, p, the stresses and the 


WMemocicies become: 





Ss = -2h me (0 < 2 < fe) 2. ae eee) 
7 1 Z x? 
Le 5 —~ h Rae (nO eo (Dots o S18), 
i 1 Te 
Mp =- 5 thers ag x hil (2.4.1.39) 
_ 1 ut-x? 
ie ao (0 < x < u) (2.4.1.40a) 
oe (u SRS 1h 
2 xX xX -U ea 
m, = log = - ={S- (2.4.1.40b) 
1 mae 2 
.— >> ->— Wee = < Hh) (2.4.1. 41a) 
ks. tog a St < x < (2.4.1.41b) 
Mg = + log = - ne? (0 x < u) oe 





263. 


v= Q 
ae 7 1 x 2 (0 aS wee EL) (22 dees 2a) 
w= wo tl 1-2 log u u2 aa a: 
5 = Sy: (a <x <1) (DA Meo) 


oa Ogu -- - 


Seaamissibility 


It can be verified that the solution obtained is 
kinematically and statically admissible. Therefore it 1s 


an exact solution according toc the yield surface adopted. 


From the velocity expressions (2.4.1.42a,b), we have: 


i. = 0, and therefore the apex 1S a regular point, and 
20. . 
t Se ° , 1 = 
(w eG 127 lee a * 0: there is a hinge at x 1 


with Aw > 0 and K 4< 0 as required 
2.4.2 Dynamic Response of Simply Supported Shells 
Le <SO lhe sere 


If the pressure is not ‘too high, the yield regime 


will be the same as in the static problem: regime 3-5. 


The flow rules and the boundary conditions are the same 
as in the static problem. Therefore the velocity expres- 
sions are the same as in (2.4.1.20) and (2.4.1.21) but now 


Oe iowa Funct On Of selme. 


For regime 3-5, relations (2.4.1.1), (2.4.1.2) still 


apply. 





NGA. 
PaOCreeGuetT om OL Ouro (2.1.2.3), Using (2.42.1 1) tend 
ieee y and y from (2.2.1.30) and (2.2.1.31), we have, 
ome (2.4.1.21): . ; 
Nee 3 Doh 
co YW (X-x ) - (p - ao 


Integrating this equation from 0 to x, we have: 


Ne Ary (Sx2- 3 =x*) - = (p = 22). ; 


Z R C274. 20) 


Paemmegucdiren (2. 22) > and using relation (254. 17205 


we have: 
'4+2m,.) = h + 
h (xm¢ 2m 4 ) h XS 


Meme smorhOMm)(2.4.2.1)5 we have, after some mani pula— 


tions: 


el 2 = ee IL 3 e, 1 5 a ae 3 
h Apes m,) = hx + yw. (5x ries ) =(p-4 A) x 


Integrating this equation from 0 to x, and simplifying 


yee , we obtain: 


a ee ee Mr 2 
st, es ot ea 
(2A) 
Doingetmie: bolindary Condition at x — 1 fom simply 


Supported shell, we obtain: 


2 lene 1 ee. 
ely Se oye 3 4h) 


From the expression of ioe in (2.4.1.8), we have: 





ISD 


=p - 4h (AE, Ome, 8) 


Then: 


ir 3 
ee (PD) (2545225) 


With YW Own Vicema yo mmiroOm (2.4271) and C2 a 2s) 


p-p,~ 8h 3(p-p,) 
7s ger are x* - Sa em De Ao on 


eae irom (2.4.2.2) and using (2.4.2.3) we have: 


-p_-8h p- 
Ly Ps ,  PrPg 


| a Eo a eee 
‘ 5 [ Téh x Téh ey) (2 426) 


iWen. rom relariroms(2.4.1.2) ,twe have: 


Mo = 3 - —ee SCX? - ee "YW (24. 2a 


VEO Grenshbll le rgshibiecteatyshe Loy 4 «Ile Ze ak ete A= sage jie 
wre. 20), and using relation (2.4.1.1) and with y” from 
(2.2.1.31) and xs from (2.4.2.5) we have, after some mani- 


pulations: 


pomalclt 3 (O15) 
dd 2 ee 22 Es Sea > 5 
ee (x are = x + > [ 1 oe 3 en] 


Integrating this equation from 0 to x and simplifying 


bye x’, we obtain: 
ag 2728 ; (p-p.) 


feat tome aetna a eee Fc: AD i PO in 
Rr | ié6 = Te * 73 


ee wo) 


Then, from (2.4.1.1), we have: 





166. 
_ _ 4 _ 22 ,P=Ps-8h 2 PF Ps _, 
ao Se SS ig 6 Se 


Iie (2. AaeeZeee 


iacediecdtcing (2.4.2.4) with respect to 1, and using 


Meni nietal Condition: at t= 0, Se O, we have: 


yw = 2 Ip (1-e %) a (2 De do) 
With this result, the velocities become: 
vv = 0 
yw a 2 ioei=er ) = Dw) lass) .Aa2 i) 


2. Admissibility 
a) Kinematic Admissibility 


From the static problem, we have seen that the solu- 


tion is kinematically admissible for: 


ee ee 


eur rom (2.4.2.10) : 
p (l-e ') -pt>0 
O Ss ° 
oP er 


a ne determined by: 


a i See Oc (2: AnD ID) 


YW, = 0, then yw 


I 
-) 
~ 


and the motion ceases. 





ibis 3 
Dime atic Admissibility 


From the expressions (2.4.2.6)-(2.4.2.9), we can 


see that we can write the stresses in the following forms: 


m, = 1/2 + P(x) (2524 Zr) 
My ee) es | (2.4.2.14) 

22 
n, = - 1/2 + =—h P(x) (24S) 

d R 

22 
eg = I = ols Oleg (a 6 

where 
p-p,7 sh P-P. 

P(x) = —~=-— x? - == x! Cra ela) 


16h Aioyel 


The conditions for static admissibility of regime 3-5 


are: 
Ox ms < i) (2 es) 
0 <™, i (Qemte. 19 ) 
ab << 1M <40 (2. ¢42,.20)) 
“i 2 (2.4.2.21) 


A gt Tyr nN, in (2% Ata? 1S) = 


(2.4.2.16) and observing that(2Z/R)h 1s small compared to 


From the expressions of m,, m n 


Mere see that conditions (2.4.2.18)-(2.4.2.21) can alli be 


satisfied 1f£ we have: 


< P(x) <= (Cree 22) 


i 
DO| 





iene 


Paciimemeome Discos tONneOrer(x)jean (20452717). welmave- 


; x 
P'(x) = g-[(p-p_-8h) - 2(p-p_)x*] (Zee) 
Ea) — eet © 14 x) = 0, and 
x2 = PoP, 8 (222d) 
2 ~ “2(p=po) ae 


ages oh ae 8h, then p-p.~8h Ss eeu n cleGeimcda es 


with tT. We have 2 cases: 


=? 8 
For p = poe ae Pp, Or T < Mere Po/P. 1 Eloy) 183) oye 
hand side of (2.4.2.24) is negative and Xo is imaginary. 
From (2.4.2.23) we can see that P'(x) is negative: P(x) 
decreases monotonically from P(0) = 0 to P(1l) = -1/2 and 


weee.Z2.22) 1s satisfied. 


= at 4 
For p = Poe < Pp, or log Po/P S les de where T 


S- ie 


is the instant when the motion Ceases, andes Gehilnecean 


Wee 2.12), the right-hand side of (2.4.2.24) is positive 


7 1 
p-p,~ 8h : 
ag 2(p-p.) 


and we have the following tables of variation of P(x) 


and Xo exists: 


depending on the position of x 


with respect to l. 


2 








6) 


Preomesnese tables, we Gan see that 1£ x. < J, 


Z 
P(x) will have a minimum at x, which is smaller than ay 
and is not admissible. If Xo > 1, then P'(x) is negative 


in the interval (0,1) and P(x) decreases monotonically 
mom P{(0) = 0 to P(l) = -1/2, (2.4.2.22) is satisfied. 


Beice to have (2.4.2.22) satisfied, we must have: 


Kee Sill 


2 


emer com (2.4.2.24) 


Pp - P. toe el < 0 


Since p decreases monotonically, this condition is 
oes tied at 0 < Tt < Te 1: see eee = Tes and we have the 


Samadi tion: 


pe lf - ee 1 Slee (DEBS) 


Oo 


(oe = Das, 2h) 


ASnumerical approach has been used for PoP. varying 
from 0.4h to 8.0h and for shells whose deptn varies from 1 
to 5 times their thickness. It has been found that 


ree? . 25) is satisfied. 


ia) it 2) loa 8h, then we have 3 cases: 
(DC eo One 0 et Log Ip /(p.+8h) | 


ies arcqmesmend Side Or. (2.1402) 2o ee is pool ei ve and 


hex SctG and is smaller than 2. 





70. 
PeoOnMmercore ~pLeSS lon. Of Pix) in (2.4.2223).) we Have 


the following table of variation of P(x) 





+ 
“a = 


Thus P(x) has a maximum at x and “fem (25 4525 Loe 


7 
this maximum is: 


2 


P(x) = yée (p-p,-8h) - (p-p,) x2] 


es 





ea with Xo aohiete 4 2 2 4s 


(p-p.—8h) * 


P(x.) = 64h (p-p_) 


Beoms(2.4.2.22) to be satisfied, we must have: 


= ae 2 
(p Ps 8h) a1 


64h (p-p. ) 
or, Since P-P. ? Ove 
(p-p..) * - 16h (p-p.) + 64h? < 32h (p-p.) 
This condition will be satisfied if we have, with 
Sh < P-P.° 


p-p, < 16(1.5 + ¥2) h 
Since p decreases with T, this inequality is satisfied 


for T > 0 if it is at t = 0. Thus, we must have: 





I idk 


PaO S WGC ss ya) in C22 DRAG) 


(B)0 < p-p. < 8h or: log po /p.+8h < Tt < log p,/p,- 
For this interval, we have: 
Dee = 0 and p-p, - 8h < 0 


From (2.4.2.24), we can see that its right-hand side 
is negative; Xo 1S imaginary. From (2.4.2.23), we see that 
P'(x) is then negative and P(x) decreases monotonically 
trom P({0) = 0 to P(l) = -1/2, and condition (2.4.2.22) is 


Satisfied. 


(Y)p-p. <0 eee lieie: p/P. <t < Ty. 


In this case, the condition of admissibility is the same 
Bassthat of the second part of i), which iS condition 


ea. 2.25). 


mia) Conclusion. 


From the previous study, we conclude that the 
solution obtained is admissible for the lesser of the 2 
values of p. determined by (2.4.2.25) and (2.4.2.26) which 


are: 


=f _ ™ 
Poe p. + 8h > 0 (eee 5) 


P, 7 Py x 16(1.5 + Y2)h (2.4.2.26) 


O 
Condition (2.4.2.25) has been found to be much more 


heooctameravye than (2.4.2.26) . 





Wee, 


The maximum value Por Cf Po for condition (2.4.2.25) 
can be obtained by solving the simultaneous system of 


sqeetions made of (2.4.2.12) which is: 


ae ye = 
Pox, (1 Ss Pott 


| 
© 


(2542 eZ ) 


tl 
OQ 


and eo Tie p, + 8h (25 Ane Sy) 


PoL 


the unknown being Te and Po.’ 


Beeeihnal Displacement 


iMieegGmattnig (2.4.2.1) from 0 to Tes with the initial 


condition: yw = 0 at T = OQ, we have: 


SS ee yee 
YW ¢ sl(p=p.) tT, 5P.te] (1 oo) (Zeca) 


The displacement distribution is parabolic and there 
are no Singularities at x = 0 (no hinge). The displacement 


1S maximum at the central point x = 0 and its value is: 


— 38 _ 2 ee: 
Wig @ Gs se > 57a, Wie! 


4. Energy Absorbed 


As in section 2.2.2.1, the energy absorbed is given by: 


U¢ 1 
= 2 ; 
E obs = 27R Ne | | ovine cl Ge 


O O 


With p = Pn.e and w from ga are Eb ate en Ka loncrstibiol 





df). 





iN PS —TF¢E 7) ‘As 
= ey = = : ae 
ees Ole 7 1 a ) Qos (tT ele ! 


Comparing this expression with (2.2.2.1.28) 
in the low pressure case of the two-moment limited inter- 


action yield surface, we see that for the same values of 


Po and Po7 we have: 
ae: 
(E aps) _ = GF Faps) 
uncoupled two-moment 
square limited interaction 


The same conclusion applies to the maximum central 


deflection. 
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(b) Velocity and Displacement Distribution 


siress and Velocity, Displacement Distributions 
Shallow Shells - Uncoupled Diamond Yield Surface 


PA Cee 
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fe comparison of the Results 


Ao Oa tOMehen ol MLeca ane Uncoupled (soMdre 
Yield Surfaces 


The static collapse pressures are the same for 
both yield surfaces, in the case of simply supported shells. 
The solutions for the dynamic problem can be applied to 
both yield surfaces, with the only difference that the 
upper limit Por Of the peak value Po of the applied pressure 
is higher for the uncoupled square yield surface. In the 
case of the two-moment limited interaction yield surface, 
aes Limit Por is solution of a simultaneous system of 
equations made of (2e222 oles) and (2.2.2. 19) Ow reheecise, 


respectively: 


—T 
~ fy) - = 
Poi, : Ps f : 


Poe Pama Iezlal 


I 
oO 


In the case of the uncoupled square yield surface, this 
eaertva Po. is solution of a simultaneous system of 
Pemecrons made of (2.2.2.1.8) and (2.2.2.1.26) which are, 


respectively: 


(1-e ‘f) 


PoL ~ Pg £7 


—TFf 


Pore - Py + 18h = QO 


Graph G.9 presents these values of Por as runes Hom 


SMEemratio depth, eal ckness . 





oe 


2.5.2 Two-moment Limited and Uncoupled Diamond 
Hore ole races 


For the uncoupled diamond yield surface, it has 
been found that the maximum value of the pressure differ- 
ence P.-P. under which the whole shell collapses in one 
regime is much larger compared to that of the two-moment 
limited interaction yield surface. In fact, for some values 
of this pressure difference, we are in the medium high 


pressure range of the two-moment limited interaction yield 


surface. 


The meridional comoonent of the velocity is zero 
according to the uncoupled square yield surface whereas it 
is of the order 0(Z/R) compared to the normal component 
according to the two-moment limited interaction yield sur- 
face and is not negligible if we want to be consistent with 


the shallow shell approximations adopted. 


Another aspect of the problem is the expressions of the 
stress resultants. In the case of the uncoupled diamond 
yield surface, all the stress and moment resultants are 
expressible in terms of a function P(x) defined in (2.4.2.17) 
and the problem of static admissibility can be carried out 
Deeestucying only this function P(x). This fact simplifies 
the whole problem considerably, for we might have noted 
that the problem of determining the admissibility of the 
Solution 1s also as involved as the problem of obtaining a 


solution for the whole problem. 





do a 


One disadvantage of the uncoupled diamond yield 
pi@e@eaece, however, 1S that it is not as accurate as the two- 
moment limited interaction yield surface as suggested by 
imeection of the bounds in (2.1.3.1) and (2.1.3.3). 

Since we have here two yield surfaces, one is more accurate 
than the other, We have also used a method provosed by 
prem Hedge and B. Paul [13] to find a better approximation 
for the linearized yield surfaces: a peak value of the 
applied dynamic pressure is chosen, which is 26h in this 
case. Then the graphs of the maximum central deflection and 
Of the energy absorbed according to the two yield surfaces 
are concurrently drawn as functions of the depth of the 
shell aint varies from 0 to 1.25 times its thickness. 
According to the choice of the uncoupled diamond yield 
surface, which 1s considered here as the approximate yield 
surface, with respect to the two-moment limited interaction 
yield surface, which is considered here as the improved 
yield surface, we have 3 cases: the approximate yield sur- 
face is made to inscribe, circumscribe the improved yield 
surface or to have the same static collapse pressure as 
that of the improved yield surface. Each choice will give a 
curve of variation of the maximum central deflection and one 
of the energy absorbed. By inspecting the position of these 
curves relative to the corresponding ones given by the 
improved yield surface, we can see which alternative gives a 


Curve that is closer to the corresponding one of the 





ees 


improved yield surface and therefore decide which approx- 


mpeton 1S better, 


From the graphs G. 7,8 we found as in Ref. [13] 
that the approximation that makes the static collapse 
pressure equal to that of the imoroved yield surface gives 
the best result (except for the energy absorbed in the case 
of very shallow shell) and that the inscribed and circum- 
scribed approximate yield surface give the upper and lower 
Beinds respectively for the maximum central deflection and 


energy absorbed. 


Because of the simplicity of the solution obtained with 
the uncoupled diamond yield surface for sage whose 
reference surface 1S a quadric, it is natural that we should 
try to extend the value of the upper limit Por of the peak 
value of the applied pressure by allowing the shell to col- 
lapse under more than one regime. This attempt has been 
Samered out: at a pressure Py larger than Por determined by 
eme system made of (2.4.2.12) and (2.4.2.25a), 1t was as- 
sumed that the shell would collapse in one regime 3-5 


merge 6a,b) for 0 < T < Ty where Ty is determined by 


Pp_e 


5 - Py + 8h = 0 


ac tT = t,, the moment Ms will reach a minimum equal to zero 
at x = 1.) | 


POtat, 2 « Mi would become negative for u < x < 1. and 
{ eer 


cee 


the shell would collapse in more than one regime. We first 





Aye 
assume the following collapse regimes (Fig. 6a,b): 


OS ok 


[A 


us regime 3-5 
eee see regime 3-8 


ie ciMeilOnetGnimhas Dbeem found to be solutionmen eam 


Ordinary differential equation of second order. 


The initial values have been found to be: 


en 1 = Ty) u(t,) = 1 which is expected, and: 


° - Sh 
u(t, ) = > 0 


3[p,(1-e *1) - p.t] 


which is not possible since u should shrink from its in- 


ietral value u(t,) = 1-to a smaller value for v > Ty) and 


this will require u(t, ) to be negative. 
Then, another combination has been tried: 


Or x uy: regime 3-5 


N 


iat) u, < X < Us? regime alse & 


ra 


mad) wu 


{A 


> <x 1: regime 3-8 


but did not give any meaningful result. 
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Ss) (GENERAL SPHERICAL SHELES 


3.1 General Relations 


feel! Strain Rate-Velocity Relations 
For spherical shells loaded symmetrically, the rela- 
tions between the generalized strain rates and velocities 


are [17]: 





Fig. 18 
é = ¥ cot ¢-w (3.1 
Kk, = -h cot b(vew') (nena 
é4 = way (a ee 
K 4 =-h (vtw') (3.1.1.4) 
with ( )' = d{ )/do 


a. seebe ls 5) 
and h MO /AN . 
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In this case, the reference length A is the radius R 


Smee ecne shell: 





A = R 
moon trom (2.1.1.5) : 
ae 
n= EN 
O 
a 
lot = aR (3 lio) 


fee. 2 BQUilibrium Equations 


Using the results from [{17] for spherical shells loaded 
symmetrically, with some modifications to include the inertia 
terms, we have the following expressions for the rates of 


internal energy dissipation d. and of external energy 


1e: 
input Bee? 
a. 
dst = | (N pe, + N62, + My K g + suse Sin @ d®¢ 
O 
e se e 
do yet = | [(p_-yw) w S130) 0) se (py-v¥) v Sin @] sin @¢$ dq 
O 
# sin b(n gv + Sw - him, (vew') 1 
d=a 


where P, and Py are respectively the normal and meridional 


components of the applied load. 


With ep ey Kgs K from (3.1.1.1)-(3.1.1.4) and apply- 


taqethe Principle of virtual velocities, we obtain the 





oie 


following equations of equilibrium: 


S Sin ¢? = h{(mysin do)’ - Mg, COS dj] - AD sin db + zyh*f,sing 
(n Son. 1) N, Cos 6 -~- s sin @ + Py SOs =o ae ee 
(ni, + n 4) Sin @ + (s sin ¢$)' + Py sin @ = yw sin 6 
2 =v +w' is the rate of Slope change or angular velocity 


As in the case of shallow shells, we observe that 9% 


p 


has a sense of rotation CON EEALY CLO etal Or Os tan ye M4 and 


4 >: oe 
the term — yh’, = -1(-Q 
a, a 


moment of inertia of unit area of the shell with respect to 


), where i is the non-dimensional 


the reference surface, represents the rotary inertia. 


The term —2hp 4 represents a momentarising from the 
fact that the dimensional meridional component P 4 of the 
applied load P is acting at a distance H from the reference 


surface. 


If the load is applied only in the normal direction, 
which is our case, we have 2 0, p= Pos and if the 


rotary inertia is neglected, then we have 


S Sila @ = h[(m,sin o)' - M, COS 6] (Stele ie) 
(n 4 siLi0) ay) Ny, COS d6- s sin ¢ = yV Sin @¢ 
es a ae) 
(n,tn,) sin » + (s sin $)' + p Sin 6 = yW Sin $6 
at (BD 8) 
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For a shell without a hole, simply supported at the 
edge >? = @ and made of isotropic material, we have the 
following boundary conditions: 

At 0 = 0: 

From isotropy considerations: 


m coca 1 


YN) see 0 


S o é p 


From symmetry considerations: 


v= 0, w'= 0 or there will be a hinge 
he) = 0s 


For simply supported edge, not free to move: 


My = Q 


w = 0, v=0 or there will be a hinge 


The initial conditions are: 


At tT = Q, Vv =ve-=Q, W =w = Q 
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peeeDynamic Plastic Response of Simply Supported Spherical 
E@elis Made of Rigid-Perfectly Plastic Material that 


Obevs the Two-Moment Limited-Interaction Yield Surface 
SSS SS SS EEE 





As in the case of shallow shells, the applied dynamic 
pressure decays exponentially with time. Its analytic 


representation in terms of non-dimensional quantities is: 


where the peak pressure Po iommlancdete tian Erne State Col 


iktese pressure. 


Bro. ls Static Collapse Pressures 
Hodge has obtained an exact solution according to 
the two-moment limited-interaction yield surface for simply 


supported and clamped spherical caps [18]. 


The static collapse pressures are found to be: 
For simply supported spherical cap: 


The yield regime is 5-7 (Figs. 4a,b) and 


>) =? | eee eee WD Tn il) 
= 1 eS eee ae 
cs COs G4 


For clamped spherical cap: 


The yield regimes are: 0 < @ < 1: Shea fl , 
(Figs. 4a,b 


_ 2h sin $1 
oe Ms [= sin a (ome 42) 
exe: ~ cos $y = Sa o7 





liejeye 


where 7 PGmoGlUL1lOn. Of; 


Sin 4 1 + log (sin a/sin 1) 
Lee eying = ~ 


ae NOG = (COs Ohm eos EO) 
log ae o Sin 4 i 


(3 02 aliens) 


These results are valid for spherical caps whose aper- 


ture angle a satisfies the condition: 


2? ® 
eos8 Slits es Sabo Sta al So 
sin? a cos a -- 
(Sie2 sc leee 
Which is required so that the value n,(a) of n, at 6 =a 


remains negative as required by the yield before profile 
selected. 
Pmeec Dynamic Response of Simply Supported Spherical 

Shells 

Sankaranarayanan has obtained an exact solution 
according to the two-moment limited-interaction yield sur- 
face for a simply supported spherical cap under dynamic 
loading [14]. The yield regime is 5-7 (Figs. 4a,b) and the 


results are: 


a = -l, Me — ay il 
Lp Paap. *: HCO SA IL 
> aw: eae Wo Singdcos log (sec a + tan a) * 
progisec ¢ + tan 6) | 0) ) } 
sin . sin ¢ cos @ 


C322 52am) 





er. 


eee (p-2)log (seco+tan>) — p 
hon a} ie Ising | 


: =i (tee SEES cern) 


p 
yw 2| x secxdx 
= {lodsecs+ tan?) - | 


hsing log (seco+ tana) log (seca+tana) 
see 2a.) 
we = — 1 + preran?¢ + Mo cos*¢[l - HOGS: Umea 
) 2 cos“ ¢ log (secattana) 
[l - yet] 
log (seca+tana) (eZee 


= [log(sec attan a)][log(sec at+ttan a)-sin oa] (p-ps) 


of 
ier es attana )]* - 2 x sec xX oe 
O 
See ers) 


: : DP uerelisSce eee cst @) 


On aor ee log(sec a + tan o Peete Zod) 


2 


log(sec ¢ + tan 9) 


log(sec a + tan qa) (3.2.2.6) 


YV = YW. elsat) 0) [cl = 


YW is the velocity of the apex and is given by: 


[log(secattana)] [log (secat+tana) -sina] [p.d-e *)-p.t] 


re OL 
2 [log(secattana)]* - 2 x sec x dx 
O 
(Cpt 7 oe) 
Mie sOlULTOnmis eenemaetca lly sacimile>s  olemrer “Oy = Te 
Were. T 


¢ Tomene “oSstant ween tneenetlonme-aces ane. 1s 
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determined by: 


p, (1-e *f) = Pate ie (Soe Zee) 


Hite OCak Value p Of the dynamic toad is such that: 
O 


me 2 + Gin sim aOog. (SeCGun + caimia) - 
log (secat+tana) [log(secattana)-3sina] + | xsecxdx 
so that N y= =] be minimum at ¢ = 0 i 


(SoA ee 8)) 
Mien, le ican be verified numerically that the solution 


is statically admissible. 


The final displacement component distributions are 


given by: 


Noea(sice tan ¢) 
log(sec a + tan qa) 


=P 
++ 


Wits Sry cise or, {Pa ] Be Pere elds) 


; DL LOG (See tte Gane? 
ywo sin od [1 loeteae wo am a. Gor. 2 oe) 


ll 


ae 
Where Ke is the finalWeentral displacement: 


log (seca+tanal[log (seca+tana) -sina] [(pop.) tT, - $P.T 27) 


lo. “5 al 
2 [log (secattana) = 2| x sec xax| 
0 
(Sere ee 
where T ¢ ovmegiven oy (3.2.2.3)% 
The energy absorbed is: 
1 oi aU > 
= —— 27 dt Ba eee 
E bs | | Paz 27TR“sin oO dd cya ara) ee) 
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Wien the non-dimensional quantities : 


tet ye 2 Gee 
7! R’' PW 
O O 
we have: 
Te. 0 . 
E abs = 2T7R°N | | pw sin > dgodt 
O O 


Pee (S.2e2e5) and (3226207) 701 Com be Snownm thaceve 


have: 


n - rR2y [eteye (sere ela mmshelen || ivere| (Vere ee eeu re) esate fey) . 
O 


abs OL 
21 [10g (sec attan a)]* - 2 x sec xax| 
O 


{3 Ds es sree) joke fab (lt) e TEI} (3.2.2:514) 
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fees Limited Series Expansion of the Solution Of the 
Dynamic Problem 


To have a better idea on the range of validity of the 
dynamic solution and on the static admissibility, we have - 
made a limited series expansion of the solution obtained in 
Wee, 2.1)-(3.2.2.7), by neglecting terms of 4th order or 
higher. 


ies olution 











With: 
i 5 7 
’ a x >.< 
= _-_ — + = 
a ee 7 Uy) 8 Sic 
zZ by 6 
x x x 
es oe Dar Tie 





ie Sr 





3 5 7 
x > 61x 
+ = —_—— —_— 
log (sec x tan of ‘Jo ce on =OA0 
3 5 7 
log (sec x + tan x) - Sin x = ae ae ee Ss 





x 





2 y 6 
a x xX 
log (sec u + tan u) du = 5 + 54 + 1a4 
O 
x 
ae FS? Blinc 


| usec udu = x ths + re90 + Ipssao 


0 


meemcsing the results from 3.2.1:1, it can be 


Shown that we obtain: 


P, = 2 + Fas se 1s ) (322 3b) 





¢ in the denominator of the second term 


The presence of @ 
comes from the fact that in this part, we have chosen the 
shell radius as the reference length A in the non-dimen- 
sionalizing procedure. If the base radius B of the shell is 


chosen as the reference length instead, then from: 


B = R Sin a, we have: 


| | 
= Sln @ 
Pe Pe 
ble eee 
Sd ee 


Then for small value of qa, where gq? << 1, it can be 
Shown that: 


Z 
' 
Po - * g t Gh 


pemeaacreeen ftound for Shallow sAells ion (2-22) 228). 





POO. 


The limited series expansion of the dynamic solution is: 


_— a2-o? — 4( 2-0?) go? | Pee, a 1562 eae Aer 
My = Ge > Bo? a 

















(282 3) 
=e _ ee, Ze 
n, = - 1+ 3h aay ne + SB 
+ oS g2y- 5 4 2 76? 4 O16 _ 3207 , 87), 
yn! 5 12 She eS 3 
(sa2eons) 
a i — Ce ee 2.3.4 
Aa 5 e 2 ) J (Si y2mmore) 
ey oe ee ee ee (B55 825) 
Lars 6 re) 6 eS 
oo olny eee eS eee (3.2.3.6) 
es 30 Po Pot Age sea 


2. Admissibility 
a) Kinematic 


The conditions for kinematic admissibility for regime 


Peeeifigs. 4a,b) are e4 < 0 and ky > 0. 
For ey: Oe econ Sei.) eae) (3,9.9.0)- -73.2.55.5) ¢ 
; b 2 o? 
25> "Wy git 9 aa 


From (3.2.3.6) and since a* is smaller than 1, we can 
see that Wo is positive for Tw < Tes where T is determined 


by (3.2.2.8). Therefore for Ft to be negative, we must nave: 





TLS) 


ee, with > = a: 


R 
A 
Oo] & 


which is satisfied since o* is smaller than 1. 


e 





For K gs we Dave £RGM (Selle). anes 2, 5.4) aoe ern 
Ww 9 242 
re: _O 3 A _— & v) 
Kg = h cot o m (1 + ad R R ) 


Row Kg to be positive, with ¢ < a < 1, we must have: 


2 a2 2 
6 





- ao? <1 
With o* << 1, this inequality is always satisfied. 


Therefore the solution is kinematically admissible. 





moe Static 
i) Ny: the conditions for admissibility of Ng are: 
-Lsny < 0 
Prom (3.2.5.3), we Raves 
Ng = > X PN(o) (32220 eo) 
with: 
= HNIOO 2 
20 P Ps ., 8h _ S 40-160 
RC = oe eee ot La ad = —(p-p,)9 
30-8a 7 25-7a* 
+ h err ase ( =) S (Sates eo) 
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peo = OF we have: 


Ng an and 
Np =<] 
Mmewertore a necessary Condition for Np to be admissible 
is that there be a minimum at ¢ = 0, or 
(np b=0 a 
Garson 5.25526), (See woo 
' 30-8a? 25-70? 
Sh “T5a2 ~ (P-P,) “gr 2 0 


Solving for p-p., we obtain approximately, with ou cele 


1 ae es 


Since p 1S monotonically decreasing with 1, thas condi- 


mem 1S Satisfied at T > 0 1f it is at t = 0 or: 
Ae 2 
Poe < 2 eta ce (S525 S510) 


If the radius B is used as reference length and if 
fame-- 1, then we find again a condition for static admis— 


Sibility for the low pressure case of shallow shells. 


Moreover, for P-P. > 0, it can be proved that for Par ee 
satisfying (3.2.3.10) the coefficients of $°, 6%, ¢ and the 
constant term in PN(¢) (3.2.3.9) are all positive. There- 


fore PN($¢) is positive for see 0. 


For p-p. < O, at can be venrliled  nunecrveaI iy seat aon 


eee Secisevyaeng (Sa2-3-10) and wi cdwOm wom We stOOk, 1/6, 
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hy 


Bence we must have @ << 1), PN($) is still positive. 


Therefore Ng increases monotonically with 9, and for 


Ny to be admissible, it is necessary and sufficient that at 





>? = a, we have: 
Np a) “ (Ng 
ie ower al 270° 
S100 oman iu emrrmmen = oy eras aeeyren) «Sie vi{ 32 2 eSraelele) 





With Dar Pe limited by (3.2.3.10) and with h small for 
Mme Shell, condition (3.2.3.11) is always satisfied and 


meeorerore Ny Po woemics) > lew 


ie) re the Conditions “HOLT pacts s Dodie 7s Os ee are: 


Cr ce 
a Om ax 


From the expression of a Teh os eo ne Arties 2) a celal alts Wet 


mi = —g x PM(Q) (epee, 
with: 
DD Oey 16 17(p- 
PM($) = ne pe Niece = p19? 


_ : j=l = 2 
p De SoMa” . Hoke s 30-144 
Lon fe) 1542 6h iS) 





(32S) 
At » = 0, we have: 
m = 0 and Mi. =e 
0) 0) 


Therefore a necessary condition for a to be admissible 


is that there be a maximum at ¢ = 0 or: 





6, 


M9 < 0 
mo) 6 
PeOtims 2 sl). we have: 


PP. 30-1407 30-807 
(m 4) Ee ee ee 
d=0 6h 15 15a? 


which yields, approximately, with a’ << 1: 
p~p 


6h a? 
g ge (bee 


tro meonoimlon wil berSatus ete anos GeO melee mld 


6h ong 
Deh. ace eee (3c Ors eeA)) 


Semparing (3.2.3.14) with "Wse2. 3.10pm we seer that 


memes. 14) 1s Satisfied if (3.2.3.10)"is: 


Since ite must be positive and since m, = 0 at ?d=a, 


Siechier necessary condition is that: 


sae em re ae Sa 6 
4 | )o< 
which yields: 
p-p 2 2 
S o} 30+8a 
Tee! (esse = © 


Since p decreases with T, ‘this ineguality is satisfied 








oie tT < Te 1f it is at tT = Te or: 
p e fa 2 2 
Oa ee Se ee ae SUG) 
6h 30 sie — 


For p. Satiseying (3.223-20)), teeniocmoec man a mines 


mumerically that this condition is satistied., However, 





Meee. 14) and (3.2.3.15) are only necessary conditions: 
From the expression of m ihm 3 8 2 ao eo eae m4 =e) 
fom.) = O and for 


AMO) 10 


where PM(6}) is defined in (3.250. 53) <r Since Ee = One, 
have 3 roots and if none or 2 of them are in the interval 
oro) , m, may still be positive in this interval. A 


numerical approach is necessary to find the variations of 


m, in the interval (0,a). This has not been carried out, 


> 


Since it has been done in [14] for the general case. 


Heels OLfscome Interest Ee mete that a2: the applved 


Meaeis a rectangular pulse, we have: 


tI 


p=p, ford<t<l 


and p 0 fom ieee 


£ 


emeeror this load, condition (3.2.3.15) becomes: 


Pp 2 2 
S a 30+8a 
TGA Scenes) ntsc 


with P. from (3.2.3.1), this condition becomes, with 


a* << 1, approximately: 


2 
a? (1 - “2-) < 3h (323.18) 
With a* < 1, this inequality will be satisfied for 


ieee, or with h<<1, a will be small.” Then af B is the 


base radius of the shell, we have, approximately: a = B/R. 





Then, we must have: 


fog hs 
[A 
Ud 
x 
| 


5 
ae | 


O 
Py 
Zales 


For thin shell, 2H/B is small, and therefore this 


inequality can be satisfied only by a very shallow shell. 
We can also write this equality in the form: 
B? 5S = x 2H x R 


We have: 


Bo = 2 (2R=%) 
Or, Since the shell is shallow, 
Bee ae 
Raemuacomed clone se 2. Slo maiks equivalent iO: 
Z S 2 x 2H 
or the height of the shell must be smaller than 3/8 its 


thickness. This result has previously been found in the 


case of shallow shell. 
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(b) Velocity and Displacement Distribution 


Biewess ‘and Velocity, Displacement Distributions 


Deep Shell - Two-Moment Limited Interaction Yield Surface 
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Sau} penane Plastic Response of Simply Supported 


meuewlcal shells made of Rigiud—Perfecely (Plaseue 
eee tal that Obeys the Uncoupled Square Yieldeoureace 





Pe beset CeCol lapse Pressures 


Mice teld regquine 156° °4—5 of Pigs. “(5a,5) for both 


Simply supported and clamped edge. 


In the case of simply supported edge, the collapse 


pressure may be shown to be the same as (3.2.1) and is 


given by: 


= 2amesin 
Bs ~ a 1 + sin a | (3.3.1.1) 
Coe Ce 


In the case of clamped edge, we still have [18]: 


At 06 





~)- br4 i i+ sin 6 _ 
) . 2h ‘Sim ob od cos ¢ I 
= ie); to = -] and we find: 
4h sin a 
Psc a 1+ sina (3.3.1.2) 
Log —————— - sin a 
Gos 


Dynamic Response of Simply Supported 
Spherical Shells 


ra 


The yield regime, for the low pressure range, is still 


moor (Pigs. 


5a,b). Then we will obtain the same equations 


ema the same boundary and initial conditions as in 3.2.2. 


Em@emeerore equations (3.2.2.1)—-(3.52.277)) cuit apoly. ihe 


menge Of pressure for which this solution applies 1s given 


mmtoe 2.9 ) 





Ica 


The solution obtained by limited series expansion in 
Were, ol) (3.2.3.6) is still valid in the same range of 


ieessure given by (3.2.3.10). 


Cond) Clonm( 3.2.3. emust be zseplaced by: 





Po PTPs 4,1 Goh 
3h(l + =) a7 5 e (= Ae 39? SS Z (Seo) 
which expresses that at >? = a, Ny S ae 


Since ms seit now be negative, condition (3.2.3.15) is 
now not necessary and is replaced by the one which expresses 
that the minimum of mg be larger than -l. This problem must 
be solved numerically and is not carried out here, since we 
know that in the range of pressure determined by (3.2.3.10), 


MMe olution 1S an Cxact one for the more restrictive evo; 


Moment limited interaction yield surface. 
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Pee Dynamic Plastic Response of Simply Supported Spherical 
coelomicade Of Rigid Pemrectiy Plastic wiatert a ee 


Seeysceene Uncoupled Diamond Yield Surtace 


So@ oily Syeehente evlalsysisie  Hicsteyc lie el 


oie yo Up OOmecc Hage 
As in the case of shallow shell, the yield. regime is 


foment Fags. (Ga,b). From table 2.1.3.3, we have the £ol- 


lowing relations: 


Nees th. “Se=all ca ee | 


ey eo oe ee i re 


Together with Ene equ libri mec lar lomo cea. 2 y= 
(3.1.2.3), in which the inertia forces are neglected, we 
have a system of five equations for five unknown generalized 
seeiesses . 

Brom equation (3.1.2.3) and relation (3-471. 7a cune —- 

(Ss Sin 6) “= (2-4) sate 

Integrating between 0 and 9, we obtain: 

S sin 6 = — (ese econ (354 ol ves) 
; = 


From equation (3.1.2.2) and relation (3.4.1.1), we have: 


n' sin ¢ + 2n, cos ¢ =- cos ¢ +S Sin 9 


® 


Multiplying both sides by sin 9, we obtain, with 


p 


mein ¢ from (3.4.1.3): 


6 hay = Sine cos. ¢ = Wp E(w e a oenm cos sa} 


dp 9 





ZO 
Integrating between 0 and 6 and simplifying by 


pao, we Obtain: 


EOS es 


_ ages 2 
ue ms (p 1) ane >| 


After some reduction, this equation becomes: 


1. pol “Pee ecos.0g 


a ee ae ee A 
"6 2 2. Jee eos % (3.4.1.4) 
Then from relation (3.4.1.1): 
ee Da bela Cos. oO 
Ny 5 + ae ae (3.24 lee 


Brome(3..6.2.1) and (824eiReZie we have: 


h({m! sin ¢ + 2m, cos 6] = h cos 6 + s sin 


> b 
Multiplying both sides by sin » and with s sin $ from 
24.1.3), we have: 


e) (m sin? Cle = hh Sine Oeicecr on. (p-1) (sin y= sin b cos) 


° , 


Integrating between 0 and » and simplifying by sin’9, 
we obtain, after some reductions: 


Cos © 


eo i =e 3-2 ee 
Then from (3.4.1.2), we have: 
ee Q.4.1.7 
At ¢ = a, we oe am = 0. Thus, we derive: 
p=p,=l+h pees (3.4.1. 8a) 
p. = 1 + h cot? ' (3.4.1.8b) 
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If the base radius B is taken as the reference length 
and if o* is neglected with Fes pece cow lt Cancers hewn 
fgat (3.4.1.8a) reduces to the expression of the shallow 


Sim@ell case. (2.4.1.8) 


eri rWeme xpress] Ohms se4 slyctd) kOlne, equates 


ero. t.4)-(3.4.1.7) become: 


ny = - = -h F (4) Ear) Ale) 

ne = - 5 th F(9) (3.4.1.10) 

my = - F(¢) (3.4.1.1) 

ae! . 
mo = 5 + F(¢) | (3.4.1.12) 
where 
_ JS cos ce E=COS o 

SWOUSSED TES a eee d Sie eat 
1 tan? s 

F($) = (3.4.1.13b) 
tan~ > 


Beem (3.4 .1.9)—-(324.1.212) and fromvtne. tact. chow neens 
small compared to 1, we can see that the necessary and 
mreeeelent condition for static admissibility is, for yield 


regime 3-5: 


hol 
hoje 


< F(o) < iene Os > <8c. 


From (3.4.1.13b), we can see that F(¢) increases from 
O to 1/2 as ¢ increases from 0 to a. Therefore the solution 


1s statically admissible. 





ZS. 


The flow rules associated with regime 3-5 are, from 


table 2.1.3.3: 
e.-e =0 (3.4.1.14) 
and bs k, = 0 (71 15) 


ean) (ome) Vand) (3.1. 123m, equcdemeonm(c.4 a). ka) 


yields: 


V 


II 


A sin 6 (3.421286) 


At ¢ = a, we have v = 0. Thus: 


A Sin a =0), Or cancers tne sO 
A = 0 
and fs 6 (3.4,0i7) 


Bom) (3.1.1.2), (3.1.164) and with 3346) ecu aeron 


eeo.).15) yields: 


w' = B Sin 6 


Integrating, we obtain: 


w=-Bcos¢+C (Sy Aleeinsh 


Using the boundary condition at $¢ = a: 


w (a) = Q, 
we obtain: 
w = B(cos a - cos 4) 
Let We = -B(1l - cos qa), we have then: 
ee ae COsme = COs 
eR gi 5 tooo) 
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The condition for kinematic admissibility are: 


a. = @ S 0 and k = k > 0 


w) 0 d 6 
From (3.1.1.1) and with v and w from (3.4.1.17) and 


3-4.1.19), we have: 


_ ° 2 COs = COs 
e = @ ee, = ee ee ee 


d s) or | = acos Va 
(OS Sono) 


Which is negative since Wa O and @ < a. 


femelle 2) and wath v and w from (3.4... 12 mand 


(3.4.1.19), we have: 


eee . COs” 
ee Kg aS i> cos. 


Sor Sr) 
With w. > 0, Kk, will be positive if a < 1/2. 
Therefore the solution obtained is an exact one for 


Seoea/2 and the static collapse pressure 1s given by 


Wome. t.Sa) or (3.4.1.8b). 


At ¢ = 0, w' = 0 and there is no Singularity at the 


apex of the shell. 


B.Clamped Edge 


iecOlLucIOn 
As in the case of shallow shells, we would expect a 
hinge at the clamped edge, and the yield regimes are 


(Bigs. 6a,b): 0< $< 8: regime 3-5 


BS 0-5 (Cee Pegi me eem. 





DOS. 
wee ) = Ps regime 3-5 
In this interval, the 5 equations to determine the 
five generalized stresses are still made of (3.4.1.1), 
Merete. 2) and (3.1.2.1)-(3.1.2.3) where the inertia forces 
are neglected. But we cannot use the boundary condition at 
fe © in this case. Therefore equations (3.4.1.3)-(3.4.1.7) 


eel remain valid. 


The flow rules remain the same, and without using the 
boundary conditions at ¢ = a, the expressions for the 


memocaties are €3.4.1.16) and (3.4.1.18). 


moe? < CO: regime 3-6 


From table 2.1.3.3, we have the following relations for 


Begime 3-8: 


Ny t+ny = Sl (3.4.1.20) 


NA = Mae (3.451220) 


With (3.4.1.20) which is the same as (3.4.1.1), the 
equation to determine s sin $9 is the same as in the simply 
Supported case and since s sin $? must be continuous across 
mec poundary ¢ = 8, the expression for s sin @ 1S the same 


een (3.4.1.3) 
Similarly, we obtain the same expressions for me and 


Np moman (3.4.1.4) and (3.42.1:5). 





Using the boundary 


wy 
we must have: 


Eaves 


ZG. 





To determine it and Ma, we Nave we Seen) anc 
fer4.1.21), which give: 
mi = CoS 9 , LS Sin $ 
¢ sind h sin 6 
With s sin ¢ from (3.4.1.3), we have: 
mt = COs o _ prt | 1 _ Ses cos 9) 
ob Sin h ‘sin 6 Sin 
Baeegqrating from 6 EO ¢, with ee = 0, Since Ms 
must be continuous at >? = 8, we obtain: 
= Dae sabia 2) db Leesa oe ces 
my (1 + h ) log Sin 8 2h (log l-cos 8 log cos 4) 
(324.1 5223) 
Then: 
7 io SBCh wl oss dl I+cos B 7 l+cos 9 
“aa an es h P=) lo “Ssin B 2h (logy cos B LOSS -~cos ) 
(3.4. 1223) 
ines low GULeS for Lregliime 3-3 aicey. ot s@lmee ccleaner roe: 
Cg = 2 (34a ia) 
Ay et 2 
Ko ns (3.45.1 325) 
mwom (3.lel.t) and (35 oie pane wii e 4. eee 
- have: 
v = D sin d 


condition at ¢ = @, we Nave: 


Dsinaz= 0 , and since sin a # 0, 
D= 0, 
v = 0 (8 < 6 < @) 





Zee. 


But v must be continuous at ¢ = 8, and from 
me4.1.16), we have: 
A sin g = 0 
Since sin g # 0, we have: 
A= 0 
Thus 
ve0 (0 <9 <a) (3.4.1.26) 
rom see lee) Se a Ga acl 3. 4 one 
(4.1.26), we have, 
i? «(are 
Sin 9 
Integrating from » to a, with w(a) = 0, we have: 
= a lo i+ cos a _ LE Lap Ses) ] (20-45 low) 
2 3T = cos o 3T — cos od ne 
c) Boundary Matching 
Ae Oo, we have, from (3.4.1.27): 
-, E 
Ww)? aan a 7 0 


The slope 


Possible if at 


With this 


1 


w' will be discontinuous at ¢ 


O. =o, Weave. 


Sondiltr1on tand ef rome. oo) 


Sin a 

9 hf{l + rect Sin 8 
al l+cos 8 

aye) Sameer l+cos a 

2 Ww COowE Iker: ST A 


= aes 1S 


, we derive: 


; Sarde 0) 


Sin og 


J - log Sin 8 


3 





Zee 
At > = B, a MUS Ee COME mNOUS . (vo inCewa Adie ace, 


the yield regime changes from 5 to 8, this is possible if we 


have: m,(g) =m (3°) = 0.9 Alse vw muse be Continuousmanc 


p p 


since m (8) | # 1, w' must be continuous as well. 


The continuity of s and n, is not in question since 


p 


ieey are expressed as continuous function of ¢ in the 
matcerval (0,a). The continuity of v has been used to find 


mac v = 0 throughout the shell. 
Eom can = OO, We Wave , iwi el (Siro lees ee ar i oe ee 
B= B sais (34. Lee) 


From wl), = 0. we haveya wii ss4 oi 1S) e242 one 


more. ..29): 


7 & ee i oF COSt iG I ee Sane 
C = B cos g + = B sin’ 8 [log aoc NOC Rane 31 


(3.41430) 


Ween thas result, (374.1.93) becomes: 


w = B cos 8 - B cos 9 


i I a Sfa2 @ los en -COSaie wa log IP SECos 8 


Z i= ces @ i COS as 
At »¢ = 0, we have: 
r( d= =w = - B 
w (o=0) w, = B cos 8 
EER eine 2 (loqe eS See a ee Osa 
2 b J [= cos CO 37 I = cos B 
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Ors: 
; (33245) 
W | 
B= - , a 
i eer nee Ia cles = Lee See ae 
I COs 6.4 5 Se B (aera ane 2 log ge ae 
With this expression of B, (3.4.1.18) becomes: 
jee GOs! 1 a Bares a 
_ (COSmeCOce TD as > sin B Llogi— os 6 i eer 
a Lee l+cos 8 l+cos a 
Jl" aleve; [Sais > ey NG eo Lex; I-cos 8 = log 12eCcueen 
OES ie 5) (374 ee) 


With E from (3.4.1.29) and B from (3.4.1.31), equation 


e4.1.2/) becomes: 


oe et COS On i-aatCcos ae 

ae - . aaa B {log i= COs. @ OS ee Ser: i 
oO. | ee: Ie eosiG | i oS eS 
1 cos 8 + 5 sin® 8 [log Some log =r coaene 
(oe Oc) (3 les Sy) 


ihe can be verified that withwts.4. 1.32 )wand seo 1 36s 


and w' are continuous at $6 = &@. 


Se 


, -, _ +. | 
From my]. =O), WR ae VE myth ) my (8 ) C4 iene. hehieol 


@oee.1.6), we derive: 


1 + cos 8 
l1 - cos &8 


Dea = Le ia (3.4.1.33a) 


1 +h cot? 6/2 (3.4.1.33b) 


he p 


By comparing (3.4.1.33a) edbove @ith (3-4-1. 28), we nave 


ero eionm to determine ©» when @ Tem enoewmeoe wate can be 





Ze: 


written as, after some reductions: 


i 1 + cos 8 _ i, COs ae! 
G(P) = 3 G+ cos 8) (log T——~cos-p 7 199 T= cos a! 


mie leo) = 2 Aes Se 0 





Sins 
Cor eee 
It can be shown that: 
Cl) Sahil sne Cae 
ea 6 
G(0) = 2 log ee 0, since cos? <1 
cos? = 
2 
emda that: 
Glo) = -(1 - cos a) < 0 


Moreover the function G(g) is continuous in the inter- 


mei (0,a) and that its derivative is: 


Al ; 1 + cos 6 1+ cos a 
' eee AI tot aes EI ec Sc es 
G'(8) = Feat 8 (log 1 - cos 8 log 1 - cos i. 
l= ) COS > 
sine Bee 


which is negative. 


Thus G(8) is continuous, decreasing monotonically in 


mae interval (0,0) and G(0) * Gita) < Q. 
Therefore there exists a unique value 8 such that: 


Ens) = 


For instance, with a = 90°, 1t can be shown that 6 = 60°. 





Zale 


With P. PEM sets eai, eCulations (3.4, .).4)— 


were. /) become: 








ae = - h F(4) (0 <<a) (3.4.1.35) 
ny = 7 5th F(6) (Oem) (3 iaeieses 
my = = - F(9) (0 < 6 < 8) (3.4.1.37) 
me = 5 + F(Q) (0 < $< 8) (3.4.1.38) 
where 
igh = 3 PSOE Fe Oke (3.4.1.39) 
and (3.4.1.22) becomes: 





Oe oghes One co omy ae 


o Sin 8 1 = cos 2 J PSO a) 


- log Lite cos By + log sin 6 


1 ees) 8 sin 8B 





(8 < @ < a) (3.4.1.40) 


Then: 


7 ind , i+ See a 
- 1 + log sh : “ J cos ; Glee IL cos 


an I= eos ae 1 - cos ¢ 


1 + cos By 4. ee sin ¢ 


ieee IL = Sisios” Sin Bg 


CG °< ao: ay Seas) 


The velocities are Givens (324 e267, (3.461 .32) 


eumee 3. 4512433) . 





Pee 
fepkunematic Admissibility 
a) 0 < $9 < B: regime 3-5 (Figs. 6a,b) 


The conditions for kinematic admissibility in this 


interval are: e.=e <0, k, = me 


Q 0) 
Mieeyl {see elle A eevevel ee-elat Vv and W Bae Sy LOR at eligie! 


moeo.1.32), we have: 


it ACOs ee Ue le cos O. 
Ik Melee as esi - cos a 


Lg tCOs. 6... lo ces Ot 
Lae COS: 26 JT = cos a 


(cos ¢ - cos 8) + 5 sin’ 8 [ives 


. I => Sale tows) 4 = Siva 8) (bie 


-— W 


La Cos 1 eres ale: 
i+ cos 6 - log —————— > 0, 


Since 8 < a, we have log in Eoome > J - cos o 


e@emrunction log eee 


whose derivative is - 2/sin x 
Beeeg monotonically decreasing for 0 < x < 7. 


Therefore Cg 1S negative since We Zt) 


Magen (3, L.2.2), andewren (354.1526 mane. (oo oc 


we have: 
a hw, cos @ 
: rae it + cos B _ 1 + cos a 
1 cos 68 + + Sin 8B [log t= see log ins ea a 
(O < > < B) 
k. is positive for 4 CS ii ac 


0 
bee < 6 < a: regime 370 (bigs cae) 


The conditions for kinematic admissibility in this 





ZS. 


interval are: eg = C4 So. Ky = -k, >? 0 


p 


Wise ee 4 el 26 lemeiices 14 les Soe We wilanto : 


Wn 1 + cos ¢ _ i + cos 4 
-_ yy OU a [Akeley ae $ 10g += saaea. 
=) ane iL aoe. IM eee! js MEP cos & 
Ve COS moe: 7 Sin B [log I - cos 6 159%: AR eel o)s) 4c! 
(8 < @ < Q) 


fech 1S negative for B < @ < a. 


» Seah sine se 


hw. sin* 8 cos 6 


a Tl cose at sin? & flog PESOS. yo, LE OSB, 
1 SO 9S Ms 5 Sin B [log 1 - cos 8B ce KEG a} 
(B <> < a) 


Mmemeesttive for B < ¢ < a@ < T/2. 


Thus the solution is kinematically admissible for 


Sam 1/2» 


Bemeotatic Admissibility 


The conditions for static admissibility are: 


= Sere (0 < 6 < a) (37.4 sdi42)) 
Ska . < 0, CO aoe ae) (eo Are 43) 
Oo<m <i, (0 < » < a) (3c Alek, Aa 
O<mcl, (0s 6 <8) (3.4.1.45) 
ale i, < OU, (B < > < a) (3.4.1. 46) 





Zake, 


From (3.4.1.39), we have: 


F(o) = 5 cot? E tan? 


No}-e 


From this expression, we see that F(¢) increases with 


¢ for 0 < ¢ <a < 7/2, and 


_t 
P(B) = 5 
and the maximum of F(¢) is at 6 =a: 
eo 2 B 2 Oo ul 
F(a) = 5 EOE 5 tan 5 (a < 5) 


Therefore, condition (3.4.1.45) 1s satisfied. For con- 


ppmmetens (3.4.1.42), (3.4.1.43) to be satisfied, we must 


have: 
hae ie) <5 or 
2 a 2 B 
h tan 5 SEU. 5 


Since tan a/2 and tan 8/2 are of the same order, and h 
is very small compared to 1 for thin shell, this condition 
is expected to be satisfied, and therefore No and Hy are 


admissible. 


Mor CONGLeELONS 91S. 4el 46) \iCmicierusiateOm nt Ge ord Oe: 


Ieee Soa 2 Ces 0 


> Sane 


for BS > < @ 





Za 


Laas, ™ 4 decreases monotonically from 0 to-l as $ 


increases from 8 to a and is therefore admissible. 
toe admissiba lity of M 6 results conseguently. 


Therefore the solution is an exact one for the yield 


Surface selected and the static collapse pressure of the 


emamoed Shell is given by (3.4.1.33a,b) or (3.4.1.28). 
ites coluct Ones Valid tor ugg / 2. 


Ht the bese radius cf the shell as taken as the reis 
erence length in the non-dimensionalizing procedure, it can 


Z 


be shown that with a~ << 1, expression (3.4.1.33a) reduces 


memeete shallow shell case. (2.4.1.34) 


feaee Dynamic Response 


A. If the peak value of the dynamic pressure is 
slightly above the static collapse pressure, the yield 
memrme will be still 3-5 of Figs. (6a,b). The flow rules 
@eee.i1.14) and (Sys eae Jk) remain valid, the boundary condi- 
tions remain the same, and therefore the velocity expres- 


sions (3.4.1.17) and (3.4.1.19) remain valid. We have then: 


v = 0 CO ct) (ee ee a 
ae cos 9 - cos a (eota4 5 een 
O l - cos a 


COPS Oe Ci) 








ZS 


Bemeese Gi e po meene relations (5.4.1.1) and (3.47182) 
remain valid. With (3.4.1.1), (3.4.1.19) and (3.1.2.3), 


we obtain: 


.. sin > cos 9? - Cos a Sin 
S Sin ' = yw — 
( b)' = yw 





== 1) Sibi 
a COS mcr 


Integrating from 0 to $9, we obtain: 


| Yo IO dp 
S Sin > = 5 oe Sin- 6 = cos a(i — cos 9) ] 
- (p-1) (1 - cos 6) (oe ee) 


ee leeee Wal te OIC 5:5 4 eaten ) re GCI CRe Cli cie Onlin ( Senden) 2 alone 


have: 
Oe Sin o> + BE COS EO) =seces (0) et es Sand 


Multiplying both sides by sin 9, and with s sin ¢ from 


morea.2.1), we obtain: 


h ax (My Si? Oo) = See Sting th 


Yw 
O dk 
* T= cos a 7 23 


> -¢ - cos a (1 - cos ) sin $6] 


=e Gey Ib MGR =e clots ip) Sia 


Integrating from 0 to ¢ and simplifying by sin’ ¢, we 
obtain, after some reductions: 
Y¥o 1 


h(l-cosa) ~ Ttcosé > 7 P< Sica ee ——(1-cos$) ] 


A 


oe o- 
¢ 2 
pot. k=cos® 


a eee he 


: eer] SEAL ae) 





Zig 


At > =a, M4 = 0, from which we derive: 
me Df 7 1s Peles 6! 
ig = ae J = aaa. 


With P. given by (3.4.1.8a), we have: 


ee oe 
is 5 (p p.). oe), 
With this result, equation (3.4.2.1) becomes, after 
peme CLearrangements: 


a Ce oa 


s sin ¢ = 5 P—Ps— (1 - cos 6) (44, 


Pei COS: 3 - cos @] 


Se CO) Se el eo Suaes) 
Using again P. from (3.4.1.8a),we have: 


_ _ ose eres 
= Ps a dl GS CIS eae 


Brom this resule, werger: 


P- Py 1 - cos 0 








S$ sin ¢gd=- (1 + 2 cos a@ - 3 cos 6) 
4 I= Vees a 
i + cos @ 
h ‘ea ese f (1 cos d) Cee 7) 2) 


Samitlariy, equation (3.4.2a2) becomes, atter some 


meaGuctions: 


dL it ees i iS COs ue 
I = a a 


d = COSme cos a 
i oes Pe COs ay EOS Wi SOS ae! 
4h IDS elers™ (3) ee OS aN 


oh e2 ao) 





PINS) 


With v = 0 and using relation (3.4.1.1), we have, from 


eemacion (3.1.2.2): 


ne Sin o + 2n 4 cos $@ = - cos @ + S sin 
Multiplying both sides by sin 6, and with s sin 9» from 


(ee. 2.4), we obtain: 


Qs 


=-(n sin* ¢) = - sin $ cos 9¢ 


oo 





32 p = 
a (Si Oh feo COSmcqskn i) 5 os ieQmeCes ma) 


4 lL - cos 


i COS: ee . 
h arr! db Sin @ cos 9) 


Integrating from 0 to ¢, and simplifying by sin* 4, we 


obtain after some reductions: 





ee a eee | 2uces es ees 
w) D 2 Siees 1 + cos 
_ Cor ediaeeos) 
¢ PYPs 1 - cos @ cos 6 = CoS «a 
4 I= sete te! 1 + cos 9 


Comparing 3.4.2.6) above with (3.4.2.5), we see that 
we Can write these expressions in the following terms: 


my = 5 = seas oy (3.4.2.7) 


Ny = = = -hF (cos od) (32 tae) 


with F(cos ¢) now defined as: 


F(cos 9) = lil+cosa,1- cos % eae sae 
21- cos ® 1+ cos ¢ 
- PPPs 1~ cos ¢ cos ¢ - cos ¢ 





4a 1 =]=Yeos & 1 + cos ? 





Then with relations (3.4.1.1) and (3.4.1 


ae 
My = aaa F(cos $6) 


NO] 


Ny = + hF (cos ¢) 


BOM S.4  2ao we Oboe amin: 


aa 
Wo 


i 


Nee 


fg. (lser ) - pt] 


Cae, 


2), we obtain 


(3524-22 e0) 


( 37Gee2 cll) 


(3.4. Zep 


With this result, the velocities are given by 


Mere .1.17) and (3.4.1.19). 


Admissibility 


ieee Kinematic 


ReOm themsOlucion Of tChetstatnme orootem, secclon so .c as 


we have found that the solution is kinematically admissible 


for 


oS WZ and 


w > QO 
Oo = 


With Wo from, (3.422.12) awe seextaas Wo will be posi- 


mee, fOr 0 < T S Te, where T 5 is determined by: 


i Ne = 
Pp, (1 e ) Pt; 0 


At T = Ter 


2. Static 


Ww = 0 and the moclonmceaces. 


(32447 so) 


FOr regime 3-5, Figs. ((6a,0)) ec sconate1ons fOr Staric 


admissibility are: 
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“ee oe ye eye SO 
0 ny, sgl ; 0 < Mp srt 
From the expressions of oe Nye Mgr Ng Ln (34267 - 


(3.4.2.11) and since h is small compared to 1, we see that 
mmemmecessary and Sufficient condition for these inequali- 


memes tO be satisfied is: 


J 
NR 


< F(cos $) <5 (3.4.2.14) 


vameme F(cos ) has been defined in (3.4.2.9). 


It can be shown that: 


19) 


O 


lake P-P. > 0, or: OS tees dHexe) >. LY 
S 


Somat ctloOne(o.4.2.14) Wali bes sacieonica elt: 


[l + Yl + cos a]? 


aoe os (3.4.2.15) 


36 Dee 01 OS ih eee: 


Gemaition (3.4.2.14) will be satisfied if: 


Z yb 
P.-P,° < 4h 


(3.4.2.16) 


Therefore the solution is kinematically admissible if 
the peak value Po of the dynamic pressure satisfies 
(3.4.2.15) and (3.4.2.16). The lesser of the maximum 


values of Py Geta rreed serOm 5.4.7 >) none Miao 2G |) Wael 





Zee. 
give the upper limit Por of the peak value p.. If the 
shell is not very shallow it has been found that condition 
(3.4.2.16) is much more restrictive than (3.4.2.15). Then 
the value Poy can be obtained from the system made of 
mre. 2.13) and the equality (3.4.2.16): 


l-e ‘f) -~ pt. =0 (3.4.2.13) 


a Syer 


(S) oe 7 An 


a eeoae 1c (3.54, Ze) 


Ps 7 POL 


te 


is given from (3.4.1.8a,»b). 


@mreconc!]usion 


Oe Ss SS Sot 


wnere Por, is determined by the system (S) above, the 


solution is statically and kinematically admissible. 


The stresses and velocities are given by (3.4.2.7)- 


wre. 2.l12) and (3.4.1.17) (3.4.1.19). 


The final displacement distribution is obtained by 


integrating the velocities from 0 to Te? 


Wien = 0 


COS 0 waecoOS me 
1 - cos a 


YWe = YW, (uae) 


where ioe is the final central displacement which is: 


= = _i 2 A 
a= ey pote > Pot ] (e422. em 





Ds 


where Ti. has been defined in (3.4.2.13) and which is the 


Mmastanc woen the motion ceases. 


iMremene nay —abserbed is - 





Te * au 
= yee 
Be = | jrat xX 21R aa bdo dt 
6) O 


Or, in terms of non-dimensional quantities: 


_ No 
Be Seep 
d oe _ 
a (Ey R 7 qdt = Todt: 
Oo 
Tp oO 
a 2 ; : 
E obs = 27R No | | pw sin o> dd dt 
O O 


with w from Coco, POE. we Slotanunn: 


2 





3 CL COcmey Py Re 
E ea 27R° NS = | 5 (l-e fy 2 
a Ay 
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(b) Velocity and Displacement Distribution 


Stress ana Veloctev, Displacement ea rst urtelutons 





Deep Shell - Uncoupled Diamond Yield Surface 
Pe gee 
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3.5 Comparison of the Results 


3.5.1 Two-Moment Limited Interaction Yield Surface 


From the results presented in (3.2.2.1)-(3.2.2.7) we 
see that the solution obtained is rather complicated, even 
in the low pressure range where the whole shell yield 


in one regime. 


The limited series expansion solution (3.2.3.2)- 
ee. 3.6), neglecting terms of 4th order and higher, could 
give a better indication on the upper limit of the peak 
value of the dynamic pressure. Although the problem of 
static admissibility has to resort to numerical méhod, 
this solution gives an indication on the maximum depth of 
the shell when the dynamic load is a rectangular pulse, in 


fe2.5.16). 


An analysis for a higher range of pressure was 
attempted,by using a series solution, the results obtained 
become too complicated, the algebra too involved and no 
meaningful conclusions could be reached. We did find, how- 
ever, that at higher pressure range, the shell would yield 
' in at least 3 regimes like the medium pressure case of 
shallow shells. The solution which assumed the shell 
yielding ee 2 regimes arrived at the same conclusion as 
in the case of shallow shells: when $6/a reaches a value of 


poome 0.5, e, becomes positive and the solution is no longer 


? 


kinematically admissible. 
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Tie wupper Limit of Py in the low pressure case 
[me .5.10) 1S comparable to that of the shallow shell, 
Wee .2.1.21), together with the fact that the shell does not 
Peetd under 2 regimes at higher pressure and other similar- 
ities found in the evChin aro LLCS considerations suggest that 
iteGe 1S a strong analogy in the solutions of neo and 
Shallow shells, and that any further development of the 
deep spherical shells could be made by adopting tentative 
solutions which are constituted of the same yield regimes 
mom the shallow shell solutions. Obviously the choice of 
yield regimes does not solve the whole problem but it gives 


some help to start the solution. 


Bosc UNCOuUo led Dilanond. y Hele stnnece 


mnie Soluticen obtained | 3245 2 ee ( Se fee eee 
(3.4.1.19) are simple enough so that the whole problem: 
obtainment of a solution and examination of the admissi- 
bility eft this SOLUEION, Ceaenwbersouved analyze ucarly) Lom 
Shallow shells as well as for deep shells, the generalized 
stresses are expressible in terms of a single function 
whose variations could be obtained without much difficulty. 
The membrane stresses are small compared to the moment 
resultants, and it appears as though they were "trapped" 
within face 3 (Figs. 6a,b). The meridional displacement 


is zero, the normal displacement has a cosine distribution, 





jd oa Ne 


the maximum being at the apex which is also a regular 
point for the displacement function (no discontinuity in 


EEE@pe) . 





eee 


4. CONCLUSIONS 


By using the two-moment limited interaction yield sur- 
face proposed by Hodge [18] and the shallow shell approxima- 
tions, we have obtained the solution for the dynamic SAenensa 
of shallow spherical shells for medium high pressures which 
Satisfy the ineguality 6h SP aaee Sen wie Bede Del come Mme 
MatiOo r = depth/thickness of the shell (e.g., for r=1, 
me 19.1, for xr = 5, A = 13.2). We have seen that there are 
three stages of motion: in the first stage of motion 
HON < 1 Sta, the shell yields under 4 regimes, in the 
second stage Coe SS ime under 3 regimes shovel, abgiieli 7 els akrel 
stage Oey es T ¢) under 1 regime. The solution obtained 
is rather complicated. The restriction on the upper limit of 
the peak value of the Anam pressure for shells whose ratio 
Gepth/thickness is larger than 3/8 arises because the moment 
resultant iy must remain positive throughout the duration of 


the motion. This requirement is most critical at the instant 


When the shell comes to rest. 


A first attempt to increase the upper limit on pressure 


would be to allow m, to be negative at some instant during 


p 
the third stage of motion (1, < 1 < t,)- Then from this 
instant the shell would collapse under 2 regimes: probably 
5-7 and 5-86 and the subsequent behavior probably woud not 


return to one regime again, since regime 5-8 would most 


likely be spreading towards the center of the shell with 





Db DES 
time. If this problem can be solved successfully, the in- 
crease of this upper limit would not be ety to be 
important since sooner or later, and sooner rather than 


later, n, would become smaller than -1l near the outer edge. 


s 
Then at least another regime will have to be added to the 
collapse pattern of the shell, so that from this instant, 
the shell would collapse under 3 regimes, separated by 2 
boundary circles, probably shrinking towards the center, and 
catching up Or widening the gap between each other. Still 
this upper limit of the peak value of the dynamic pressure 
would be limited by the fact that this pattern may occur 


emmy during the third stage of motion (T) °S ce Te) when the 


shell starts to collapse in one regime. 


Thus the complete solution for shallow shells of prac- 
tical geometry (ratio depth/thickness larger than 3/8) sub- 
jected to any magnitude of pressure would be the one that 
starts with at least 6 regimes separated by 5 boundary 
circles all moving towards the center. The mathematics would 


be complicated. 


For deeper shells, the problem is still more involved 
and only a solution for low pressure range has been obtained. 
The limited series expansion of this solution has allowed us 
to gain a better understanding of the upper limit of the 
peak value of the dynamic pressure and to draw some analo- 


gies between solutions for deep and shallow shells, and thus 
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to decide on the direction to start a solution for higher 


pressure range if one has the desire. 


One noteworthy aspect of the solution for the shallow 
shell according to the two-moment limited interaction yield 
Surface is that the meridional component of the velocity 
1S of the order of 0(Z/R) of its vena component so that it 
cannot be neglected, since the approximation is that only 
(278)- is negligible with respect to 1. Another aspect 
Concerns the variations of Nos which are quite small so 


0 


of trouble if one desires to go to higher pressures. 


that n, remains very close to -l, and which may be a source 


With the uncoupled diamond yield surface proposed by 
Jones, we have seen that for spherical shells, the mer- 
idional component of the velocity is zero, and No and am 
have been brought within a narrow interval around the mid- 
point of the corresponding yield regime. The solutions 
obtained, for shallow and deep shells, are quite simple and 
the problem of admissibility, which is usually as difficult 
as obtaining the solution itself in the case of the two- 
moment limited interaction yield surface, can be examined 


without much difficulty. It is unfortunate that attempts to 


obtain solutions for higher pressures did not succeed. 


We should note, however, that the uncoupled diamond 
yield surface does not monopolize in denying us a solution. 


Im the limit analysis of spherical caps with the two-moment 





Zour 
iamLtced interaction yield surface [18], the cap angle 
has been limited to a value smaller than 90°; otherwise 
the value of N 4 at a would become positive. We have tried 


to extend this solution to higher value of the cap angle 


by assuming the following regimes: 


AN 
O- 
A 


Hom 0s Se: pecmbero— 7 (macs. 4 ayo) 


IMO) BIS 


A 
© 
A 


SC icieve abuts (927 (Ualepse leiylo;), 
At @ = 8, we obtain: 


and 


re (IB) 


| 
= 


| 
, 


ng (B) 


The first equality comes as a definition of the angle 8 


and the second as a consequence of the continuity of Dae 
9 should have a minimum até. 


— wit ie 5 ae Z 
We have found that at @9 = 8 , Ng = O and Ny = Sane 07, 


resulting rather in a maximum at @¢ = g, Helse pele pyeis. (exons — 


Mrs in the interval (8,0) n 


Se . 


Other combinations made of regime 5-7, 6-7 and 6-8 


were also not successful. 


The common factor between the Hodge's and Jones' yield 
Surfaces 1s the absence of iuBeyaeesee between moments and 
membrane stresses. Perhaps for general shells made of 
rigid-perfectly plastic material, the interaction between 
moments and membrane stresses would not be negligible and 


more sophisticated yield surfaces should be used. 
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In the case of static plasticity, She general theorems 
allow us to construct bounds on the collapse load obtained 
from an approximate yield surface. In dynamic plasticity, 
no such general theorems exist, although Hodge and Paul [13] 


gave some useful hints in that direction. 


From this study, we could see that the problem of 
dynamic plasticity of spherical shells loaded symmetrically, 
even with the use of approximate yield surfaces, is rather 
complicated, and a complete analytical solution for a wide 
meamec Of pressures is not likely to be possible in the near 


future, if at all. 


It might be mentioned that the only previous studies on 
mac rigid-plastic behavior of Shallow shells are those of 
Hodge and Lakshmikantham [23] who obtained lower and upper 
bounds to the static collapse pressure of shallow parabolic 


and conical caps with central holes. 


It is hopeful that the analytical solutions presented 
herein for the dynamic and static behavior of shells would 
be useful in checking the validity and accuracy of numerical 
programmes and to assist in the development of approximate 
methods such as those developed in [24]. It appears that a 
numerical procedure would have to be used in order to extend 
the existing solution to higher pressures. Then, relatively 
more complicated yield surfaces can be used, such as the one 


proposed by W. Flugge and Tsuneyoshi Nakamura [22]. 





Zea 


Thus, it is believed that numerical methods and more 
Seonisticated yield surfaces would be the direction for 
future researches on dynamic behavior of shells of revolu- 


eon . 
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Maximum central deflection according 
to the two-moment limited interaction 
yield surface. 


Maximum central deflectm according 
to the uncoupled diamond yield sur- 
face when the static collapse pres- 
sure is made to be equal to that of 
the two-moment limited interaction 

Surface. 


Maximum central deflection according 
to the uncoupled diamond yield sur- 
face inscribing the two-moment lim- 
ited interaction yield surface. 


Wonc = Maximum central deflection 


according to the uncoupled 
diamond yield surface circum- 
scribing the two-moment lim- 
ited interaction yield surface. 


(Z/R)/h 
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p 
110 h 
100 
90 
80 
70 
60 
50 
40 
30 ope = static collapse pressure according to the 
Ls two moment limited interaction yield surface (ss) 
and the uncoupled square yield surface (SS) 
Pon = static collapse pressure according to the 
20 uncoupled diamond yield surface (DS) 
Pory’ Pois? PoLD = upper limits of peak values 
Ot CYNaAMIG prescuUrcs (or Sol- 
utions using (HS), (SS), (DS) 
{ °F respectively. 
e 4 8 2 16 20 (Z/RVh 


G.9 = Static Collapse Pressure for Simply Suso0rted Shells and 
Upper Limits of the Peak Values of Dynamic Pressure 


(Shallow Shel!) 
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APPENDIX Al 


MEDIUM PRESSURE 
ea Sone. ope (She. 


betermination of the Unknowns : 
mero 7e «© s Aaa SMe no een Ge 
A, coy / Gar Cor R2 Cy! Ky pi Us Re E crane o G, a SS 


J ——— 


To determine these 12 unknowns, there are 12 equations 


: tt pe 
Expressing the conditions of continuity of yw, Yw,Yv, De 


Ss, as at each of -the boundaries X=U, y X=Us- 


Let ele be the jumo of the function f at x=u, we have: 


ROM. BU2 ste 2 oie 2 ee Oe ANON U2 eee oe 


yw'] = 0 > C, =7— cs 
i O Ire 


Bron 242:5.2 6 2623) Pane 2 2 ee a 


Aj 
pz “1 


-|- 


Wy, =0-> C= Q 


PMO Makt2 2. 2. 2 « 2 2 cm? oe eee a 


e 27,8 

pes ee) 
: Di Dyn 1: 7 R 22. 73) ae 
oa, = Oo? Mey eypecn 9 a Soh) ey RU (C7 R Buy) 


tL 


Prom (2.2.2138) andw24252.2-- > =. 


. 2, 3 _ 2 


' = So 
yw' J Q > R 


From (2.222.2238), Ancien 


Pu, * QO >+ C = =K 


From (2.2.2 02237) "onede2  2ece a). 





Za, 


° Zo 2s 
eee 
: R 22 : Zs 22, ° 
nl Oe ceria, eee ae ee je Set — 
Pu, > Sa7R as u5(C = Bu. a K (u., 1) 
There are 6 equations for 8 unknowns: A, eos ec 


R 
2 
4 We ° ° 
pe Oy K, Uy, U.- Therefore we can express 6 of them in 


O 


terms of the other ‘2 which are C and Uys and we obtain: 


K=-C (Al. 1) 
Vs e = e 
= 8 — ae (Adie ey) 
Co = -C Giese) 
47? . - a 
ar C, = C- 0 (Al. 4) 
ie? “ 277% 4 
A = -C(l => Sl ) VA oS) 
R 2 
BC=10 sme, 
2 =a ewer ~_ — 
us ei uy Z Uy CON IL ia) 


So far we have used the 6 conditions of continuity of 


e e e e 27,8 e 472 e 
wr, YW and Yv at Uy and Us to express A, a By Cor Rey! 


K, in terms of C, Uy and uu. and obtain a relation between 


Uy, Us and C. Now we will go on to determine the 4 remain- 
2 Dee Le ; 

ing unknowns =5-D, — ; =o F, (GT alig\  eehienns, ene (C ysheva Uy, Us 

Epaecing the 4 conditlons Of geen time wo recs ond NxM at 


Uy and U5. 


Differentiating (Al.5) with respect to tT, we have first: 


oe = ee 27, 2 ee ”) e e 
A = C Rn? (Cus a5 2c ,uou,) thie 4) 
From XS], = Oye Using ( 2aeece 2. eee 2 22592, 2 A) 
1 : are aN 
with x = Uy, having A, <3, Cyr =e, by differentiating with 


mesoect to t/(Al.5) 5) (AZ) (oi ee ee ees peel vely 
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and substituting the results obtained into (2.2.2.2.21) 


BiG me. 2.242 4 ONe, Wenn omen 

















y ee 2 
42"; E-Q ae et ee 2+ (Gu2 +2Cn0 U5) 7 2. 
R R 3 ti! 42, 2 
RO—_ - RR 5 - a) 
Bg RO 47° /R* 
, (Ades) 
I. 7 ee oe 3 a 
+ %m (35 D - Ga) a = 0 
From xs], = 0) USING 6222222. 45) eon 1002 oe sO) 
Z 
with x = Vos and making the necessary substitution, we 
obtain: 
Dar Ses a 
Went e-o oe 22” (Gu2 +200,u,) 
. R2 R3 2 “2 
SP S222 2 Rt 4 
Af Re AZ Jae 
+ £(9-42°8) v2 + a =e) ue 0 (Al. 9) 
Z R2 2 6 2 a 
From hxmel = 0, using (25252. 2224) Seng (722i oe 
i 
ny AL AW, ; pees ; 
with x = Uys having A, np B: Cor eee by differentiating with 


respect to Tt (Al.5), (Al.2), (Al1.3), (Al.4) respectively and 
substituting the results obtained into (2.2.2.2.24) and 


ee. 2.2.41), we arrive ac: 














162° 2k) Eee SVC DAS | c0i a 
Ra F Tr Rs D-C A Re ae -Q - eae = Re (Cu, aE 2Cuou.) oe: 
16Z°/R* AZ ie ol 8Z°/R? : 
+ z te) = 42 * 33) u2 + _ 8275 C) ut = 0 (CAeiies@)) 
6 R2 J. 24 °R3 1. ; 
From namg]y = 0a WSC mlycreZ ace eciO) ccc. 2a 272.4 1) 
with x = Uos and making the necessary substitution, we 


@mtain: 
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loo as alee Ss F + ne pee b =f ag ©) 
gee tS > 0 Fe a 6 
1677 Ae ie 
82 5 c= 22" éu2 DOE 
i R? R Z oy: me 
87°/R° Z 
Tk yok 42, IL eee ho 
+ EQ - B3 24 ‘R: Pei = 8 
(Ae ee) 
Let: 
4278 5 2 i (Al. 12) 
Pa 1 
a7 8 os 
so C= Dy (A1l.13) 
162% 82%s os 
eC iene GC = Py LS Abel) 
Then equations (A1.8)-(Al.11) above become respectively: 
E D 2 ) .. 
IL alg Z as 2 Hie 3 
= eu, - ELS + SD LUT = 0 
2 
4Z Paya, 1-2 ae geet a 
R 
7 E Dowie ie rag ole +2Cu, ie) 
i ale ae 2 1:: aT 
pee = Raa ee eee ee = 
a 472 2 2 Z 6 lee 
sy /R 
R 
F E D 22° (Gy24 ret u.) 
1 + 1 R? Z 22 Ol Spe a 
16Z4 472 /R2 Ii. 972 /R? Dy One ee lie al 
R Marg 
: : ERY: 
Fy E D,- : ) 
C+126-36 wel ae: Ee ie ce 
i ae H(n-Gta75) Uy BS) Bao) bon new 





Late 


Pliminatronwof D, Pome A ee Sp) eee ici flr) 2 Laie: 


F, i oe 
GET RZ = = 5 u By (Allee aon) 
SUbSti Lucmumar enema pres si6n on F,/[427R] into 
(lle) and multiplying beta Sides Of tne Lesuiting equaeilen 
mye 2, we Obtain: 


of “ 97° maeene 27 = 2 
uU5D, eat) Ey ca rR? (1- “Us) G= pea [ (Cu a5 a 


2 - i. 
a (C-20) PE, 20) 


Lat No 


2Cu,u,) Ut 4h |} Us +. 


Elimination of G from (Al.20) and (Al.16) yields: 


oe 7 _ 272 . 
u,(2-u,)D, - (2-u,)E, = ora (Cus eeu nu. ) (2-u 5), 
2 72% oo a 
+ Zaye (C - 20 + 12h) (Al 2s) 


Solving for D, and E, from (Al.15) and (Al.21) above, 


we obtain: 


‘i C=06s12h 


DEL oe eens 22 
Dy = Az (Cu5zt2Cusus) =" 3 R? (a5-u (2=(a) Fu) 1 CXS 2D) 
‘ 2 2? C-204+12h 
E, == =7 - i) (Adee) 
iL Be 1s (U,-Uy, [2-(u,+u,)] “1 
From (Al.19), we have: 
Fy ee Hig _ QOS re: (Al. 24) 
A ie rei eee +U Sy 1 : 


aa ELOM “(A12.16) 3 we weve: 





j3) 0) 


AGe: 7 972 ee 47? ee ee 
ue a ae —(Cus+2Cuou,)u, Sr h + Ey ~ u,D, 


With D, and Ey from. (Al. 22). end alec ss: We have: 


472s 42°? Z> C-20+12h 
2-- 


2 
ps a pen ro 3 R7 (a. Tail) (Aol 2 





Thus far, we have expressed all the unknowns in terms 
of C, Uy and us. To determine these 3 unknowns we have the 
metlacion (Al.6), and the 2 7condi tions, of continul Gy vor ng at 


eee ana x =! 








i Zz 
PIS » GCM 202 «ee 2 6S 9) echleme ) nig E) Dy + F,. defined 
meinen 2)—{Al.14), and with £33 Eon (Ae ee) ea ea ave 
from sa er 
It 
Dy 22° (ey2 +2Cu.u ) E D 
1 R Z Lo, nN — i we 
Bi 7 Re 22Z/R i 42/R I 
,i2Z Ta, SZ ee 
3 ee O- E ))u uy + 1D Rid, 5C) uy, = (0 


With D, and Ey From (Al.22) and (Al.23) and grouping 


Similar terms, we have: 


© Bye Ss 


—~ - i LN Le eee a ie eon OL oa 
Wag" O)) 2 kay tMold L378 ee ee re Be 
zZ oa ee - ny 9 AG 7, 2 rh 
+5 (Cu, + 2Cu,u.) [ > (Ue 6 pry ] 
L Zz eee —_ oan Z ee aa 
7 os ee Q) uy tm ee &) = © 


Using the shallow shell approximation, we have: 





Pasi 


oa 7 Pr 7 72 - - ws > ome > 
e 26.4: 12h = ae u,) [2 (u,tu,) ] [6 (Cu5+2Cu,u,) uj 
= 4(3C-9) uf i 5Cut | CnieeoG) 


Then: peu obiwedane 


; 2 i. 
{ c[{l + 0(S2) } = 20 - 12h 


Al 
= 


Or, approximately: 
C = 2(0-6h) (Al.27a) 


With Q defined in (2.2.2.2.35), we have: 








C = 2(p-p.) (Al.27b) 
iiaeim: 
e = =" Z 
Cn 2[p 1 e §) pot] (A loan) 
‘Swing sil omer Uv Po mmerese @ptmg.e a = 0, we have, after a first. 
reduction: > 9 
BD, - ee Gu2e2@hiu.) aE p 
it R 2 Oe? 1. he 
et Uy — us 
SIZ / Pe 2Z2/R 4Z/R 
ii Z ae : oe a ia eo 
+ 5 B(3C @) E,)u, + 15 5 (D5 5C) u, = 0 


Then with D, and E. from (AL 22 cand (24.23) 3 wernage: 


Je if 
at aT +28 mae er ws4, 
ib ee ; 
- 2 (5(Cuz+20,u,) 4 a a 
+r = = [3C-O - = z, Taye t= (a, FIT u,] us 
+ Ts 2 22 (eu +208, u,) a CRN SCAR - 5C]us =0 


Using the shallov shell approximation, we optain: 





Zoe 


my ey = 72 Ley) e 2 
C-20+12h = pe (Ug) ) [2-(u,tu,)] [6 (Cu5+2Cuju,)u, 


~4 (36-9) u5+5Cus] (Al,29) 


By comparing (Al.29) with (Al.26), we obtain a 


relation between C,u, and u, which is: 


|) 2 
(u,-u,) {6 (Cu5+2ch,u, ) (a, tu, ) - 4(3C- eye pt Ustu, *) 
+5C(u?tu2) (u,tu,) } = (Adss3 0) 


where C and C have been determined in (Ful. 2 (is eect eo 


Another relation between C,u, u has been found in (Al.6). 


2 
Then (A1.30) and (Al.6) arethe system of equations to 
determine Uy and Uys C and c being known. 
From (Al.30) if we take: 
WU, = Uy =u 
jmem £rom (Al.6) we will find: 
15 Po 


apa kA 


Gece p.t ~ 1.2ht 





p, (1-6 a ese 


But we know that this solution is not kinematically 


- admissable, therefore, we must have 


ee) 


6 (Cu,+2Cu,u, ) (uy tus) - 4(3C-Q) (uj +u5) 


bpp! ane? 
" 9 @Q 

ee ea la = (0 (Al.31) 
meet Al. 31) CONSELEVCES Wien (ALO) eene ss hahaa of 


equations to determine Uy and u é iC and 0 being known. 


De 





Asa .B) 


APPENDIX A2 
MEDIUM PRESSURE 
6h < PoP. See 
Determination of the Unknowns: 
2 3 Be 4 
: 47, : eee 47: Tai 
Ae TROT \) Soe ee ee ee ee 


a ar, 








As has been mentioned, to determine these 13 unknowns we 


Mave 14 conditions of Continuity Us él. designate the 
jump of the function f£f at x = u, then these conditions are: 
AE x =a 
S 
ae = 0 yv], = 0 
O O 
OP Uy and x = Us 
dha | = 0 hs = 0 5; | = 0 
w) u, od us u, 
Be = e _ : ® = 
wy ee = 0 ode = vane ae 
1 al 1 
(ao ened) 


From (2.2.2.3.4) anew 272225 emnarre, 


Y wl = 00> eG eee he 220 <0) LA) 
O 


ErOM (2.262 cule Wena ger 


is 2 
gs 2° pL eee 
ere a ae aoe Oe ie ane a | ae Cig t FG ug dt 0 


(A2.2) 





254. 


Integrating by parts we have: 








T 
Weis 2 5 oly SS TI 5S 
O 
1f 
AZ2 . : 
= | (“ag ©, + Coudu_u dr 
O 
Or with (A2.1) above: 
i 2 
Zhi se 2 1 2 3 : = DU e 5 L a 3 
| (az Cius + 3 C ous) dt = “Re Cus se 5 Cus (A2.3) 
O 


Waeth this result, (42.2) becomes: 


as e 2 _ Dip = e Es e 9 
4 5 Cus (“Ez Cu. + 3 Cus) 


Za @ 
a 620) 


Wl 
© 


O 


thus; the Condi e10n Of sone rma, son vv 1S a conse- 
quence of that of yw and from these 2 conditions, we have 


only one equation: (A2.1). 


From (292.2.3.9)) and M2224 2472 260 ee ome 


ails 


yw'] =07C =- B (A2.4) 


4 ®) 
Promemtcs.2 <2. 320) tami: 22a es eee cies eee) 


above: 





vw], =") =o C. =C (A2.5) 


1 








FLOM (252.2 eo eeamcln ee ee ane 1th (Al) ee 
v] ~g +. 24 22% 4 ,»t@y-t (22° 6 
Y u, R R* i cece uy ee eile 
e 27, « 
ls 2 Be PMA : as 
Sola) Un) ae + =u, (C - = Bu,) (Ago) 
3 
Ore:. 36 22,/R R i R L 
POM Ut 2Ze2e2e2 oO) anole eo ene oe 
yw") =0>+- es a CA Zier 
u R 
2 
Slew Chime Chil Guerosm lt: 
yw] , =e Gh = onc (2c) 
ye, 
Pieeom) (202,.2.2.37)s anc m2 ee en ete or 
| Sea ee 
vv), = 0>- —— u. (Cc - = Bu.) 
2 22/R a R 
_ 24 3 
== Ku (uo) (AQ) 
As bef ill first x, 223, @ , 4278 
ore, we Wl 1rst express A, Ree oy: R2 ey! K 
Br terns Of C, Use Uys U. WSLS) SCI SNE Sys! (OAR ND) 5) 
and (A2.7)-(A2.9). 
KEOm (AZ. Sm 
K = -C U2 10) 
From (A2.7) and using the result above: 
27, = e 
2 po = (AZ a) 





IANS E 


From (A2.4) and using the result above: 





ees iG (A2.12) 
O 
From(A2.5 
Ae ss ‘ 
pe yp =O 7 QD: (A2.13) 


Pom {A2.9) and wsang (AZ 710) and {A2. 11) ,above,, we have: 


A+ 228 -_ 
27,/R | 
ome, with (A2.11): 
ee Clink 22° us) (A2.14b) 
R2 “2 | 


pmi@eamwe are left wath two Unused relations (A230 ang (42.6) 


and which will be used later in the determination of C, u 








ips 
Uy, U- 
Pog@me(2.2.2.3.7) and (2.2.2.2.40),we have: 
2 
Ze 2 ee li 3 3 
= > = ae a = 
Ps Oy ar ol) ee ey 
2 ee See 
aie E-O ae ane D _ 9 
= Re 6 a ee ae eo 
AN, 122 4Z2/R2 I Rage lL 
deers hhc 6 3 
2 a a: ae 
g bee 2 eT 
22: AJ . nae 
Having 3B, Cy, —pxC, and A by differentiating RO ie 
(A2.13) and (A2.14b) respectively, we have: 
42" 3 O 82° pee 22° (Cu2 + 2Cu.u 
R _ —R _R Z LOD, n 
aaa | 472/R2 eal 
Ae ieee C-O 1 aah 
= i 2 n ae ee. 13 Sale=s 2 = - 3 = 
= (6 Ree uy i = | ee 8, C) uy as 5 US BCU. 0 








Zo) 


2G AIBIER BS eae © 


























meom {(2.2.2.2.40) and 
2. Wee ee a 
228-60 A + 820 
R2 R3 
xs], = 00> ———— - > + 
2 Aira A158 47°/R* 
l Ar 472 ‘6 3 1 Q73 Py _ 6% * 3 
+ 3 tC - “pr EB) un + & (“Re D - BB) uy = 
s@ == K + 9 ve +2K uv 
2 a a 
Substituting again A, <5, K by their expressions in terms 
of. C. and U5, we have: 
ae 3 a 2 - ee 
ae E—O 2. oe ee aoe (Cus = 2Cu,u,) 
———+h-G- u. 
1 A A727 Re 
e246 - 422 gy? 4 L (82. Fee (A2.16) 
Do Pe Ge Eee aan 
Bieemy(2.2.2.3.13) and (2.2.2.2241)) wer have - 
hxm_ ] = 0-7 hu, + 2 2° C (u2 - 3u2u. + 2u2) + 
d uy alt 3: Re ee ik ©. 4h O 
es (uy - 4u3u, + 3u%) = 
27 *O ile cet O 
3 
a7, 4 oe 4a 27, we 407, 2 ae x 4 rend D 
at ee D- 5—B y Ee R 3 2 
_ oR R3 R fe ha as eS Ju, - uy 
16Z°*/R* AT, JiR SYA OI 
pe AM : LCC ae ‘i 
use seers 2) an ae PS 
Havin Za C 42° and A by differentiating (A2.11)- 
JR '§o!’ "R271 


(AZ2.13) and (A2.14b), 


respectively, we have: 




















ZG: 








eZ 8Z° See BU ee 2a eh 
Ro F + el pe R2 ja) - Ey C R2 (Cu, ae 2Cu,Uo) uy 
162°/R° 47,7 /R? : S7e/ Re 
lL. 42? ea) =. C-020 C—O 4 
+ e (OQ 52 E) uy + aq | ae D-C) uy + Uy xz US 
= ee 3 is lh _ 
ACU Uy + Zou, = 0 
CAD Ly) 


feeom, (2.2.2.2.41) and (2.2.2725465) 


mron Of a 4s, A and K in terms of C 


as before, we have: 








and making the substitu- 





and U5 wherever possible, 
Oo meals) ) 
AZ = 
Wy oo oe 
Be) 
oe 2 es 
-0+6+4 (h+—4 —GJu, 
BY, = Re 
ade 7} Q7 3+ . 
Maen ate = — 
See!) Bie 2 oe De) te 0 


defined in (Al. 12)—(Alw4 Sires eeorively, 


osha = 00> 
IGA Pe GVA 
gti ae ee 
75 ae 
P62 e 
87 3s 27a io 2 ne 
eke Ne ee 
Se Z 6 
Using Ey Dy, By 
equations (A2.15)-(A2.18) can be written as: 


Memo {A2.15): 





al 


é-8 


oe 27,2 oe D) ee 
Ey - Di ma (Cus = 2Cu5u,) 
MOR Ae J) 
- a u2 ot a 
2 il 6 


3 ) 
us + == 
1 2 ao 


Multapolying both sides by A ee and grouping similar terms, 


we have: 





Zoo 








DAES : Die 
a By ar Uj)u,D, ca R2 uj,)E, = 
7 27, 5 node dei o8 2 2+ 3 
ae lt (Cu, +. 2Cu,u,) uy +(C Q)ue 5 Cus] 
(A2.19) 
mom {(A2Z.16): | 
r 2 a : oar 
E - D. = =s5 (Cu5 + 2Cu,u,) 
1 ino ee ee 2 a ae 
4727R2 47 2/7R? i 2 
E D 
eee ee ae 
eg Ua 8 


Multiplying both sides by 42Z7/R* and grouping similar terms, 





we have: 
ca ee 20t oe AZo 
am el 3 R2 >) UD) + (1 = jae u,) BE, = a2 G= 
= - 22) eu? + 2Cu u.)u + 2h] (2220) 
R Z 2 2 : 
feo (A2.17) : 

i ee Bb, - 22% (Gu2 + 2¢4.u,) B 
ia + iol -- to RA ee i 3 + 
16Z%/R* ° 472/R2 8Z2/R2 ail Gia 

Dr et 6-55 1: 1 
+ 4 ao i €30 ee 3 as anda 
+ sq YY + x UGYy ZUG 3 Cusu, + 7 Cus 0 


Multiplying both sides by 42*/R? and grouping similar terms, 


we have: 
i 
1 AG ee 7 7 ee es IL mr 
@ Re) Ue ee ee ee 
_ Z ii@ue  2@u ie yu ee ene ae 
RZ 2 22a ol 
_ 4 Cpe Oy) 


¢e ee 4 e¢¢ e¢ 
3 €6 Q) us 3 Cusuy + Cull 





fieon (A2. 18): 





ak 1 . 
Rh ey on ee (h - 
<2. one a 2Cu,u,) ; 
8Z2/R2 =p 


Zee. 


Multiplying both sides by 42°/R* and grouping similar terms, 











we have: 
= i _ % he v4 Dees 2 2 ig 2 
(5 - € Ry Ug)ugd) + (1 - F Ry ug)uyk, + 
F 
Ae = i 
ors (1-u,)G + Taz R 
ae 7, 2 oe —_ 5 C226 
Seas Cu, cera Ce ie 
With the shallow shell approximation: 
a7 < AGie ATs 
Se oe Eee 
equations (A2.19)-(A2.22) become: 
ee 7 _ 22° ee > 
u,D, EL = az (Cus 5+2Cu,u 2), * (C Q)us 
. ee ee 47, ee - 27% 
u,D,-E, + oe G = sr ( (Cus +2Cu un) U. + 2h] 
1 es F, 72 Pe, oe eee 
han _ - n D 2 = 2 
5 1), u,E 4Z2/7R2 pzl(Cu5+2cu, 5) By + 2(C Q)u-u, 
_ 4 ee an 
3 (C Qu. ecu ay + Cus] 


VA2Z 22) 


(A2.25) 











Zid 


F 
ioe ae Aly = oa eal 7 
Ft ie weak oy (a nae 


72 “ y up ; 
RZ [ (Cu, 4. 2Cu,u,) uy + 4hu, + 


C06 


3 ] (A2.26) 


With the equations (A2.23)-(A2.26), we can solve for D E 


ee a 
24 and G in terms of C and Uo, Up, Ud- We proceed by 
SamMinaecitOn? Multiplying (AZ 223) sy Uy amen e( AZ. 25) 5 Oia 
and adding the resulting equations, we obtain: 
2F ; 
os a = 27, 2 es > 4 oe 3 De > _ ee b 
uj,E, + Tez7RE ae | kG Q)usu, + z(C Q)u- + cucu, Cul] 
From this, we have: 
Lie = on 
4Z°/R Papere se gal 
4 Zz" [ = (c-oyu2u, + cd (C-O)u? + 2 Gu2u. - Gu! ] 
R* ol 3 O Si hatere wae O 
(A222 7) 


Bwostituting this expression jor F)/(42°/R*) 


into (A2.26), we have: 


iL bi 7 “a 7 Az? _ ze iL sie 
FUDD) — UE) — ~pzil-us)G + 5u,F, 
7 7, 2 ae 2 ee 2 C-20 _ pes 2 
= Re [(cu, + 2Cu,u,) us, + 4hu, + =e (C Q)ucu, + 
4 ae 3 Des > 7 oe 4 
+ 3 (C Q)ue + BCucu, Cuo] 


Multiplying both sides by 2 and grouping similar terms, 


we obtain: 





DAS PAE 


> ee a9 87, 2 ee 
a = = ic 
u5D, (2u.,, u,)E, ez Gh u.)G 


_ 27? es 2 a @ y) G26 aa 2 
= ays [ (Cu, + 2Cuju,) U5 + 4hu, ia has ic Q)umu, a 
+ S(6-G)u? + Séu2u, - Cu (A2.28) 
om Oo Be Oro O 


Duletiplvang (A2.24) by 2(l-u,) and adding the resulting 


equation to (A2.28) above to eliminate cr we have: 


ee 7 _ AP 7 97,2 oe 2 e°¢ a 
U, (2 uy) D, A u,)E) = “pz [ (Cu, “| 2Cu,u.,) (2 Uy) U5 + 4h 
C90 ts aot g 2", ae 
+ ae > (C Q)usuy, 4. ZC QO) us + Bcusuy Cucl 


(EP IS) 


With (A2.23) and (A2.29) above, we can solve for D, and E, 


To eliminate E we multiply (A2.23) by —(2-u,) and add the 


nice 


result to (A2.29). The result, after some operations, 1s: 


os = 27,2 os > e16 Di eee 
D, - yz (Cuzt2Cusu5) “a (u,-u,) [2-(u,tu,) J * 


@-202 12h 7 ee es 4 ee 3 2 Gs 
x [ —s ZAC Q) ue z (2c Q) us Cus] 
AZ op 
From (A2.23), we have: 
227 — 
- ay _ Pacer 
a R 1 is OR LAY & 2(C-6) u2 + 
I hee a, PSs is 3 O 
ig | (ae. 
Le a 27,2 See 
+ 3 (2c Q)ue Cus “RZ [ (Cc Q)ue 3c) 


COR Se) 





20S: 


otha E, Given by (AZ. S13 yexerness One 222) tem a 
becomes: 
ah he ee ee 
.i Sl ee ee ee 
AZ“/R (u,-u,) [2-(a,+u,) J 2 O 


4 Ee 3 72 4 sos 3 
ar z(2C Q)u- Cul] te RZ L=(C Q)u- Cun] 


(AZ oSZ) 


To obtain G, we'can use equation (A2.24), with 7 and Dy 


from (A2.31) and (A2.30) respectively: 


472+ 47? 222/R? C-20+12h og 
"REG = “REP a= {a,+0,)) ‘ae cee eMS* 
4 afwe: 3 oy _ 27,2 ae 2 Des > 
+ q(2C Q)ue Cua] pr [(C Q)ue 3Cus] 
(A2.33a) 
aaen : 
ey ee > | 1 C-20+12h _ 
a ee ee EC 
- 2 (C-O) u2 + 4 96-6) u3 - Cu®] (A2.33b) 
O 3 O O ° 


We observe that when us = 0, expressions (A2.30)-(A2.32) and 
(A2.33a) reduce to expressions (Al.22)-(Al1l.25) of the medium 


pressure case. 


So far we have expressed all the unknowns in terms of e; a 


Uy, Up. To determine these 4 remaining unknowns, we have 


LreWwemmonc (AZ (42.6) and the 2 Gondreions of Ccontinulty 





264. 


of ns Bee uy and x = Us. BOM Ge 32 ee ee oe. Mery ing 


A and ““n by differentiating (A2.14b) and (A2.11) respec- 


tively and using 23 and E, defined in (Al.12) ,(Al.14), we 


Seocain: 
b, - 220 (Gu2 + 26h,u,) 
1. ket eee ace 
d Uy BWA fie 
D WS 2) E 





 w 
I 
(nN 
NOTQ: 
7 
= 
ll 
© 





With on anda E} from {A2.30) and (A2.31) respectively, and 
using the shallow shell approximation, this equation 


becomes: 


____ i eee pen 20e ten ~ 2(C-0) u2 + § (90-6) u3 — Cu® | 
ey laa) | 2 i . 
Se ee . 





de hdd ” e 8 2 = BE20 3 5C fh pee > 
[= (Cu. “L 2Cu,u,) uy 3 uy -- Ts Uy + (C Q)ucu, 


AS 


7 ie e¢ 3 
pemtlarly, at x = U5, we have: 


xn), ane 7 af p—2ge tah _ 
2 —lu,-u,) [2- (a, tu,)] 


7 oe > 4 een 3 _ ee 4 n a aa oe > ee > a 
2{C Q)us + 3 (2c Q)u, Cua] =F [= (Cu, + 2Cu,u,) us 
- 3C-9 \: f 5C aC“) uu. — 26a ow ] (A2.35) 
3) 2 12. ae Oo °-2 Cio ° 





AOOe. 


Comparing this equality with (A2.34), we obtain: 





= (Cuz + 2cu,u,) (ug-us) : 36-8 (y3 3) + SC (yh —u) 
+ (C-Q) u2 (u,-u,) = 26 2 (u,-uy) 0 
er: 
(u,-U,) = (Cu? + 2Cu,u,) (u,+u,) = 370 (y24u,u,tu2) 


5c a oy Xe ; 7 
+ 7a (4, tus) (uytus) se WUE Q) us Be | = 0 


From the result of the medium pressure case, we have found 


that Uy 7 Uo and therefore Uy, “Uy of O. .Thus a relation 


between Uo, Ups Us and C is, after some reductions: 
6 (Cus + 2Cu,u 5) (a, +u >) _ 4(3Cc- 0) (u2 eee : 
+ 5C(uy tuy) (ujtus) + WA e= Q)u2 ~ gCu; = 0 (A2.36) 


Comparing this relation with (Al.31) we see that the former 


will reduce to the latter when QO. 
Also, from (A2.34), we have: 


C = 20 do _ nee 2 4 2 3 oy 
R= Ro = 2(C Q)u- z (2c Que as Cu, 





3 
MS 2 Vs 2 3C-0 3 5C \ 
+ ar [= (éu2 r 2Cu,u,) uy 3 uy 4. ia 
4 
+ (C-O)u2u, - ent ] (A2. 37) 
Ok S} ell ; 


or, using the shallow shell approximation: 








DEES 
C = 2(0-6h) + 6 (C-Q) ue zs 4(26-0) u2 et 3Cut (A2. 38) 


mom (A2.1),wnhich is: 


i : 7 
me S| Sots Ss ¢ 
ee Be 
and with pz Cy ieee (O82 aL S| Co LPeOmeh 2.22 )e we navee 
- 9 
C= ra (A2.39) 


from which we derive: 


_ (l-u_)O + Qu 
CG = ee oe aes (A2.40) 
ia) 


We can also compute C from (A238) 


7 A (GE I6 = + 2u5)Q - 12h 
Se ———————— (A2.41a) 
1 - 6u* + 8u? - 3u’ 
O O QO - 


Ors 
2 eee ae) eZ 


(lor) be su.) (A2.41b) 


Comparing (A2.40) with (A2.41b), we obtain, after some 


reductions: 


= g a ee _ e@.hU68 
a us) (1l+u.)O Gl uo) (1+3u,) Qu, Wye A (ene 


Iba" 


d = 2) : 
ae ee uo) (1+u_)Q] 
Integrating, we obtain: 


is 2 a 
(1 us) (L+u)0 LZ2na ese 8 
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Since at tT = Q; Q = Q, we have L = 0 and: 
(l-u_)*(l4+u_) = a (A2.42a) 
O O 6 
With O defined in (2.2.2.2.47), we have: 
7 > - L2hT 
(1 us) (1+u.) Sliienens rae oe eam {(A2.42b) 


p,(i-e eas pot + 6ht 


Using a formula in the solution of an equation of 3rd 
degree which has 3 Leal, LOOES wil pierre. Cam=oe showme chat Uo 
can be expressed in the form: 


=, ae T+ 
us = 3 z COs ( 3 ) (hea eS 


where ¢ is determined by: 


$d = cos ¢ ——____ 2 | lL) tower > © 
4{p.(l-e je p.t + 6h7] 
(A2.42d) 
and 
¢ = cos”* [ ee - 1] fee ieee (0 


ee Det 6h) 


At t = 0, the right-hand side has the undeterminate form 


0/0. Using the L'Hospital's rule, we have, at Tt = Q: 


7 ; 7 12h 
(1-u_.)* (tu...) = Bob. + 6h (A2. 43) 


From this, we see that if PD. Pa 6h, then Ue = 0 and the 
Shell collapses only in 3 regimes as has been found in the 


medium pressure case. (u = u_(0) ) 
exe) O 
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Thus we have determined C, C and Us: TOMODp cain uy and 


U5, we have 2 relations: ({A2.6) and (A2.365). 


. 47,7. : 27° °° 
From (A2.6) and with rly Cor RB: De eis@ynnae oaks ee 


(A2.12), (A2.11) and (A2.14b), respectively, we have, after 


some reductions: 


voor 2 =— : 3 =— - 2 — 70 2 ° 3 ou 
By SIG Q) us, SCuy 3Cu5u, Sw AG ol wep 7(Gro) 0 
(A2.44) 


With C from (AD. 39) 7 we bave: 
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With these results, (A2.44) becomes, after simplification 


oo 
by y=, * 
O 
2s ae ok 2 _ 3: : 
3(2+u,) uy SU 3u5u, ae 0 Ones 
3(2+u_)u2 - 5u3 - uv? 
v2 ee ee ee (A2.45) 
2 , su, 
We observe also that when ws 0, (A2.44) wall reduce for- 


mally to (Al.6) of the medium pressure case. 


Thus relation (A2.6) becomes (A2.45) above and relation 


(A2.36) remains the same. 


lee rom 


CP Oy Uo being known, we can determine then Uy, Us 


(A2.45) and (A2.36) 
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With C and C known in terms of Yor U u we can express 


1 eg 


all other unknowns in terms of oer Uzr Uys and hence 


express the velocities and stresses in terms of Uo, U Uu 


ee 
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APPENDIX A3 


SAMPLE OUTPUT FOR MEDIUM PRESSURE RANGE 
Cle il SD Agoe < 6h) 


The symbols used in this output are: 


ZETA = ‘(Z/R)/h 
DPI = (P.-P.) /A 


J5ae 


Po/P. 


oo = To (initial value of g) 


TSTP = Te 

TAUF = Ty 

RMAX = maximum value of 1+n , 
eT Ae) 


RREST, REST1, REST2 are the approximations with which the 
solutions of the equations to obtain O00, TslP and PAUP 
have been made (the closer to 0 the better) 

IR, IFl, IF2 are indications of possible error in the eval- 
uations of QO, TSTP and TAUF. The value O means 
no error. 

Mis the number of points in the numerical integration of 
equation (2.2.2.2.74) 


TAU = tT, the independent variable in equation (2.2.2.2.74) 


Q=g 
ee — Uy 
ie — u, 


DQ = 1/g 





Ze 


P(U1) and P(U2) are values of P(x)/h where P(x) is 


iSO ei 


DD2 


ADL 


defined in’ (262.2425 eee ores sana = a 


a 2 
respectively. 
= by S/2 
= ea S72 


where S is the sum of the two root of the equation 
P(x) = 0 

mS Cond1 tion, (2,292.22 50 ie Wiener: tenn Cr menue meinS 

right-hand side has been transferred to the left-hand 

Side. 

PomeCOnGLt hom, «(2.23202 qo bots P=P. 7 0 and c@ndi tion 

eee eS tox Pape = oO” i fone DD ae 0, it has been 


proved that n, is admissible). 
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APPENDIX A4 


SAMPLE OUTPUT FOR MEDIUM HIGH PRESSURE RANGE 
(6h < pomp, < Ah) 


The symbols used in this output are: 


Depth Ratio = depth/thickness 
Pressure Difference = (po~P,) 7A 


Bees siren at To: — Po/P. 


TMU TE T > 

Leo = TS 

ZAUO = T 
O 


Reotl, REST2, REST3 and IF], IF2y) 28S) ave related co tre 


ae aly T% 


and have the same significations as in Appendix A3. 


Subprogramme subroutine that solves for T 


RPO = (pe °° - p.)/h 
RPl = (p.e ‘1 - p.)/h 


RPF = (poe ‘f - p/h 


TAU = T 
U0 = u. 
Ul = uy 
U2 = u. 


Pigewia  Conameron (2.2.2.3.43) 

DAM and DBM are coefficients of the x? term and x? term 
respectively of -the polynomial PM(x) inside the curly 
brackets of the m! expression defined in (2.2.2.3.46). 


p 
DPM2 = PM(1) value of PM(x) at x = 1. 





ey ce 





Zi. 


DNPL = Uy 


DNP2 = U5 


where S is the sum of the roots of theequation: 


oye 


oye 


eCoN = 0 
P(x) being the polynomial of 2nd degree inside the 
curly brackets of the =D Expression ini X22 oe oye 


PNP] = P(u,) 


PNP2 = P(u,) 
DNTL = R' (uy) 
DNT2 = R’' (1) 


where R'(x) is the derivative of R(x), whch is the 


polynomial inside the curly brackets of the n, expres- 


6 
Galea clic. i 2.5 2) CAR eel ley 


ENT i 


R(u,) 


IE Naa 


I 


Ra) 


Other symbols in the second stage of motion have the same 


Signification as in Appendix A3. 
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APPENDIX Ad 


INSCRIBING AND CIiRGU ee (ete eee ae eres Cyee Vea 
SURFACE BY THE UNCOUPLED DIAMOND YIELD SURFACE 


fie Cl rCumscrioing 
The uncoupled diamond yield surface fh = ]1 (Table 
2.1.3.3) is a convex polyhedron. To have this surface cir- 


cumscribe the exact Tresca yield surface f, = l, it is 


Zh 
necessary to expand the polyhedron until it is exterior to 
the surface Eo = 1. This will be the case if each of the 
eight planes of ee =K {(K > 1) istextente eer ean ne ero 
the surface at the point where the normals to the two sur- 
faces are the same. Because of symmetry and since the 
planes ny7n, = a IL tins my-m, = + Jntow Ey = iyare already 


exterior to the surface Ein = 1, 10 1S Sue riciaent to Comeucde: 


the two planes: n,+n, = 1 and m,+m, = i 


Using a procedure analogous to Section 3.7 of Reference 
{17], we find that the corresponding circumscribing planes 
are respectively: 


(Pl) = n,+n, oe (P2) = m, +m, = 


Pai meOlmmesomommcomeact being (1,1,0,0) =benme wand (070,171) 
mor (P2)e. 
Therefore if the surface fy = 1 is subjected to an 


CxPanslonm O@eraero, 2 1 will be exterior te the surface 





28s 


im, Inscribing 

Because of the convexity of the yield surfaces and 
Since the surface f,=1 LS 7a polynedron, to Wave ice anscribe 
the exact Tresca yield surface fi = oe Der So Uae C1 Cle 
to reguire the vertices of fy = — (00s -. <elys eo leon or 
within fin = 1. 

Furthermore, since both surfaces are symmetric with 
respect to moments and direct stresses, only 2 of the 16 
Wertices need be considered. They are: P(0,f£,0,£) and 


Som, 2: ,0). We find that P iS on Ene surrace mM. = Ws and 


with M, = Ny = f, we find: £ = mane = 0.618 and 0 is on 
= —_ = = 2 ] ° = — 

the surface m4 -M., 1 (ny no) ana wi em: ny oe no gear 

m,=£, m,=0, we find f = at em = 0.618. 
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